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Abstract 

We review the theoretical developments and conceptual advances that stemmed from the gener- 
alization of QCD to the limit of a large number of color charges, originally proposed by 't Hooft. 
Then, after introducing the gauge-invariant non-perturbative formulation of non-Abelian gauge 
theories on a spacetime lattice, we present a selection of results from recent lattice studies of 
theories with a different number of colors, and the findings obtained from their extrapolation to 
the 't Hooft limit. We conclude with a brief discussion and a summary. 
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1. Introduction 



In 1974, two seminal papers were published in the theoretical physics literature: on April 18, 
the article titled "A planar diagram theory for strong interactions", by Gerard 't Hooft, appeared 
in Nuclear Physics B [1], while, on October 15, Physical Review D published "Confinement of 
quarks" , by Kenneth Geddes Wilson [2] . The former paper studied a non-Abelian gauge theory in 
the double limit of an infinite number of color charges N and vanishing gauge coupling g, at fixed 
g'^N, showing that, in this limit, Feynman diagrams could be arranged in an expansion in powers of 
1/N, in which the dominant contributions come from a limited class of diagrams with well-defined 
topological properties, and revealing a close analogy with the string model for strong interactions. 
The latter paper, on the other hand, introduced a regularization for non-Abelian gauge theories, 
by discretizing the continuum Euclidean spacetime on a grid, and proved the confinement of color 
charges in the strong-coupling limit of this model; this formulation laid the foundations of modern 
lattice field theoryQ which provides a mathematically rigorous, gauge-invariant, non-perturbative 
definition of non-Abelian gauge theory]^ and makes it amenable to investigation by statistical field 
theory methods, including Monte Carlo simulations. 

Both papers have been highly infiuential during the course of almost four decades, but the 
developments they inspired followed somewhat different historical paths. In particular, the elegant 
mathematical simplifications of the large- limit, the appeal of its qualitative phenomenological 
predictions, a record of successful applications of large- A^ techniques in deriving analytical results 
for statistical spin systems, and the (semi-) analytical solution of the spectrum of large- A^ QCD in 
two spacetime dimensions [8], led to early expectations that this approach could soon disclose the 
key to understand the low-energy dynamics of strong interactions. These expectations, however, 
turned out to be delusive, and, although this approach led to a number of fruitful theoretical 
developments (including, for example, the discovery of a closed set of equations for physical gauge- 
invariant operators [9], and, later, the formulation of a systematic 1/A^-expansion for baryons [TO]), 
it did not provide an exact solution for QCD in four spacetime dimensions. 

On the contrary, the early results (both analytical and numerical) in lattice field theory were 
mostly limited to a domain of parameters (gauge couplings, quark masses, lattice sizes, as well as 
simulation statistics) characterized by large discretization effects and systematic and/or statistical 
uncertainties, and provided valuable, but only qualitative, information about the continuum the- 
ory. Although QCD (with its full continuum symmetries) is expected to emerge as a low-energy 
effective description for the lattice model when the lattice spacing a tends to zero, the suppression 
of artifacts of the lattice discretization in numerical simulations relies on the separation between 
the typical energy scale fj, relevant for a given physical quantity, and the intrinsic lattice cut-off 
vr/a, so Monte Carlo computations have to be performed in a parameter range for which fi ^ a~^; 
in addition, at the same time the lattice sizes should be sufficiently large, to suppress finite- volume 
effects. The combination of these requirements poses a non-trivial computation challenge, which 
is made even tougher by the technical complications that arise when fermionic fields are included. 
For these reasons, until the end of the last century, many lattice QCD simulations were still per- 
formed on relatively coarse lattices, and/or within the so-called quenched approximation (i.e., the 
effects of dynamical quark loops were neglected); the difficulties in approaching the continuum 



^ While the origin of lattice QCD can be identified with Wilson's article [5], ideas related to the regularization 
of field theories on a spacetime grid were also discussed in earlier works by Wentzel j3] , by Schiff [4j , by Wegner ^ , 
and in unpublished works by Smit and by Polyakov (see ref. [6^). 

^An alternative computational framework for defining non-Abelian gauge theory in the strong coupling regime 
is provided by the formulation in light-cone quantization — see ref. [71 and references therein. 
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and large-volume limits in simulations including the dynamical effects of sufficiently light quarks 
were summarized by Ulcawa at the international Lattice conference in 2001 [H], in a presentation 
including a famous plot of the computational costs, which looked like an unsurmountable barrier, 
and was henceforth dubbed "the Berlin wall" of lattice QCD. 

Fortunately, dramatic algorithmic and machine-power improvements during the last decade 
have radically changed this scenario. Today, large-scale dynamical lattice simulations, at physically 
realistic values of the parameters, have become the routine, and can provide theoretical predictions 
for quantities such as hadron masses from the first principles of QCD [121 [13] • At the same time, 
there has been a similar progress also for other challenging lattice field theory computations, 
including, in particular, simulations of large- gauge theories. This is particularly timely, given 
that during the last fifteen years, the interest in large-A^ gauge theories has received a further, 
tremendous boost, with the conjecture of a duality between four-dimensional gauge theories and 
string theories defined in a higher-dimensional, curved spacetime [HI [15], HE] , which is often called 
AdS/CFT, or holographic, correspondence. According to the "dictionary" relating quantities on 
the two sides of this correspondence, the strong coupling limit of the gauge theory corresponds to 
the limit in which the spacetime where the dual string theory is defined is weakly curved; moreover, 
when the number of color charges in the gauge theory becomes large, the string coupling in the 
dual string theory tends to zero — so that the string theory reduces to its classical limit. In 
fact, it is only in this limit that string theory can be studied with analytical methods: hence all 
holographic computations rely on the approximation of an infinite number of color charges for 
the gauge theory. For this reason, understanding the quantitative relevance of the large- limit 
for real-world QCD becomes a particularly important issue — one which can be reliably addressed 
in a first-principle approach via numerical computations on the lattice (at least for a large class 
of observables in non-supersymmetric theories). In addition, the level of control over systematic 
and statistical uncertainties in present lattice computations is also sufficient for reliable tests of 
some of the phenomenological predictions that have been derived using 1/N expansions and other 
(non-holographic) large- methods. 

This motivated us to write the present review article, whose purpose is twofold: on the one 
hand, we would like to summarize the results of recent lattice computations in large- gauge theo- 
ries, and to communicate them to the broader high-energy physics community, presenting them in 
a concise format, readily accessible also to people who are not lattice practitioners. In doing this, 
we point out some important fundamental aspects of lattice field theory, and its mathematical ro- 
bustness, but also some of the non-trivial technical and computational challenges, that one has to 
cope with, in lattice studies of certain physical problems. On the other hand, we would also like to 
present an overview of the theoretical aspects of large- A^ gauge theories, and to draw the attention 
of our colleagues working in lattice QCD to the many phenomenologically interesting implications 
of the large- A^ limit — with an encouragement to compare the results of their computations (not 
only those for theories with more than three colors, but also, and especially, those for "real-world 
QCD" with A^ = 3 colors) with the quantitative predictions for many physical observables, that 
have been worked out using 1/A^ expansions. 

In order to make this review as self-contained as possible, we decided to include an introduction 
to the large- A^ limit of QCD (and to its main phenomenological implications) in section [2| as 
well as an introductory overview of the lattice regularization in sec. [3] These sections are very 
pedagogical, and (each of them separately) should be suitable as a general introduction to the 
topic it covers. Then, section |4] discusses some aspects related to the property of factorization 
in large-A^ gauge theories: these include volume independence, loop equations in the continuum, 
the lattice Eguchi-Kawai model, and the interpretation of these properties in terms of large- A^ 
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"orbifold" equivalences; the results of various numerical investigations of reduced large- models 
are also critically reviewed. Finally, the main lattice results for physical quantities extrapolated 
to the large- limit are presented in section [5} while the conclusive section |6] summarizes the 
present status of this field, and outlines those that, in our view, are potentially interesting research 
directions for the future. 

The ideal target of this article encompasses a broad audience, including people working in 
various fields of high-energy physics, and with a varying degree of expertise on the topics that are 
presented. Throughout the various sections, we tried to keep the discussion at a pedagogical level, 
which should also be suitable and easily understandable for graduate or undergraduate students. 
In particular, we hope that some sections (e.g. sec. [2] and sec. [3]) could be useful on their own, 
as introductory material to the topics they cover, and could provide a concise overview for the 
readers interested in large- iV QCD or in lattice field theory, respectively. 

We conclude this introduction with some disclaimers. First of all, the implications of the large- 
limit for QCD and for QCD-like theories have been studied in literally thousands of scientific 
papers, for almost four decades: while it would be impossible to cover all of the relevant literature 
in the present article, we encourage the readers to integrate the discussion presented here with 
some of the many excellent reviews that are already available. The most significant early works 
on this subject are collected in ref. [E], and introductory lectures on the topics discussed therein 
can be found, e.g., in refs. [IHlIIHlEOlEIlEaEHlEl]. In addition, we would like to mention 
the proceedings of the "Phenomenology of Large Nc QCD" conference held at the Arizona State 
University in Tempe, Arizona, US in 2002 [2S], and the slides of the "Large N @ Swansea" 
workshop held at the University of Swansea, UK in 2009, which are available online Besides, 
there exist various review articles focused on specific topics, like, for example, the interpretation of 
the large- iV limit in terms of coherent states [27] , baryon phenomenology [28] , loop equations [29] , 
Eguchi-Kawai models [30], lattice results from simulations of Yang-Mills theories with two, three 
and four colors [31] and aspects related to the dependence on a topological ^-term [32j. Most of the 
lattice results in large- lattice gauge theories, that we discuss here, have also been reported in 
plenary presentations at various recent editions of the International Symposium on Lattice Field 
Theory [33l[3l|35]. 

Finally, although we tried to summarize the contents of the articles cited herein as clearly and 
as accurately as possible (compatibly with the tight constraints of a review article), we apologize 
with the authors of those works, which we may have presented in an unsatisfactory way. If they 
feel that we misinterpreted their results, or presented them in a misleading fashion, we encourage 
them to contact us to provide constructive comments and suggest corrections: we shall be willing 
to amend the preprint version of this review accordingly, before proceeding to journal submission. 

2. The large- A?^ limit 

Expansions around the large- A^ limit are a mathematical tool to study statistical models and 
quantum field theories characterized by invariance under a certain (local or global) group G of 
transformations of their internal degrees of freedom, whose number is related to a parameter A^. 
In contrast to the naive expectation, that the dynamics might get more and more complicated 
when A^ becomes arbitrarily large, in many cases it turns out that the opposite is true: often, 
the theory becomes analytically more tractable — or even solvable — in the limit for A^ — )■ oo, and 
corrections accounting for the finiteness of A^ can be arranged in an expansion around such limit, 
in powers of the "small" parameter 1/A^. 

In statistical spin systems, early applications of this idea date back to the seminal works by 
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Stanley and by various other authors EHl ESI SD]: in particular, this approach showed the 
deep connections between the nature of critical phenomena and the spacetime dimensionality, and 
laid a theoretical basis for mean-field computations. In addition, when applied to matrix models, 
large- techniques can be used to derive exact analytical solutions [IT] , and reveal the connection 
of these models to discretized random surfaces and to quantum gravity in two spacetime dimen- 
sions [121 1121 im SSI SSI SZl SH|- Generalizations to higher dimensions can be formulated in terms 
of tensor models, or group field theory |l9l ISDl EH |52l [53] . Although explicit computations in these 
models have been hindered by the difficulties in formulating a viable, systematic generalization 
of the expansion, significant progress has been recently achieved, with the introduction of 
colored tensor models [S^ — see the very recent ref. [55] for a discussion. 

For non-Abelian gauge theories — in particular, quantum chromodynamics (QCD) — the large- 

expansion was first studied in the 1970's by 't Hooft P, who suggested a generalization of the 
theory in which the number of color charges is the parameter A^, which is sent to infinity (in an 
appropriate way, as discussed below). The large- limit of QCD is particularly interesting, because 
in this limit the theory undergoes a number of simplifications — which, in certain cases, allow one 
to derive exact solutions with analytical or semi- analytical methods. Well-known examples can 
be found in two spacetime dimensions, with the computation of the meson spectrum [|8j, the 
identification of a third-order phase transition separating the strong- and weak-coupling regimes 
in a lattice regularization of the theory [561 \57\ 158] , and the determination of the properties of the 
spectral density associated to Wilson loops in the continuum theory [59] . 

In the physical case of four spacetime dimensions, the large-A^ limit does not make QCD 
analytically solvable; nevertheless, it discloses a considerable amount of information not captured 
by conventional perturbative computations around the weak-coupling limit. In addition, it also 
offers a simple interpretation — often based on elementary combinatorics arguments — for some 
empirical facts observed in hadronic interactions. Finally, a persistent Leitmotiv in the scientific 
discourse on the subject, has been the idea that the large-A^ limit of QCD may correspond to 
some kind of string theory. 

In this section, we first introduce the 't Hooft limit of QCD in subsection 2.1[ and discuss its 
main implications for the meson and glueball spectrum (subsec. 2.2) and for baryons (subsec. 2.3). 
Then, after presenting some implications of the large- A^ limit for the topological properties of the 



theory in subsec. 2.4, we review the expectations for the phase diagram at finite temperature 



and/or finite density in subsec. 2.5 Finally, in subsection 2.6 we briefiy discuss the role of the 



large- A^ limit in the conjectured correspondence between gauge and string theories [TH [151 [16] . 
which has been widely regarded as a powerful analytical tool to study the dynamics of strongly- 
coupled gauge theories for the last fifteen years. 

Although, due to space limits, we cannot cover these topic in the present review, we would 
like to mention that there also exist interesting studies of the implications of the large-A^ limit 
for hadron scattering amplitudes in the limit of large invariant energy s at fixed transferred 
momentum t [601 [61] and physics at the partonic level [62l[63l[6il[65l[S6l[67l[6Hl[69],for QCD 
evolution equations at high energy [70l[7Il[7a[73l[HIZ5l[Z6l[ZZl[Z8l[79l[80l[^ (see 
also the discussion in the reviews [851 ES]), for the hadronic contribution to the muon anomalous 
magnetic moment [HZl [HHl [HH], for the physics of electroweak processes like, e.g., those relevant for 
kaon physics [1011211 [121 [131 [Ml [15] , for studies of the "large- A^ Standard Model" , in which large- A^ 
QCD is combined with the usual electro-weak sector of the Standard Model and implications for 
Grand Unification are derived [961 [13 [SHI [99]. Ref. [100] proposed to use mathematical tools 
relevant for the large-A^ limit — specifically: random matrix theory — to study neutrino masses 
and mixing angles, considering the number of generations in the Standard Model as the "large" 
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parameter N. Finally, we would like to mention a very recent article discussing large-A^ gauge 
theories from a philosophy of science point of view |101j . 



2.1. The planar limit of QCD 

In the Standard Model of elementary particle physics, strong interactions are described by 
quantum chromodynamics (QCD): a non-Abelian vector gauge theory based on local invariance 
under the SU(3) color gauge group. The QCD Lagrangian density in Minkowski spacetime reads: 

rif 

£ = -^Tr {F^pF"^) {z^D^ - rrif) ^Jf, (1) 



where g is the (bare) gauge coupling, = — igA'^^{x)T"- is the gauge covariant derivative 
(with the T'^'s denoting the eight generators of the Lie algebra of SU(3), in their representation 
as traceless Hermitian matrices of size 3x3, with the normalization: Tr (T"T^) = 5"^/2), Fap = 
{i/g)[Da, Djs] is the non-Abelian field strength, the 7"'s are the Dirac matrices, while iplx) and 
ip{x) = ip'^'j^ respectively denote the spinor associated with the quark fields and its conjugate. 
Quark fields are in the fundamental representation of the gauge group, and occur in Uf different 
species ("flavors", labelled by the / subscript in the equation above), with generically different 
masses m/. 

The 't Hooft limit of QCD [T] is a generalization of the theory, in which the gauge group 
is taken to be SU(iV), and the number of color charges N is assumed to be arbitrarily large. 
Standard weak-coupling computations show that, in order for this limit to make sense at least 
perturbatively, it is necessary that at the same time the coupling g be taken to zero, holding the 
't Hooft coupling A = g'^N fixed. As for the way the numberof quark flavors nf should be scaled 
in the large- limit, 't Hooft's original proposal was to hold it fixed. Then, for example, it is easy 
to see that for A^ — )■ oo the two-loop perturbative QCD /3-function: 



„dl l_ f nN-2nA 3 1 f 3m'-13N'nf + 3nA 

^dfx ~ (47r)2 1^ 3 (47r)4 ^ 3N ) 9 + U{g 

turns into the a renormalization group equation for A with finite coefficients: 



(2) 



Note that, since the coefficient of the first term appearing on the right-hand side of eq. ([3]) is 
negative, perturbation theory predicts that the 't Hooft limit of QCD is an asymptotically free 
theory. While higher-order coefficients of the /3-function are scheme-dependent, inspection of the 
results obtained, e.g., in the minimal-subtraction scheme |102j shows that it remains true that the 
coefficients appearing on the r.h.s. of eq. (|3]) are finite in the large- A^ limit. It is then natural to 
assume that the Aqcd scale parameter of strong interactions is held fixed for N ^ oo. 

Another interesting feature of eq. (|3| is that it does not depend on n^: this is simply due to 
the fact that the number of quark degrees of freedom is 0{nfN), i.e. 0{N) in the 't Hooft limit, 
and hence sub leading with respect to the number of gluon degrees of freedom, which is 0(A^^). 
In fact, a different large- A^ limit of QCD (Veneziano limit) can be obtained, if rij is also sent to 
infinity, holding the nj/N ratio fixed |103j : however it turns out that this choice leads to generally 
more complicated computations, and, hence, has received less attention in the literature. Another 
inequivalent large- A^ limit of QCD (Corrigan-Ramond limit) is obtained, by assuming that ^/ is 
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fixed, but (some of) tlie fermions ( "larks" ) are in tlie two-index antisymmetric representation of the 
color gauge group |104[ll05j (which, for = 3, is nothing but the antifundamental representation), 
whose number of components scales like 0{N'^) at large A^. Finally, yet another type of large- 
limit, in which A grows like A^^ (with c > 0) for A^ — )■ oo, has been recently proposed in 
refs. |1U6[ llU7j . In the following, unless where otherwise stated, we shall focus on the 't Hooft 
limit. 

The properties of QCD in the 't Hooft limit can be studied in terms of so-called large-A^ 
counting rules, and are determined by the combined effects that arise from the number of colors 
becoming large, and the coupling becoming small. Since this double limit is taken at fixed A, it 
is convenient to write all Feynman rules by replacing g with ^J\/N . Furthermore, an easy way 
to keep track of the number of independent (fundamental) color indices appearing in Feynman 
diagrams is based on the so-called double-line, or ribbon graph, notation: since quarks are fields 
in the fundamental representation of the SU(A^) gauge group, a generic quark propagator of the 

form {■ip\x)tp'' (y)) , (where i and j are fundamental color indices, and we assume a suitable gauge- 
fixing, like, e.g., Feynman gauge) is proportional to a Kronecker delta 6^^ , and can be associated 
to a single oriented line. By contrast, the properties of the SU(A^)algebra generators imply that 
a propagator for gluon field components (y4^j (x)A^ ;(?/)) is proportional to {Sl6j — 6j6f /N), and, 
hence, can be denoted by a pair of oppositely-oriented lines. Here and in the following, we 
neglect the 6j6f /N term appearing in the expression of the gluon propagator; strictly speaking, 
this amounts to replacing the SU(A^) gauge group with U(A^), and induces a subleading, 0(A^"^) 
relative correction to the results derived in the large- A^ limit 

With this notation, it is easy to see that, in the 't Hooft limit, the amplitudes for physical 
processes are dominated by diagrams which are planar in index space, and which do not contain 
dynamical quark loops. This can be clarified by the example depicted in fig. 2A which shows three 
different types of diagrams contributing to the gluon propagator at three-loop order for generic 
A^: the planar diagram on the left panel, which contains virtual gluons only, is the one with the 
largest number of independent color lines, and (if the color indices of the external gluon are not 
fixed) is proportional to g^N^ = N'^X^. Note that this diagram can be drawn on the plane (or, 
equivalently, on the surface of a two-sphere) without crossing lines, and is thus called "planar". 
By contrast, replacing an internal gluon loop with a quark loop (central panel) removes a color 
line, reducing the overall power of A^ by one, down to A^A^. Finally, the diagram on the right 
panel of fig. 2J_, which includes a line crossing and is proportional to ^f^A^^ = A'^, is suppressed by 
two powers of A^. 

Part of the mathematical simplification of QCD in the 't Hooft limit is due to the fact that 
in this limit the dynamics is only sensitive to planar diagrams (whose number grows only expo- 
nentially with the order of the diagram, in contrast to the total number of diagrams, which grows 
factorially — see ref. |109j and references therein). In spite of this simplification, the summation of 
all planar Feynman diagrams of QCD in four spacetime dimensions is a daunting task, which has 
not been possible to complete — see refs. [HI \IV2\ \11'S\ I114[ 1115] for a discussion. 

Another interesting simplification taking place in the 't Hooft limit is that virtual quark loops 
can be neglected. Note that this feature, which corresponds to the so-called "quenched approxi- 
mation" in lattice QCD, or to the "probe approximation" in the context of string theory, arises 
dynamically at infinite A^, and, in particular, does not imply any of the fundamental problems that 



^The subleading corrections due to this difference between the U(A^) and the SU(A^) gauge theories can be 
accommodated, by introducing an unphysical U(l) "phantom" field, which cancels the extra U(l) degree of freedom 
of U(7V) [108\ . 
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Figure 1: Different types of three-loop Feynman diagrams contributing to the gluon propagator; in the 't Hooft 
hmit, only the diagram on the left panel yields a non-negligible contribution. 



occur in quenched lattice QCD at finite (which include, e.g., the lack of a well-defined Hilbert 
space, unitarity violations and negative- norm states, et c). In relation to the suppression of sea 
quark effects, we note two observations, suggesting that the 't Hooft hmit may be a "good" (i.e. 
"quantitatively accurate") approximation of real-world QCD. The first is the phenomenological 
success of quark model spectroscopy: the masses of experimentally observed hadronic states can 
be described surprisingly well, in terms of just the fiavor symmetry patterns for valence quarks, 
neglecting sea-quark effects. The second is the remarkable success of quenched lattice QCD calcu- 
lations , which are performed neglecting the dynamical effects (virtual fermion loops) induced 
by the determinant of the Dirac operator — see section |3] for a discussion. 



Note that the diagrams appearing in the central and right panel of fig. |2.1| can be drawn 
without crossing lines on the surface of a punctured sphere and of a two-torus, respectively. This 
is a manifestation of a general property of the classification of Feynman diagrams in the 't Hooft 
limit: a generic amplitude A for a physical process can be expressed in a double series, in powers of 
the coupling A and in powers of where the latter expansion is of topological nature: the power 
of associated to a given diagram is related to the number of "handles" (h) and "boundaries" 
(6) of the simplest Riemann surface on which the diagram can be drawn without crossing lines 
(and with quark lines along the boundaries), and equals its Euler characteristic: 

oo oo 

A= Y.N'-'^-'Y,c,,,,r.\-. (4) 

h, b=0 n=0 

Since at large N the leading contributions in eq. ^ are 0(A^^)and correspond to planar {h = 0) 
diagrams with no quark loops (6 = 0), the 't Hooft limit is also called "planar limit". Obviously, 
in the study of quantities involving gauge-invariant fermionic bilinears, there exists at least one 
closed quark line (which is conventionally drawn as the exterior boundary of the ribbon graphs), 
hence in this case the leading contributions come from terms corresponding to b = 1, h = 0, which 
are 0{N). 

We conclude this subsection with an observation about generating functionals for connected 
and non-connected planar graphs: in general, due to some subtleties related to cyclic ordering of 
operators, exponentiating the generating functional for connected planar graphs does not yield the 
generating functional for planar graphs only. The correct definition of the generating functionals 
for planar diagrams can be formulated in terms of appropriate non-commuting sources |117] . 

2.2. Phenomenological implications of the large-N limit for mesons and glueballs 

Under the assumption that the 't Hooft limit of QCD is a confining theory (see the discussion 



in section 3|, the large- counting rules defined in subsection 2.1 lead to a number of interesting 



phenomenological implications, in particular for the lightest physical states: glueballs and mesons. 
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To see this, following the discussion in, e.g., ref. [H], it is convenient to rescale the fields 
appearing in the Lagrangian density defined by eq. ([T]), according to: 

Alix) ^ -^A^ix), ijfix) -> VN^fix) (5) 
so that £ can be written as: £ = NC, with: 

rtf 

^ = -jx (C^'^'l + E^/ (^^'^M - mj) iJf. (6) 

/=1 

Naively, one could then imagine that the quantum theory defined by the functional path integral: 

Zo = y VAVi^V^Pexp !^iN j dt d^x (7) 

reduces to its classical limit pS], obtained as the stationary point of C This, however, is not 
the case, because an "entropic" term |118j . of the same order in iV, also arises from the field 
measure, and, as a consequence, the large- limit is not equivalent to the classical limit (see also 



the discussion in subsec. 4.2) 



Connected correlation functions of gauge-invariant (local or non-local) single-trace operators 
Oa, built from gauge fields and possibly fermionic bilinears, can be studied by adding corresponding 
source terms of the form (NJaOa) to £: 

Zj = j VAV^V^Ijex^i^N j dtd^x {t[A,%lj,i>] + JaOa 

and taking appropriate functional derivatives: 

(0.(x,) . . . 0,.(xj>. = m-'j^^ . . . ^ l.>^^.„ . (9) 

As discussed in subsection |2.1[ the leading contributions to the sum of connected vacuum graphs 
in the 't Hooft limit is 0{N'^) (or 0{N), if fermionic bilinears are considered): as a consequence, n- 
point connected correlation functions like the one appearing in eq. (|9| are dominated by diagrams 
of planar gluon loops, and scale like 0(A^^~"), in the pure-glue sector, or like 0(A^^~"), for the 
case of correlation functions involving quark bilinears. 

Let Qi be a purely gluonic, gauge- invariant, Hermitian operator, with the appropriate quantum 
numbers to describe a given glueball; as a two-point connected correlation function of the form 
{QiQi)c is 0(1) in the large- iV limit, the operator Qi creates a glueball state with an amplitude of 
order one. Glueball-glueball interactions (including decays, scatterings, et c.) are then described 
by higher-order connected correlation functions {Qi . . . Qn)c, with n > 3: as they generically scale 
like 0(A^^~"), they are suppressed by at least one power of 1/N relative to the free case. 

Similarly, meson states can be described by gauge-invariant, Hermitian operators Aii, involving 
quark bilinears, with the appropriate quantum numbers. In this case, (A^jA^j)c scales propor- 
tionally to in the large- iV limit, so that the \fNM.i is an operator creating a meson state 
with amplitude 0(1). Note that, as a consequence of this, the pion decay constant /^r, defined in 
terms of the matrix element describing the overlap between the isovector axial current and a pion 
state: 

(0|A^V(x)7^75Tj?/'(x)|7r,(A;)) = -if^k^Sjiexp{-ikx) (10) 



(in which the factor in the definition of the axial current is necessary, in order to express this 
quantity in terms of the rescaled quark fields) is proportional to \fN in the large- limit. 
Meson-meson interactions are described by connected correlators of the form: 

N^'\Ml...Mn)c, (11) 

with n > 3: they scale like 0(A^^~"/^), hence they are suppressed at large N. Finally, meson- 
glueball interactions and mixing processes are described by correlators of the form {Qi\fN Aii)c 
(or of higher order), which scale (at most) like N~^^'^, and thus are also suppressed. 

To summarize, if the 't Hooft limit of QCD is a confining theory, then its low-lying spectrum 
consists of stable, non-interacting glueballs and mesons. Exotic states like tetraquarks, molecules, 
et c. are absent in the 't Hooft limit, because the leading-order contribution to their propagators 
comes from terms which correspond to the propagation of mesons or glueballs. At finite but 
"large" values of A^ (or, better: at finite but "small" values of 1/N), hadron interactions are 
suppressed by powers of 1/\/N, so that, in the 't Hooft limit, QCD is a theory of weakly-coupled 
hadrons, with interactions described by a "coupling" which vanishes for A^ — )■ oo. It is worth 
noting that, although the dependence of this "hadronic coupling" on A^ is of the form I/a/ZV, the 
parameter determining the actual numerical accuracy of large- A^ expansions may be proportional 
to a power of it, and possibly be further suppressed by numerical factors. 

Combining the two pieces of information, that at low energies the large- A^ spectrum consists of 
non-interacting, infinitely narrow hadrons, and that at high energies the theory is asymptotically 
free (see eq. ([s])), one can prove that the number of stable glueball and meson states is infinite 
by reductione ad ahsurdum: since two-point correlation functions in momentum space can be 
expressed as a linear combination of hadron propagators with "sharp" poles, proportional to 
{p^ — im?)~^ (where denotes the four- momentum) , if the number of hadrons were finite, then 
it would be impossible to reproduce the functional form expected for meson correlators, which 
involves a dependence on the logarithm of p^. 

Another interesting phenomenological implication of the large- A^ counting rules relevant in the 
't Hooft limit of QCD, is that they provide an explanation for the empirical rule due to Okubo, 
Zweig and lizuka |119[ 11201 1121] . i.e. the suppression of strong-interaction processes described by 
Feynman diagrams that can be split in two by cutting only internal gluon lines. Fig. |2] shows 
an example of this, for two different Feynman diagrams describing the propagation of a meson 
(denoted by a blob): the diagram on the l.h.s., in which the quark/antiquark lines originating from 
the initial state propagate all the way to the final state (getting dressed by two virtual gluons 
on the way), is proportional to AA^, whereas the one on the r.h.s., in which the same process 
goes through the annihilation of the initial quark and antiquark, the intermediate emission of two 
virtual gluons, and then the creation of a quark/antiquark pair that ends up in the final meson, 
is proportional to A, and therefore suppressed by one power of A^ relative to the latter diagram. 

Large-A^ counting rules also provide phenomenological information relevant for the meson 
sector — which, in particular, allows one to discuss a low-energy effective theory for the light 
mesons |122[ 11231 11241 I125j (see also ref. |126j and references therein). In particular, the very 
existence of light, pseudoscalar mesons is related, by the Nambu-Goldstone mechanism, to the 
spontaneous breakdown of chiral symmetry. In the 't Hooft limit, this phenomenon was first 
studied by Coleman and Witten in ref. [127] . 

As it is well-known, in QCD with a generic number A^ of color charges and massless quark 
fiavorsFlthe Dirac operator i'y^D^ anticommutes with the chirality operator 75, hence for n/ > 2 



^Here and in the following, we restrict our attention to > 2. Note that the Uj — 1 case, albeit of limited 
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Figure 2: Large- iV counting rules offer a simple interpretation for the Okubo-Zweig-lizuka rule: the figure shows 
an example of Feynman diagrams relevant for the propagation of a meson: the diagram on the right, in which the 
process goes through an intermediate stage where only virtual gluons appear, is suppressed by a power of N with 
respect to the one on the left. 



the Lagrangian in eq. ([T]) is classically invariant under a global, non-Abelian chiral symmetry 
U_R(n/) X Uiin-f). This symmetry describes the independent rotations of the right- and left- 
handed components ipR^L = Pr,l4' of the quark fields (and the corresponding components of the 
conjugate spinor ip), where the right- and left-projectors are defined by: 

Pk,. = ^4^. (12) 

The global Uniuf) x Ui(?2/) symmetry of the classical Lagrangian for massless quarks can be re- 
expressed in terms of vector (V) and axial (A) subgroups as: SUv(^/) x Uy(l) x SUa(^/) x Ua(1). 

At the quantum level, the Uy(l) symmetry remains exact, and corresponds to the conservation 
of the baryon number B in QCD. Also the SUv'(n/) symmetry remains exact, and, in the real 
world, is only mildly explicitly broken by the small mass differences between the light quarks; 
this symmetry manifests itself in the approximate isospin and strangeness degeneracy patterns 
observed in the light hadron spectrum. 

On the contrary, the SUyi(n/) symmetry is spontaneously broken: the QCD vacuum is char- 
acterized by a non- vanishing chiral condensate (ipip). Correspondingly, {nj — 1) light (pseudo- 
)Nambu-Goldstone bosons appear in the spectrum: the three pions (if one assumes n/ = 2 light 
quark flavors), or the three pions, together with the four kaons and the t] (if one also considers 
the strange quark as approximately massless). 

Finally, the fate of the Ua(1) symmetry is particularly interesting, as it is related to the so- 
called Ua(1) problem: were the QCD vacuum invariant under this symmetry, then the spectrum 
of light hadrons should consist of (nearly) mass-degenerate partners of opposite parity, which are 
not seen experimentally. On the other hand, if the 11^(1) symmetry was spontaneously broken, 
then there should exist a light non-strange, isoscalar, pseudoscalar meson — but such particle is 
not observed either: the lightest meson with these quantum numbers is the t]' meson, with a 
mass of 957.78(6) MeV |129j . much heavier than the pions, kaons and the t]. The solution of 
this problem is the following: the Ua(1) is explicitly broken at the quantum level, due to the 
non-invariance of the path integral measure for quark fields |130[ I13H 1132] . The corresponding 



anomaly, however, is proportional to (for details, see the discussion in subsec. 2.4), i.e. it is 
suppressed like 1/A^ in the 't Hooft limit, so that, for vanishing quark masses, the t]' also becomes 
massless in the N ^ oo limit. As a consequence, the spectrum of large- iV QCD with nf massless 
(or light) quarks features nj- light pseudoscalar mesons, as the (pseudo-)Nambu-Goldstone bosons 
associated with the spontaneously broken axial symmetries. 



phenomcnological interest (because in nature both the up and down quarks are very light), is somewhat special jl28| . 
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Since, in general, the masses of other hadrons will be separated from these by a finite gap, the 
low-energy dynamics of QCD can be modelled by chiral perturbation theory (xPT) |133l I134j — 
see also refs. \135\ \136\ I137j for reviews of the topic — , which is an effective theory for the fields 
describing the light mesons. Assuming that the theory has Uf = 3 light quark flavors, the latter 
can be conveniently packaged into the components of a U(3)-valued matrix field U{x) defined by: 

U = exp (^iV2^/F^ , (13) 

where the entries of $ are linear combinations of the meson fields (see, e.g., ref. |126] for an 
explicit expression). According to the usual construction of effective field theories, the Lagrangian 
C^ff describing the dynamics of U must be compatible with the global symmetries of QCD, but, a 
priori, it is otherwise unspecified; in addition, it can also include couplings to external Hermitian 
matrix fields La, Ra, S and P, to account for finite quark mass, electromagnetic- and weak- 
interaction effects]^ Defining dJJ = dU — iRoU + iULa and x = ^Bq^S + iP), it is possible to 
organize the most general expression for C^g according to the number of derivatives it contains 
(i.e., in powers of momentum). At the lowest order in derivatives and quark masses, the expression 
for C^ft reads: 

C^, = ^Tr (a°Wt dj^ + + ^ (14) 

and involves two low-energy constants: F can be interpreted as the meson decay constant]^ while 
-Bo is proportional to the quark condensate |138] . Note that ~ = 0{N), so £<,ff is 0{N), 
hence for — )■ oo — see also refs. |122[ 1123] I124[ [125] — the tree-level approximation for the effective 
theory becomes exact Moreover, since the definition of U in terms of the meson fields in eq. (13) 



also contains a 1/F factor, the expansion of U generates meson interactions which are suppressed 
by powers of A^~^/^ in the 't Hooft limit. Hence, large- A^ xPT can be interpreted as a consistent 



combined expansion in powers of 1/A^, of the (squared) momentum, and of the quark masses |139 

m\m\. 

The expression of C^g at the next-to-leading order involves ten new low-energy constants 
their numerical values, fixed using experimental information on various decays |136t 1137] , agree 
with the expectations from large- A^ counting rules — which predict some of them to be 0{N), 
and some other to be 0(1). Quantitative large-A^ predictions, however, can also be obtained, by 
working in the so-called single-resonance approximation |144j . i.e. by neglecting the contributions 
of resonance states above the nonet of lowest-lying mesons with a given set of quantum numbers, 
and are in very good agreement with the values determined from experimental input (see, e.g., 
table 1 in ref. [126J). In addition, the large- A^ limit also allows one to get a better matching 
between the perturbative QCD predictions at high energies, and the xPT predictions at low 
energies |145j ,^ Further details about the applications of the large- A^ limit in xPT can be found 



^In particular, the quark masses are encoded in the diagonal entries of S. 

^Note that, with the conventions of eq. (131 and eq. (14), the phenomenological value of F (or, more precisely. 



of /^) is about 92.4 MeV. 

^Note that — in contrast to our discussion of the functional path integral in eq. ([t]) — in the present case there 
are no subtleties related to the large-A^ scaling of the functional measure, because the degrees of freedom of the 
effective theory are hadrons, whose number does not grow with N in the 't Hooft limit. 

*In principle, at this order one could also include a contribution reproducing the effect of QCD fermion anomalies: 
the effect of non-Abelian anomalies can be accounted for by a Wess-Zumino-Witten term [1421 1143j . with no free 
parameters, whereas the axial U(l) anomaly can be mimicked by a 6'-term |139j . 

large-iV effective model for quarks, gluons and pions, which could also provide a quantitatively accurate 
description of hadronic phenomena at moderate and intermediate energies, has been recently proposed in ref. |146j . 
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in ref. [139] . Other implications of the large- iV limit, which are relevant for mesons, have been 
derived for electroweak processes at low energies |147l 11481 11491 1150[ I15H I152j . There are also a 
number of works (e.g. \153\ 11541 1155j and the works reviewed in ref. |126] ) discussing the sense 
in which large-A^ QCD provides a theoretical justification for computations of various physical 
quantities in the context of the vector meson dominance model [156j . A very recent example is 
given by the study of generalized form factors carried out in ref. |157] (to which we refer, for 
further references on this topic). 

Finally, in the literature there exists a very large number of articles studying the properties of 
mesons in QCD-like models in the large- limit, using holographic methods. Since it would not be 
possible to discuss all of these works here, we refer the interested reader to the review [158 J. How- 
ever, we shall present a brief discussion of the role of the large- limit in holographic computations 
in subsection 12.61 

2. 3. Phenomenological implications of the large-N limit for baryons 

The basic features of baryons in the 't Hooft limit of QCD have been known for more than 
three decades |159] . Baryons are defined as completely antisymmetric, color-singlet states made 
of A^ quarks, thus, in contrast to mesons, the Feynman diagrams that represent them (not just the 
associated combinatorial factors) depend on A^. Another difference with respect to mesons is that 
baryon masses scale as 0{N): the mass of a baryon receives contributions from the constituent 
quark masses, from their kinetic energies (both of which are proportional to A^) and from the 
total potential energy associated to quark-quark interactions (for each pair of quarks, the inter- 
action strength is proportional to g'^, i.e. to A^~^, while the number of distinct quark pairs scales 
proportionally to A^^, so that their product is, again, 0{N)). 

To get a quantitative picture of the ground state baryon energy in the 't Hooft limit, one can 
study the large- A^ baryon in a Hamiltonian approach, and resort to an approximate Hartree-Fock- 
like solution (or to a relativistic generalization thereof): this is based on the idea that the total 
potential felt by each quark can be interpreted as a background, mean-field-like potential resulting 
from the interaction with all other quarks in the baryon. Although this approach does not lead 
to an exact analytical solution for the large- A^ baryon spectrum, it reveals that, while the mass of 
a baryon is 0{N), its "form factor" (i.e., its charge distribution) is independent of the number of 
colors in the 't Hooft limit. 

Arguments based on large- A^ counting rules also show that the amplitude for baryon-baryon 
(or for baryon-antibaryon) scattering scales as 0{N), i.e. is of the same order as the baryon mass. 
Baryon-meson scattering, on the contrary, turns out to be 0(1), so this physical process does not 
affect the motion of the baryon, but only of the meson (whose mass is also 0(1)). 

Finally, there exist physical processes involving baryons, which are missed at all orders in the 
large- A^ expansion. One example is baryon-antibaryon pair production from high-energy collisions 
of lepton-antilepton pairs. This process proceeds via an intermediate virtual photon, which then 
creates a quark-antiquark pair. In order to form the baryon and the antibaryon in the final 
state, however, this process must take place A^ times, and the probability for this to happen is 
exponentially suppressed as . This dependence on A^ is never captured at any order in the 
expansion in powers of A^~^, hence the corresponding physical process is absent at all order in 
large- A^ QCD. Similarly, the processes of baryon-antibaryon annihilation into a pair of leptons, 
and lepton-baryon scattering (which are related to the previous example by crossing symmetry) 
are also exponentially suppressed with A^. 

These features reveal an intriguing analogy, suggested in ref. |159] . between the 1/A^ expansion 
of QCD, and the usual weak-coupling expansion of quantum electrodynamics (QED), with mesons 
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and baryons playing a role analogous to electrons and 't Hooft-Polyakov monopoles, respectively. 
Indeed, electrons are the physical states in QED, their interactions can be successfully studied in a 
perturbative expansion in powers of Oqed = e^/(47r), and their mass is independent of e, i.e. scales 
like 0(e°). 't Hooft-Polyakov monopoles, on the other hand, do not appear as elementary objects 
or as simple bound states in the theory (neither do they appear in any Feynman diagrams), but 
rather as particular topological solutions; their masses are inversely proportional to Oqed, and, 
for small Uqed, their structure can be studied by classical equations (in which Oqed factors out, 
yielding results for the monopole size and monopole structure, that are independent of aqEo)- 
Moreover, both monopole-monopole scattering, and the scattering of an electron off a monopole 
are processes with well-defined, non-vanishing limits for a — )■ 0. Processes like the production 
of monopole-antimonopole pairs via electron-positron annihilation are not captured at any order 
in weak-coupling expansions of QED. As discussed above, all of this bears a close qualitative 
resemblance to the description of mesons and baryons in the expansion of QCD. 

An interesting prediction for the large-A^ spectrum of baryons of different spin J was first 
suggested, in connection to the Skyrme model (see ref. |160j for a review), in ref. jl61] . and later 
discussed also in ref. |162j : the "hyperfine splitting" between the masses of these states should be 
described in terms of a rotor spectrum: 

M = ciN + C2^^^^ + 0(l/iV2). (15) 

Other important large- implications for baryon phenomenology that were discovered during 
the 1980's include the determination of the group structure for baryon multiplets |163lll64j and its 
relation to the SU(3) chiral model and to current algebra Lagrangians and to the quark model |143l 
I165[ I166[ I167[ I168[ I169[ 1170] . Large- A^ baryons can also be interpreted as chiral solitons of the 
Nambu-Jona-Lasinio model [ITl] 1172] . 

Later, it was further realized that, combining the large- A^ counting rules, that characterize the 
't Hooft limit, with unitarity conditions, it is possible to derive quantitative predictions for several 
baryonic observables, and to organize the large- A^ expansions in the baryon sector in a systematic 
way [ini 1173] 1174] 1163] 1175] I162j — see also refs. [21], |28] for a review. This possibility arises from 
the fact that a contracted SU(2nj) spin-flavor symmetry for baryons emerges in the large- A^ limit. 

A simple way to see how these predictions can be derived is the following. At leading order, 
the coupling of a meson to a baryon is obtained by insertion of a fermion bilinear onto one of the 
valence quark lines in a baryon — which can be done in A^ different ways. Since, as discussed above, 
fermion bilinears create mesons with amplitude 0(A^^/^), it follows that the effective baryon-meson 
coupling must be 0{N^^'^). The process of baryon-meson scattering, on the other hand, can be 
obtained by insertion of two fermion bilinears onto the valence quark lines of the baryon: if they 
are inserted on the same line, then the process can occur in A^ different ways; alternatively, there 
are 0{N'^) possibilities to insert them on different lines, but then, in order to transfer energy from 
the incoming meson to the outcoming one, a gluon has to be exchanged between the two baryon 
quark lines, and this involves a factor, which reduces the overall amplitude of the process by 
one power of A^. Finally, taking into account that each baryon-meson coupling is 0{N^^'^), one 
obtains that the total amplitude for baryon-meson scattering is 0(1). 

From a kinematic point of view, however, in the scattering of a meson (with mass 0(1)) off 
a baryon (whose mass is proportional to A^), the baryon can be considered as static, and the 
process can be interpreted as the absorption of the incoming meson, followed by the emission of 
the scattered meson. Since the couplings of the baryon to both the incoming and the outcoming 
meson are 0(A^^/^), the total scattering amplitude for this process should be proportional to 
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A^. This, however, would violate unitarity, and would be incompatible with the fact that the 
baryon-meson scattering amplitude is 0(1). This apparent paradox can be resolved, if there are 
cancellations among the baryon-meson scattering diagrams, and this leads to a set of large- 
consistency conditions for baryon-meson couplings |176l 11631 1175j . 

The process of pion-nucleon scattering at low energies in the chiral limit provides a concrete 
example: the absorption of a pion by a static nucleon, A/^ + vr — )■ A/^', can be described by the 
formula: _ 

./v^^^rY-^ with- (x^-^ {B'lii^Y^.T'^i^m 

gJ\l — [X with: [X )bb' = ^7 (loj 

Itt gN 

where an overall factor (together with an arbitrary, A^-independent, normalization constant 
g) has been extracted from the isovector baryon axial current, so that X*" has a finite limit for 
A^ — > 00. At large A^, the leading contributions to the baryon-pion scattering come from the 
difference of two tree-level diagrams (corresponding to the fact that the nucleon can either first 
get into its excited state A^' by absorption of the incoming pion, and then decay to the ground 
state with emission of the outgoing pion, or vice versa). Requiring the total scattering amplitude 
to be 0(1), leads to the constraint: 

N[X'^,X^'']<0{1). (17) 
By expanding the X*" operator in powers of 1/N: 



eq. (17) yields a sequence of constraints between the different terms in the expansion; in particular. 



at the first order in 1/N, it implies that: 

[K\Xi']=Q. (19) 



Eq. (19) can only be satisfied if, besides the nucleon, there exist other intermediate states, with the 
same mass. Since a pion in the p-wave state coupled to a nucleon can also produce a baryon with 
spin and isospin both equal to 3/2 (i.e. a A baryon), the large-X consistency conditions impose a 
constraint relating QT^j^fj^f to g-j^jVA- In turn, the same argument can be applied to the A + vr — )■ A + tt 
scattering process, implying the existence (at least in the large-X limit) of a baryon with the same 
mass as A^ and A and with quantum numbers J = / = 5/2, and fixing a constraint among the 



respective couplings, and so on. The conclusion is that the consistency condition eq. (19) implies 
the existence of an infinite tower of degenerate baryons with J = I = k + 1/2, for any K G N, 
and fixes their couplings to the pion (up to the overall normalization factor g). 



Eq. (19) can be interpreted as a subalgebra of a contracted sm(4) (or su{2nf), in the case of 
Uf light quark flavors) spin-flavor algebra: the Xq" operators can be seen as the large-X limit 
of the su{A) generators that mix spin and flavor (rescaled by N). Note that the emergence of a 
spin-flavor symmetry, mixing internal and Lorentz degrees of freedom, is possible because in the 
large-X limit the baryon field is infinitely heavy, and hence static. The irreducible representa- 
tions of the contracted spin-flavor algebra can be classified using standard techniques for induced 
representations |177j . 



More generally, consistency conditions of the form of eq. (17) can be used to systematically 
organize the 1/N expansion using group-theoretical methods and operator identities, and to fix 
constraints on the subleading corrections. In particular, one finds that the first corrections to 
the mass of low-spin baryons, as well as to the baryon axial couplings, are 0(X~^) relative to 
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the respective leading terms. This framework can be studied both for rij = 2 or 3 hght flavors, 
and allows one to derive predictive theoretical expectations for a number of observables, including 
axial couplings [IZS] , masses [I62l HSOl HHH |l82l [I83] , magnetic moments [1761 HH MB [IZHl 
[T86] . masses of excited states [Ml [1881 [ffl [M [Ml [Ml [Ml [Ml [M], their couplings and 
decays [Ml [Ml [Ml [Ml 12001 12011 [Ml [2021 [2031 [203, as well as those of baryons containing 
a heavy quark [20SI [2001 [203 I20H1 [20SI I2I01 [211] • In this approach it is also possible to study a 
large- generalization |212] of heavy-baryon chiral perturbation theory |213l I214j . the nucleon- 
nucleon potential and scattering [213 MB [213 MB EISl 12201 [223 [2221 12231 [221] (see also ref . [225] , 
for work on related topics), and various quantities related to the baryonic structur^^ |2261 1227] . 
Further phenomenological implications for baryons in the large- limit are discussed in refs. |2281 
fim fim fim |2S3 12221 1233] and in the review [231]. 

Other articles on related topics include refs. |2351 12361 1237] , which used large- iV consistency 
conditions to show a flaw in the theoretical prediction of the 0"*" pentaquark state |238] (various 
other exotic baryonic states were studied in refs. |239l I240[ 1241] ). and refs. |242[ 1243] . in which a 
relation between the large- approach to baryons and the quark model was pointed out. Finally, 
large- QCD expansions based on the contracted spin-flavor symmetry also show the connection 
between the large- A^ limits of the Skyrme model |161l 11651 12441 1245] and of the quark model |246l 
1247] : this formalism enables one to derive model-independent predictions, purely based on group- 
theoretical relations. 

2.4- Topological aspects in large-N QCD 

Various types of gauge field configurations characterized by non-trivial topological properties 
are believed to play a prominent role in determining the vacuum structure in QCD and in other 
confining gauge theories: center vortices |248[ I249[ I250[ 1251] , Abelian monopoles |252j and instan- 
tons |253l 12541 113U] 12551 12561 l^^j — but possibly also more exotic objects [23S]. In general, the 
analytical understanding of the way these objects affect the dynamics, however, is limited to some 
form of semi-classical approximation. 

The topological objects, for which the large- A^ limit has the most interesting implications, are 
instantons: they are solutions of the Euclidean equations of motion with finite action, which can 
be interpreted as "tunneling" events between the topologically inequivalent classical vacua of the 
theory, classified by a Pontryagin index taking values in the third homotopy group of the gauge 
group: 7r3(SU(A^)) = Z [259]. Their relation to the large- A^ limit of QCD has been investigated 
since the late 1970's [2601 EMl [2621 [263], and is still raising continued interest [261] . 



As mentioned in subsec. 2.2, in QCD with n/ massless quarks, the global Abelian axial sym- 
metry is anomalous, due to the non-invariance of the fermion measure]^ correspondingly, the 
conservation of the Ua(1) current is violated: 

dajlix) = -2nfq{x), (20) 

by a term proportional to the topological charge density q{x): 

q{x) = ^Tr (f^,(x)F'^'^(x)) , with: F^,{x) = ^e^.^^F''^ (x). (21) 



-'^^For a different type of approach to hadronic structure functions using the large- iV hmit, see refs. [551 [751 [77115^ 
and references therein. 

^^An alternative, equivalent way to see the origin of the Uyi(l) anomaly, is from the perturbative analysis of a 
one-loop triangle diagram of virtual quarks, where the anomaly arises from the ultraviolet divergences of a one-loop 
triangle diagram of virtual quarks. 

15 



Note that q{x) can be written as the divergence of the Chern current K"{x): 

q{x) = d,K^{x), where: K^{x) = l^e^^^M^ {d,Al - p^^'^A^Al) , (22) 

which is proportional to the Hodge dual of the Chern-Simons three-form. Although q{x) is given 
by the divergence of a current, its spacetime integral (the "topological charge" Q) can take non- 
vanishing, integer values in the presence of instantons. In general, the value of Q for a given gauge 
field configuration is equal to the index of the Dirac operator, i.e. to the difference between the 
number of left- and right-handed zero-modes of the Dirac operator |265t \'2bb\ 12671 12681 [2B9j|. 
In the 't Hooft limit, the right-hand-side of eq. (20), being proportional to (?^, is — at least 



naively — suppressed like 0{1/N): as a consequence, m^, should vanish like 0{1/N) for — )■ oo. 
Thus, in the 't Hooft limit one expects the meson spectrum to include nj^ — rather than (n^ — 1) — 
massless (or light) Nambu-Goldstone bosons. 

By studying the Ward-Takahashi identities for j^(a;) + 2nfK°'{x), it is possible to derive the 
Witten- Veneziano formula |26H I262j : 



mj + m;, = 2m^ + ^^^, (23) 

J n 



(see also ref. |270j for a discussion), relating the pseudoscalar meson masses, the pion decay 
constant, and the topological susceptibihty Xtopoi, which is defined by: 

Xtopo, = -i [ d^a;(0|T{g(a;)g(0)}|0). (24) 



Combining the Witten- Veneziano formula eq. (23) with the expectation that the square of the 
rj' mass vanishes proportionally to in the 't Hooft limit (and with the fact that scales like 
jyi/s N — )■ oo), it follows that the topological susceptibility should have a non- vanishing, 0(1) 
limit at large A^. This is at odds with traditional pictures of the QCD vacuum as a semiclassical gas 
of dilute instantons, according to which the probability of tunneling events between topologically 
different vacua of the theory (i.e. the "statistical weight" of instantons) should scale proportionally 
to exp(— 87r/(7^); in particular, the number density d of instantons of (small) size p is expected to 
be proportional to |260[ 1257] : 



d{p) cx exp 



(25) 



Hp) 

i.e. to be exponentially suppressed with A^ in the 't Hooft limit. A naive dilute instanton picture, 
however, is invalidated by infrared divergences; furthermore, one should note that the exponential 
suppression which holds for instantons at the cutoff scale, may not be valid for instantons of finite 
size, see, e.g., refs. |27H 1272] 1273] for a discussion. 

The non-trivial interplay between 1/A^ suppression of hadronic mass differences, the anomaly 



equation (20), and topologically non-trivial non-Abelian gauge fields (in particular instantons) 
have also been discussed in ref. |274] . Further implications for the so-called proton spin prob- 
lem [273] were discussed in refs. [27S1 EZZl 1275] . 

It is interesting to consider the issues related to the topological aspects of large- A^ QCD in the 



context of a generalization of the theory, including a topological ^-term (see, e.g., ref. j279] for a 
discussion): in the Standard Model, there exists no fundamental symmetry principle forbidding a 
term proportional to the topological charge density q{x) in the QCD Lagrangian: 

Ce = -^Tr {Fo.pF'^^) +J2^f (^^"^" - ^/) V-/ - Oq{x) (26) 

/=i 
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(so that, in particular, Xtopoi can be expressed in terms of the second derivative of the corresponding 
effective action, with respect to 6). 

Since the topological charge density is a pseudoscalar, it explicitly breaks the discrete parity 
(P) and time-reversal (T) symmetries. For massless quarks, a 6'-term in the QCD Lagrangian can 
be reabsorbed via a redefinition of the fermion fields, by a global U^(l) rotation. For massive 
quarks, however, such rotation would make the quark masses complex, leading to phenomeno- 
logical consequences: these include, in particular, a non-vanishing electric dipole moment for the 
neutron, proportional to 6. Experimental results on this quantity |280l I281j pose very stringent 
bounds on the magnitude of the 6 term, indicating that — if non-vanishing at all — it must be ex- 
tremely ("unnaturally") small: |^| < 10~^°. The smallness of 16*1 is usually called the "strong CP 
problem" [282] of the Standard Model (see also ref. [52] for a discussion). 



If the QCD Lagrangian includes a ^-term, as in eq. (26), then the quantization of Q in integer 



units implies that the vacuum energy density of the theory is periodic in 9, with period 27r. In 



addition, the large- counting rules discussed in subsec. [271] imply that the vacuum energy density 
should also be 0{N'^) at large N . Finally, the form of eq. (26) shows that, in the 't Hooft limit, 
the natural scaling variable is 9/N] note that, in particular, this implies that the dependence on 
9 in quantities like, e.g., the string tension or the mass gap, should be suppressed like l/N"^ in 
the 't Hooft limit. A possible way to reconcile these observations consists in assuming that the 
ground state energy is given by the minimum of a multi-branched function of 9/N |124j : 

E{9) = N\ny.nf(^-^^\. (27) 



0<k<N \ N 

The potential barrier between states corresponding to different k values is expected to be propor- 
tional to A^, so that tunneling effects should be suppressed in the large- A^ limit |283l 12841 1285^ 1286] : 
at large A^, different candidate ground states become stable, but not degenerate. In addition, for 
A^ ^ 1 the energies associated with the lowest few A;-values at ^ = can be approximately 
expressed as: 

EiO) + 27r2x,op.„ (28) 

hence the gap with respect to the actual ground state energy is independent of A^. Since the 
relevant scaling parameter of the theory in the large- A^ limit is expected to be 9/N, it follows 
that the dependence on 9 in quantities like, e.g., the energy density per unit area of domain walls 
separating different vacua, or the mass gap, should be suppressed like 1/A^^ in the 't Hooft limit. 

Finally, note that, for any A^, the theory is not only CP-invariant for ^ = 0, but also for 
9 = 71. In the latter case, if two degenerate, CP-violating vacua exist, CP symmetry may get 
spontaneously broken |287] (a phenomenon that cannot occur at 6' = 0, if the assumptions of the 
Vafa-Witten theorem hold |288] ): the interesting phenomenological implications of this possibility 
have been discussed in a number of works, including refs. |l25l UMl |290l |29ll MB 12921 1293| 

MM- 

2.5. Large-N theories at finite temperature and/or density 

The focus of this subsection is on QCD at finite temperature and/or net baryon number 
densitjr^in thermodynamic equilibrium (which is relevant for the description of heavy-ion collision 



^^Since the baryon number is a conserved charge in QCD, it is possible to associate with it a corresponding 
chemical potential /i, whose value determines the net density of quarks over antiquarks in a thermodynamic 
system. 
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experiments |295j ). and on the implications for its properties that can be derived in the large- 
limit— a problem which has been studied since the 1980's [2961 EHS EHHl EHS EODl IM]. 

First of all, if the large- iV limit of QCD is confining at zero and low temperatures and densities, 
then the spectrum of physical states consists of color-singlet states: mesons, glueballs, and baryons. 
Their number is 0(1) for iV — )■ oo, and the masses of mesons and glueballs are also 0(1), while 
baryon masses are 0{N). Furthermore, meson-meson interactions scale as 0{1/N), while meson- 
glueball and glueball-glueball interactions are even more suppressed, 0(1/A^^). As a consequence 
of these properties, the equilibrium thermodynamic properties of large-iV QCD in the confining 
phase are expected to scale hke 0(1), and to be described in terms of a gas of non-interacting 
mesons and glueballs. Thus, in the thermodynamic limit, the pressure p at a finite temperature 
T can be written as: 

where the summation is done over meson and glueball states, with the index s labeling a species 
of particles with mass Ms, spin Jg, isospin Ig and = (2/^ + l){2Js + 1) physical degrees of 
freedom; Ki,{z) denotes the modified Bessel function of the second kind of index u. Note that the 



sum in eq. (29) does not include baryon states: since their masses are 0{N), their contribution 
to equilibrium thermodynamic quantities is exponentially suppressed as in the large- limit. 



On the other hand, due to asymptotic freedom |302[ 1303] , in the high-temperature limit T — )■ oo 
non-Abelian gauge theories are expected to be described as a gas of free, massless quarks and 
gluons. Since the number of gluon degrees of freedom is 0(A^^), while for quarks it is 0{nfN) 
(i.e. proportional to A^ in the 't Hooft limit), in this limit the pressure scales like 0(A^^): 

(see also ref. |3U4] for a discussion on the counting of degrees of freedom at finite temperature). 
Approximating the thermodynamics of a deconfined non-Abelian gauge theory in terms of a rela- 
tivistic gas of free quarks and gluons is accurate at high enough temperatures, where the physical 
running coupling becomes negligible. At lower temperatures (but still in the deconfined phase), the 
physical coupling is finite |305j , so that corrections to the Stefan-Boltzmann limit can be expressed 
in terms of thermal weak-coupling expansions [SOSl EHZl IMl M EIB EISl EH MB ESI 
I317[ 13181 1319] . It is important to point out that, even at high temperature, and for any number 
of colors, the deconfined plasma retains some aspects of non-perturbative nature: in particular, 
the low- frequency modes of the plasma are strongly coupled at all temperatures |320l I321j| Via 
dimensional reduction |324[ I325j . the long- wavelength degrees of freedom of the deconfined phase 
of an SU(A^) gauge theory at high temperature can be mapped to a confining effective SU(A^) 
theory in three dimensions, possibly coupled to an adjoint scalar field |313[ I326[ I327| 1328] : the 



^•^The non-perturbative aspects of the deconfined plasma are closely related to the so-called "Linde prob- 
lem" |320j . affecting perturbative expansions in thermal gauge theories, namely to the breakdown of the familiar 
correspondence between expansions in powers of the coupling, and expansions in the number of loops. This effect, 
which is due to the existence of infrared divergences in finite-temperature gauge theory (see, e.g., ref. [322] ). is 
responsible for the non-trivial mathematical structure of perturbative expansions in thermal QCD, which, include, 
for example, terms proportional to odd powers and to logarithms of the coupling, as well as terms whose coefficients 
receive contributions from infinitely many Feynman diagrams, of arbitrarily complicated topology. A pedagogical 
introduction to these topics can be found in ref. [323] . 
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latter effective theory is known as "electrostatic QCD" (EQCD), and captures the physics at 
"soft" scales 0{gT), while the former goes under the name of "magnetostatic QCD" (MQCD), 
and describes phenomena at "ultrasoft" scales, parametrically 0{g'^T). 

The deconfined phase at high temperature is expected to be separated from the low-temperature, 
confining phase by a deconfinement transition (or cross-over), taking place at a finite tempera- 
ture [3211 ESDI which for the 't Hooft limit of QCD was first discussed in ref. [2SS]. In 
SU(A^) Yang-Mills theory, the standard interpretation of the deconfinement transition is in terms 
of spontaneous breakdown of the exact global Z^r symmetry associated with the center of the 
gauge group: in the Euclidean formulation of thermal field theory, bosonic fields obey periodic 
boundary conditions along the compact Euclidean time direction, and the Yang-Mills dynamics 
is invariant if gauge fields satisfy periodic boundary conditions up to a global transformation in 
the center of the gauge group The corresponding order parameter is the Wilson line along the 
compactified Euchdean time direction (or Polyakov loop) |330[ I333| I334j : 



L(x) = — Tr 



( f^'^ 
Vexp ig / dry4^(a;, r 



(31) 



(where V denotes path ordering). In the thermodynamic (i.e., infinite- volume) limit, the thermal 



Wilson line defined in eq. (31) has a vanishing expectation value in the confined phase at low 



temperature, while it acquires a non-zero expectation value at high temperature,^^ It can be 
interpreted in terms of the propagator of an infinitely massive external color probe source, implying 
that the free energy associated to an isolated color source is infinite in the confined phase, and is 
finite in the deconfined phase 
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Although in QCD with = 3 colors, the center symmetry is explicitly broken by the quarks, 
in the 't Hooft limit the theory becomes quenched; as a consequence, the finite-temperature 
deconfinement in large- iV QCD is expected to be a genuine transition ^296j In particular, 
arguments based on the large mismatch between the number of degrees of freedom in the two 
phases suggest that the large- iV transition should be of first order (see, e.g., ref. |341j and references 
therein), with a latent heat 0{N'^) — but the possibility of a second-order transition has also been 
discussed in the literature |342j . 

An analytical study of the deconfining transition in large- A^ gauge theories on compact spatial 
manifolds was presented in refs. |343[ 1344] , by reducing the Yang-Mills partition function to an 
integral over an appropriate matrix model, which could be tackled with analytical techniques in 
the large- A^ limit. The results showed the existence of a high-temperature phase with a free energy 
0(A^^), separated from the confining phase by one or two phase transitions. Analytical studies 
based on a similar approach include refs. |345[ 1346] . and ref. |347] (which considers the case of 
orientifold gauge theories). Other interesting implications of the large- A^ limit for the thermal 
properties of non-Abelian gauge theories can be found in effective matrix models |348[ 1349] . in 
quasi-particle models |350|, 1351] , and in computations in the weak-coupling limit, when the theory 
features various types of matter content [352j . On the lattice, there are analytical studies based 



^^An alternative view on the role of center symmetry and on the physical relevance of different center sectors in 
continuum Minkowski spacetime, however, was discussed in ref. 332 . 

^^Strictly speaking, this holds for a proper definition of the Polyakov loop, which requires renormalization :335j- 

"'^^Some expectations about the Polyakov loop temperature dependence in the large- limit have been recently 
discussed in refs. |336| 1337] 1338] . using arguments related to the gauge/string duality. 

^^The existence of a finite-temperature deconfinement transition can also be argued on the basis of an exponential 
growth of the density of hadronic states p{m) as a function of their mass m |339| (see also refs. |299l 1340] ). 
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on approximate renormalization group transformations |353] . on effective actions for Polyakov 
loops |354] , and on strong-coupling expansions |355t 1356] — see also ref . |357j for a review of early 
works. 

Besides deconfinement, another important phenomenon taking place in QCD at finite temper- 
ature is the restoration of chiral symmetry [358J — which, depending on the dynamics of the model, 
can occur at the same temperature as deconfinement or at a different temperature |359[ 13601 1361] . 



As discussed in subsection 2.2, in QCD with nj massless quarks, the SUa(^/) symmetry is sponta- 
neously broken by the existence of a non- vanishing quark condensate, while the Ua(1) symmetry 
is explicitly broken by an anomaly — which, however, is suppressed in the 't Hooft limit. At 
high temperatures, both these symmetries are expected to get restored (271J and, as discussed 
in refs. [362^ 1363] . the restoration of the 11^(1) symmetry may trigger the appearance of parity- 
and CP-violating effects. The qualitative features of the finite-temperature phase diagram of 
SU(A^) gauge theories (for generic) with n/ massless quark flavors have been recently discussed 
in ref. |364j . In the large- A^ limit, one can get more quantitative insight, e.g., by calculating 
the dependence of the topological susceptibility on the temperature semiclassically, in terms of 
instantons |365[ 1366] , or by considering some phenomenological models [367J . 

Asymptotic freedom implies that QCD should also deconfine under conditions of large net 
quark density |368j , where interesting novel phases may be realized |369l 13701 13711 13721 13731 13741 
1375] . The implications of the 't Hooft limit for QCD-like theories at finite temperature and finite 



quark chemical potential fi have been reviewed in ref. [3 76 J , in which (under the assumption of 



fermionic baryons — i.e. for odd values of A^) it was also suggested that, at temperatures below 
the deconfinement transition, baryons could form a dense phase of "quarkyonic matter", with 
pressure O(A^) and with interactions between baryons near the Fermi surface (see also refs. |379l 
13801 138n M 13831 l38il 13851 13861 13871 12251 1388] . for discussions on related topics). At large A^ 
(with A and Uf fixed) the deconfinement temperature is expected to be independent of /i: this 
can already be seen from the leading-order perturbative expression for the trace of the gluon self- 
energy tensor at vanishing momentum (which is gauge- invariant), which yields the square of the 
Debye mass [3901 l39T] : 



3 2A^ V 3 7r2 



AT 



2 



+ 0(1/Ar), (32) 



and holds at all orders in perturbation theory. In the confined phase at finite values of fi, the 
contribution to the free energy due to baryons is exponentially suppressed, because their mass is 
0{N). When the quark chemical potential exceeds the critical value corresponding to the mass 
of a constituent quark rriq = Mb/N, which is an 0(1) quantity, a Fermi sea of baryons is created; 
if fi is just slightly larger than m^, the baryons (including, possibly, resonances) form a dilute 
ideal gas, with pressure 0(1/A^). For larger values of /i, however, their interactions become non- 
negligible, and yield an 0{N) contribution to the thermodynamics of the system. By going to 
much larger values of /i (in particular: /i ^ Aqcd), the free energy can be reliably computed 
perturbatively, in terms of quarks; however, the degrees of freedom in a region of width of order 
Aqcd near the Fermi surface are baryonic, i.e. non-perturbative. In addition, in such quarkyonic 
phase the chiral restoration transition is expected to take place at values of /i comparable with 
rriq, implying that, for a larger chemical potential, quarkyonic matter should be confined but 



chirally symmetric, with a parity-doubled hadron spectrum, This problem was studied using 



For a similar discussion in the Veneziano limit, see refs. |377l 1378] . 

The arguments discussed in refs. [3921 1393] may not apply in the presence of a Fermi sea. In addition, at finite 
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Schwinger-Dyson equation methods in refs. \S95\ 1396] . The imphcations of the 't Hooft hmit for 
a possible color-superconducting state at large values of the chemical potential were discussed in 
refs. |397[ 1398] . There are also analytical studies of the large- iV limit of QCD at finite chemical 
potential and at weak coupling, when the theory is defined on a compact space |399l 14001 1401] , or 
in the strong-coupling regime on the lattice |402[ 1401] . 

Various properties of the QCD phase diagram (including, in particular, the coincidence or 
splitting of the chiral and deconfinement transitions |403j . p-meson condensation |404j . the possible 
existence of second-order critical endpoints |405j . or of the chiral magnetic effect |406j ) may be 
influenced by the effects of strong electromagnetic fields: it has been estimated that the magnetic 
field strength realized in heavy- ion collisions at LHC energies could be as large as 15m^ j407j. 
The properties of the QCD plasma under strong magnetic fields in the 't Hooft limit have been 
recently discussed in refs. |408[ 1409] , in which it was found that the magnetic fields decrease the 
deconfinement temperature — at least if the quarks exhibit paramagnetic behavior. 

Finally, there exists a very large number of articles discussing the properties of QCD-like 
theories at finite temperature and/or density, using holographic methods, which implicitly rely 
on the large- limit. Like for the holographic computations of mesons mentioned at the end of 



subsec. 12. 2[ it would not be possible to cover all the relevant literature in the present manuscript. 
Hence, we recommend some recent review articles on this topic |410l 14111 14121 14131 14141 1415] . and 
limit our present discussion to a brief summary of the role of the large- iV limit in the gauge-string 



correspondence, which is reviewed in subsec. 2.6 



2.6. The role of the large-N limit in the gauge-string correspondence — A brief summary 

An expansion analogous to eq. Q also arises in open string theory, where the expansion 
parameter of the topological series is the string coupling Qg. Were the "planar" and "non-planar" 
diagrams contributing to eq. ^ to be interpreted as defined in real space (rather than in the 
internal index space of the theory), then one could imagine that, while the perturbative expansion 
is organized as a series of diagrams which, at any finite order, look like "fishnets", in a fully 
non-perturbative formulation of the theory the holes may close, and the surfaces of the Feynman 
diagrams may become real surfaces, or the world-sheets spanned by strings. This observation has 
led to speculations that string theory might provide a reformulation of non-Abelian gauge theory 
since the 1970's — a possibility which, in fact, was already discussed in 't Hooft's original work [1]. 

More recently, this idea resurfaced again in the context of the holographic correspondence, 
i.e. in the conjectured duality between gauge theories and string theories in higher-dimensional, 
curved spacetimes [TH [15], HE] (see ref. |416] for a review of the early work on the topic). This is 
clearly exhibited in the most studied example of this correspondence, namely for the maximally 
supersymmetric A/" = 4 U{N) Yang-Mills theory in four spacetime dimensions, which is believed to 
be dual to type JIB superstring theory in a ten-dimensional spacetime, with a stack of coincident 
D3-branes, resulting in a 3-brane black hole as the metric background. Approaching the horizon, 
the geometry of the spacetime in which the string theory is defined tends to the product of a 
five- dimensional anti-de Sitter spacetime times a five-dimensional sphere, AdS^ x S^, with metric: 

ds' = i-dt' + dx^) + ^dr' + R'dnl (33) 

The conjecture of a gauge/string duality arose from the observation [14J that the low-energy 
dynamics of this system can be described in terms of TV = 4 super- Yang-Mills theory. In a nutshell, 

temperature the implications related to anomalies cannot be immediately translated into a constraint relating the 
deconfinement and chiral restoration temperatures [3945 . 
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the gauge/string correspondence is a form of open-closed string duality: the excitations of the 
system can either be described in terms of open strings ending on the D3 branes (which describe 
gauge interactions in a four- dimensional spacetime), or in terms of closed strings propagating in 
the bulk (associated to gravitational interactions in the higher-dimensional spacetime). It turns 
out that the parameters of the two theories — namely, the number of color charges and the 't Hooft 
coupling in the gauge theory, and the string length and coupling, and the spacetime radius R 
on the string side — are related by: 

A = ^, ^ = 47r<7,. (34) 



As a consequence of the second relation in eq. (34), the large- limit of the gauge theory corre- 



sponds to the limit in which loop effects on the string side become negligible, i.e. to a classical 
string limit. If, furthermore, the 't Hooft coupling of the gauge theory is large, then the string 
theory reduces to its classical gravity limit (which, typically, can be studied analytically). 

The gauge/gravity correspondence has a geometric interpretation, relating the symmetries in 
the two theories in a non-trivial way: 

• the 7^-symmetry in the gauge theory is SU(4), which is isomorphic to the SO (6) symmetry 
of the five-dimensional sphere 

• the conformal invariance group of the gauge theory is isomorphic to S'0(2,4), which is the 
symmetry group of AdS^. 

The meaning of the correspondence can be formulated in a mathematically more precise way 
by stating that the generating functional of connected Green's functions in the gauge theory 
can be related to the minimum of the supergravity action, with certain boundary conditions: a 
deformation of the field theory (defined in a flat spacetime of dimension D) via a source term of 
the form: 

' A^xO{x)J{x) (35) 



(where (9 is a local, gauge-invariant operator), corresponds to a dual string theory with a bulk field 
J which reduces to J (up to a factor dependent on the conformal dimension and tensor properties 
of O) on the boundary of a (D + l)-dimensional AdS space. For this reason, the gauge/string 
duahty is also called "holographic correspondence" , as it relates the description of dynamics within 
a volume of space to information encoded on its boundary |417[ 1418] — see also |419] for a review. 

The gauge/gravity correspondence has later been extended to gauge theories qualitatively 
more similar to QCD, including non-supersymmetric theories with a linearly rising confining 
potential |42U] in the so-called "top-down approach". On the other hand, a different, more 
phenomenology-oriented, "bottom-up" approach has been pursued in a series of works [4211 14221 
Wm Wm. Wm Wm. Wm Wm Wm HHOI li3Tl gSSl iS] trying to reproduce the main features of 
QCD theories by constructing an appropriate five-dimensional gravitational background for the 
dual model. 

In principle, these gauge/string models allow one to get non-perturbative information about 
the strongly coupled regime of a field theory, by studying an appropriate limit of the dual string 
model. While both the top-down and the bottom-up approaches have led to some qualitative 
and quantitative success, it is fair to say that a completely realistic holographic dual of QCD still 
remains missing. In particular, as nicely summarized in ref . |435j , the main difficulty appears to be 
the inclusion of non-negligible a' string corrections, which account for the finiteness of the physical 
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coupling in the dual gauge theory: this comes from the relations between the parameters of the two 
theories expressed in eq. (34), revealing that at fixed 't Hooft coupling, the string coupling constant 
is inversely proportional to the number of colors of the gauge theory, so that quantum corrections 
on the string side of the correspondence are suppressed if is large. Furthermore, if in this limit 
the 't Hooft coupling is large, then the string length Ig becomes negligible with respect to the 
radius of the space on which the theory is defined, and hence string theory essentially reduces to 
classical gravity (in an appropriate higher-dimensional spacetime). Since this limit is particularly 
interesting, because it can be studied analytically or semi-analytically, the correspondence is often 
called "gauge/gravity" dualityp°| In addition to the vast literature investigating the gauge/gravity 
duality with analytical methods, it is worth mentioning that, especially during the last few years, 
there has been a tremendous boost of studies of related problems using numerical tools: for an 
overview, see, e.g., ref. |437j and the links to the online talks therein. 

Finally, the large- iV limit plays a crucial role also in the integrability of the A/" = 4 super- 
symmetric Yang-Mills theory: the planar sector of this can be "solved" (meaning that scaling 
dimensions of local operators can be expressed analytically as a function of the coupling constant) 
in terms of a system of algebraic equations, which can be derived as a particular limit of the 
integral equations obtained from the thermodynamic Bethe Ansatz |438l 14391 14401 1441] — see the 
review |442] and references therein for a discussion. 



3. S\J{N) Gauge Theories on the Lattice 

For energies above 1 GeV, QCD is perturbative. Hence, a controlled expansion in as = g"^ / (47r) 
can be used to compute physical quantities. This programme is successful at explaining experimen- 
tal data obtained at particle collider facilities, and is indeed one of the key theoretical ingredients 
for the analysis of collision events observed at the LHC However, as the energy is lowered below 1 
GeV, as rapidly becomes of order one, and the perturbative expansion looses its predictive power. 
Hence, in order to derive from first principles phenomena that are typical of the non-perturbative 
domain of QCD, like chiral symmetry breaking and confinement, a different approach is needed. In 
the previous section, we have discussed the large- iV expansion as a possible systematic expansion 
to study non-perturbative QCD and the gauge-gravity dual as a way to compute observables in 
the zeroth order in that expansion (corresponding to the N ^ oo limit of the theory). In this 
section, we shall discuss a complementary approach, based on numerical simulations of the theory 
discretized on a spacetime lattice using Monte Carlo techniques. Once the continuum limit is 
taken, numerical results for QCD obtained in the lattice framework can be compared to obser- 
vations. Thanks to recent technical advances in the field of lattice gauge theory, not only do we 
have now a reasonable numerical proof that QCD is the correct theory of the strong interactions, 
but we can also trust the predictive power of lattice QCD calculations. The progress in the field 
has encouraged further development targeting a more theoretical understanding of gauge theories. 
Among possible directions, numerical studies of the 't Hooft large- limit have generated a large 
volume of quantitative results that on the one hand significantly advances our understanding of 
the theory and on the other hand can be used to inform analytical approaches. In this section, 
we lay the foundations for understanding Monte Carlo results for SU(A^) gauge theories in the 
't Hooft large- limit. While this part brings only limited benefit to the reader already familiar 
with numerical calculations on a spacetime lattice, it could provide a useful reference for those 



^°Note that, strictly speaking, this is actuaUy a misnomer [436], since the correspondence (in its "strong" form) 
is expected to hold at all values of the parameters of the string theory. 
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who are interested in using lattice results to inform or inspire analytical calculations. The expo- 
sition is pedagogical and aimed at underlining the main conceptual steps and the technicalities 
needed for a critical understanding of the numerical results and fixing the notation for later chap- 
ters rather than at providing detailed derivations of all the basic lattice results needed later on. 
The reader interested in this latter aspect is referred to the excellent textbooks on lattice gauge 
theories |443l I444[ 1445^ 14461 1447[ I448j . The works reported in the early original literature on the 
subject are by now considered classic papers in the field. Although we have mentioned the main 
original references, there are many more to which we can not and we are not making justice. Once 
again, we refer the reader to the textbooks, which contain a better account of the early original 
contributions to the field. 



This section is organized as follows. In subsec. 3.1 we discuss the discretization of the free scalar 



field and some general aspects of the recovery of the continuum limit. How to construct lattice 
variables for the gauge fields so that gauge invariance is respected in a discretized spacetime will 



be the subject of subsec. 3.2, which will also deal with subtleties connected with the discretization 



of fermionic fields. Numerical calculations will be discussed in subsec. 3.3, while in subsec. 3.4 
we shall present how continuum large- physics is extracted from numerical simulations. For the 
sake of definiteness, unless otherwise stated, in this section we only deal with four-dimensional 
theories. 

3.1. The free scalar field 
3.1.1. Path integral approach 

The lattice approach is based on the path integral quantization. For the sake of definiteness, 
let us consider a free scalar theory, described by the Lagrangian 

= ]^d''<P{x)d,cP{x) - \n?<\>{x)\ (36) 

where is a spin zero scalar field of mass m. The easiest way to quantize this theory is to use the 
canonical approach. With x = {xo,x), with xq being the temporal component and x the spatial 
component of the quadrivector x, we define the momentum n(x) as 

n(x) = (37) 

and we impose the usual equal-time (xq = yo) commutation relations 

= [n(x),n(i/)] = o. 

This approach naturally leads to the Fock space, in which the base states are multi-particle 
states labelled by the momenta of each single particle. 



If we now consider an interaction that can be written as^^ V{(f)) = af{(j)), with / e.g. a 
polynomial of degree k, for small a we can still start from the canonical quantization and treat 
the effect of the interaction perturbatively in a. It is then possible to compute systematically the 
n-point correlation function for arbitrary n 

C„(0(xi), = (38) 



^"'^We assume that V{(f)) is bounded from below. 
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as a power series in a. Since all the observables can be expressed using the C„, all physical 
processes can be accessed in this way. The basic physical assumption underlying the perturbative 
expansion is that multiparticle states of well-defined momentum are still a good approximation of 
the eigenstates of the interacting theory, with corrections that can be accounted for by a systematic 
expansion in powers of a. However, in a quantum field theory a is not constant, but depends on 
the momentum. Hence, in order for the perturbative calculation to be valid, a needs to remain 
small for the relevant range of energies. 

An alternative approach for computing the correlation functions is the path integral. In this 
formulation, they are given by 

where has to be interpreted as a classical field and the integral has to performed over all 
possible classical field configurations (this is indicated by the expression (V(j){x))). The integrand 
is weighted by the factor e*"^*^'^-', where 



S{(f)) = J d^x C{(f)) (40) 
is the action evaluated over a field configuration (f){x). The denominator 

Z = [ {V<p{x))e'^^^\ (41) 



needed to normalize the correlation functions, is what is referred to as the path integral. The 
advantage of the path integral formulation is that it can be used also for the interacting theory, 
irrespectively of the value of the coupling. 

3.1.2. The scalar field on the lattice 



At this stage, the path integral (41 ) is still a formal expression: in order to be able to use it, we 
need to give a prescription on how to perform the integration. To this purpose, we first perform 
a Wick rotation, which consists in the change of variable r = ixq. In the new variables, up to an 
overall minus sign, the metric is Euclidean. For this reason, the space of the vectors xe = {ixQ,x) 
is called Euclidean space. In Euclidean space, we write the Lagrangian as 

^^(0) = ld'^<PixE)d^,E<P{xE) + ^m^(pixEf, (42) 

where the metric is the identity (and in fact we could disregard the convention of lower and upper 
indices), and the action as 

Se{^) = j (fxE Ce^cI)). (43) 

Then, the path integral becomes 

{V(t){xE))e-^^^'^\ (44) 



Over eq. (41 ), Ze has the advantage of having the dumping factor e ^^{4)) replacing the oscillating 
factor e^'^^^Mn the integrand, which improves the convergence of the integral. Real-time correlation 

25 



functions can be obtained from Euclidean-time correlation functions by analytic continuation. 
From now on, we will work in Euclidean space, and for convenience we will drop the subscript E 
from all expressions. For instance, we will write the Euclidean path integral as 

{V(f){x)) e~^'''''\ (45) 

where it is understood that (p is defined in Euclidean space and S is the Euclidean action. 

We still need to give an operational implementation of the measure {V(j){x)). To this purpose, 
we can consider a spacetime grid of spacing a (the lattice) and define only on sites of the grid, 
i.e. on points x such that x = (nga, riia, n2a, n^a), with the rij all integer. If we also impose that 
< Hi < Ni — 1, i.e. that the lattice size in the i-th direction is aiVj, then the path integral 
measure becomes a multidimensional integral. In the remainder of this section, we focus on the 
discretization of the theory, while in the next section we shall see how the continuum theory can 
be recovered from its discretized version. 

For convenience, we rescale the field by a, obtaining the dimensionless combination (p{i) = 
a(j){ai), with i = (no, ni, n2, ns) = (no,n), in terms of which we formulate the free field scalar 
theory on the lattice. To this purpose, we replace the integral / d a; with where i runs over 

all lattice points, and the derivative of (p with the finite difference 

+ fi) - 

d^(j){i) , (46) 

(X 



where (x is the unit vector in the direction /i. The lattice action becomes 

(47) 

where m = am, and the discretized path integral reads 

-s{^)^ (48) 



i I \iJ.=0 J 



with S given in (47) and the product running over all lattice points, n— point Green functions are 
easily computed in momentum space, where it has to be considered that the momenta are cut-off 
at Pmax = vr/a (in crystallography, this corresponds to the first Brillouin zone). For instance, for 
the two-point function (the lattice propagator), we have: 

{v>{lMm)) = I l2.^. , ,^2 - (49) 



n (27r)MV sin^(p;/2) 



m 



3.1.3. Continuum limit 

The lattice path integral (48) is formally identical to the partition function of a statistical 
system with degrees of freedom and Hamiltonian S. The problem of recovering the 

infinite volume limit of the original system is then mapped into the problem of performing the 
thermodynamic limit of the associated statistical system. This mapping allows typical concepts 
of statistical mechanics to be carried over to lattice field theory. Although we do not pursue the 
analogy further, important progress has been performed exploiting this correspondence. 
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Once the infinite volume fimit has been performed, the continuum hmit can be recovered by 
taking the lattice spacing a to zero. At the classical level, this implies that the discretized action 
should reproduce the continuum action in the limit a — )■ 0. This request is easily fulfilled, since 
the lattice action has been constructed as a simple discretization of the continuum action. For 
the quantum theory, we need to systematically compute all lattice n-point functions and show 
that they converge to the corresponding continuum functions when a — )■ 0. Since a quantum field 
theory is uniquely specified by its n-point functions, this would suffice to prove that the lattice 
field theory reproduces the wanted field theory in the continuum limit. The above procedure 
can look tautological. However, there is no guarantee that the lattice theory obtained by naive 
discretization of the path integral describe the wanted field theory in the continuum limit. 

Since for the free theory the issue is trivial, let us consider the case of the interacting theory. 
Once the theory is discretized on a lattice, observables depend on the dimensionless couplings 
(m, a, . . . ), which in turn are functions of the lattice spacing. If we now take the lower mass of 
the physical spectrum, M, if a continuum limit exist, we must have 

lim M/a = M, (50) 

a->0 

where M is the dimensionless mass determined on the lattice, and a reinstates the dimensions. 
The previous equation implies 

liml/M = oo, (51) 

a— >-0 

where 1/M = ^ can be interpreted as the correlation length (in dimensionless units) of the 
statistical system associated to the regularized quantum field theory. Hence, in the language of 
Statistical Mechanics, reaching the continuum limit means finding the values of the couplings 
for which the system is critical. Physically, this means that the correlation length of the system 
is much larger than the lattice spacing a, which implies that the system looses memory of the 
discretization. In general, whether values of the coupling exist such that the system is critical 
is a dynamical problem. If it happens, it might not be immediate to identify the corresponding 
quantum field theory in the infrared, since this is determined by the value of the n-point functions 
at the critical point. For SU(A^) gauge theories, thanks to asymptotic freedom, the system is 
critical at the ultraviolet fixed point, where the theory is pcrturbative. Hence, the lattice provide 
a way to regularize the theory in the ultraviolet, and removing the ultraviolet cut-off a corresponds 
to a well-defined renormalization procedure, in which an infrared scale (e.g. the string tension in 
the pure gauge case or the mass of the p meson in the presence of dynamical quarks) is fixed, and 
all spectral quantities are determined in terms of this scale. In fact, once one has proved that 
the lattice discretized version of QCD has the right continuum limit, one could invert the logic 
and define QCD starting from the construction on the lattice. This would give a more rigorous 
way to specify the quantum field theory than canonical quantization, which has the problem of 
defining the theory starting from the pcrturbative Fock vacuum, which is very far from the real 
QCD vacuum. By the same token, using the lattice we can rigorously define the large- A?" limit d 
la 't Hooft of SU(A^) gauge theories as the theory defined by all the n— point correlators when 
this limit is taken. In this approach, the lattice could be used to prove the existence of this limit, 
even in the non-perturbative regime. 

In passing by, we notice that difi'erent prescriptions could have been used to construct the 
lattice action from the continuum one. For instance, we could request that the discrete derivative 
be defined as 

^ ■ (52) 
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The corresponding action would have differed from (47) by corrections that go to zero in the 
continuum hmit. This ambiguity in defining the lattice action can be exploited to improve at the 
quantum level the convergence of the lattice theory to the continuum one as a — )■ 0. 



3.2. Discretization of SU{N) gauge theories 
3.2.1. Gauge fields 

A good regularization of a quantum theory respects the crucial properties of the original theory, 
with the others recovered when the ultraviolet cut-off is removed. On the lattice, the property 
that we absolutely need to preserve is gauge invariance. This can be accomplished through parallel 
transports. Consider for instance continuum scalar electrodynamics (in Euclidean space). If the 
(complex) scalar field is coupled to an Abehan gauge field A^{x), the continuum derivative is 
replaced by the covariant derivative: 

D^ = + igoA^{x), (53) 
where 5^0 is the gauge coupling. Under gauge transformations defined by the function A(x) 

A^{x) ^ A^{x) - d^A{x), <j){x) ^ e-'^'^'^'^^Uix), (54) 
so that the Lagrangian 

Ci(P,A^) = D^(l){x)D^(l)*{x) + mmx)(t)*{x) + ^-F^,{x)F^,{x) , F^,{x) = d^A, - d,A^ (55) 

is invariant. It is immediate to see that a naive discretization of this Lagrangian will not preserve 
gauge invariance, the problem being that the finite difference mixes fields defined on different 
lattice points. This can be remedied as follows. We introduce the parallel transport along the link 
joining the sites i and i + fi as 

f/^(^) = e^5o/.^^'^'^-(^)'i^ ^5Q^ 
where x = ai. On the lattice, we make the replacement 

D,<f>i^) U.i^m - <^(i) _ (57) 
Gauge transformations act as usual on (f): 

m ^ e-^^(^-)0(j); ^ e''(^-)0*(j)- (58) 

For the link variables, we have 

U,{j)^e-^'^^^^U,{j)e'^^^+^\ (59) 

Using the fact that U-^{i) = U*{i — fi), it is easy to see that the terms involving in eq. (55) are 
invariant under discretized gauge transformations when using the prescription (57) for the covari- 
ant derivative. The part of the Lagrangian involving the field tensor F^^y can also be expressed in 
terms of the link variables. The simplest possibility is given by the Wilson action [2] 
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where (5 = l/g^ and 



u,Aj) = u,{j)u,{j + fi)u;{j + i))u:{j) (61) 

is the parallel transport of the gauge field around an elementary lattice square (plaquette). It is 



worth noticing that 5* defined in eq. (60 ) is manifestly gauge invariant. Expanding the exponentials 



defining the links at the leading order in a, we get 

^^tE^^'^(^')^^^(^')' (62) 



4 

which is the naive discretization of the gauge field action. 

The generalization of the above discussion to SU(A^) is immediate. Now, the link variables 
U^{j) are defined as the path-ordered exponential 

where now U^{j) is a matrix in the SU(A^) group, and the plaquette is the path-ordered product 
of links around the lattice plaquette: 

U,.U) = U,{j)U,U + fi)Ul{j + 0)Ul{j) (64) 

where as before U^^{j) = f^j_/i(j)- Under gauge transformation implemented by the SU(A^)-valued 
function Q{x), the t/^(j) transforms as 

u,ij)^nij)u,ij)n\j + fi), (65) 

which is the lattice version of the continuous gauge transformation 

A^{x) -> Q{x)A^{x)Q\x) - -Q{x)df,Q\x). (66) 

9o 



In terms of U^^{j), the Wilson action reads 



5 = /3 5^ - lReTrt/^,(j)) . 



(67) 



where ReTr indicates the real part of the trace and P = 2N / g^. Gauge invariance is guaranteed 
by the trace. The path integral reads 



\{dU,U)]e-', (68) 



where each factor dU^{j) in the measure is the Haar measure of SU(A^) associated to the link 
Ufj_{j). For our purposes, the most important property of the Haar measure is that it is uniform 



in the group. The action (67) yields the naive discretization of the SU(A^) Yang-Mills theory at 
the leading order in a as a — )■ 0, as it should. The lattice theory can also be shown to reproduce 
the correct continuum limit when the lattice spacing goes to zero. As we will see in more detail 
later, this is a non-trivial consequence of asymptotic freedom. 
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3.2.2. Fermions 

Lattice discretization of fermion fields is not immediate. To see the origin of the problem and 
the possible solutions, we start by following closely the strategy we have used for the bosonic field. 
Given a continuum field ip{x), we define the (dimensionless) discretized fermion field as 



and its discretized derivative as 



4^{i + fi) jl) 



(69) 



(70) 



Note that, differently from the bosonic field case, we have used a symmetric definition for the 
derivative. By substituting these equations in the Euclidean Dirac action, we get the lattice 
fermionic action 



V'/3(j), 



(71) 



where m = am, with m the fermion mass. In this equation, Dirac indices are expressed with 
Greek letters, while Latin letters indicate lattice points. The 7^ are the Euchdean 7 matrices, 
satisfying the anticommutation relations 



In terms of the Minkowskian 7 matrices, the Euclidean 7 matrices are given by 

7^ = 70, lk^-il\ i = l,2,3. 
It is convenient to define the Dirac operator D as 



(72) 



(73) 



(74) 



so that the fermion action can be written in a more compact form as 



(75) 



Fermionic correlation functions can be expressed in terms of the inverse Dirac operator. For 
instance, for the correlator, we have 



(4(^)^MJ)) = Dl'M- 
Using the momentum representation, we can rewrite this expression as 



(V'a(i)V'/3(i)) = 



d^P -^E^7jsin(p^)+m 



(76) 



(77) 



where p = ap and p is the continuum momentum. This expression reproduces the continuum 
Euclidean propagator of a free fermion when a — )■ 0, as it should. However, the same continuum 
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form is obtained also when at least one of the components of the momentum = n/a. Hence, 
each corner of the Brillouin zone contributes equally to the propagator, which means that the so- 



called naive discretization, given by eq. (69), eq. (70) and eq. (71 ), gives rise to 16 fermion flavors 



in the continuum limit. This is the problem of fermion doubling. A well-known no-go theorem 
due to Nielsen and Ninomiya [449 [ l45Uj implies that doubling is an unavoidable consequence if 
one requires a discretized fermion theory that preserves chiral symmetry and is ultralocal (i.e. the 
action only involves couplings between fields in a localized region of space). Hence, in order to 
avoid doubling, one has to relax the request of chirality or the request of ultralocality. The Wilson 
solution to the problem of doubling was to relax the request of chirality |451j . In his approach, an 
irrelevant operator in the limit a — )■ provides the doublers with an infinite mass in the continuum 
limit A different approach was proposed by Kogut and Susskind |455l 14561 1457] . who were able 
to reduce the number of doublers from 2^ (in D spacetime dimensions, assuming D to be even) 
to 2^/^ by spreading the four components of the Dirac spinor on the corners of the Brillouin zone 
(hence the name of "staggered fermions" for this discretization). More recently, a series of fermion 
discretizations, that avoid the doubling problem, by satisfying a modified form of chiral symmetry 
on the lattice |458t 1459] (which goes over to the usual continuum chiral symmetry for a — )■ 0), 
have been proposed: these include the domain wall formulation |460l 14611 14621 1463] . the overlap 
formulation |1MI ESSl HSEl SSH ESHl 1169], and the fixed-point formulation EH I1Z2]0 For a 
general discussion about lattice fermions and chiral symmetry, see refs. |474] 1475] I476j . We stress 
that the physics in the continuum limit is independent of the lattice discretization used. However, 
a particular formulation could be more suited for a particular problem. In lattice calculations of 
SU(A^) gauge theories with fermions, mostly the Wilson formulation has been used. The main 
motivations for this choice are the following: 

• simulations using Wilson fermions are much faster than simulations using non-ultralocal 
fermions; 

• unlike in the staggered fermion case, a generic number of flavors can be simulated in the 
Wilson approach; 



chiral symmetry can be recovered by tuning the bare quark mass to a critical value 
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The Wilson discretization starts from a modiflcation of the Dirac operator with an additive term 
that goes like the Laplacian of ijj with strength controlled by a parameter r, conventionally taken 
to be equal to 1. Explicitly, this term reads 

This yields the following Wilson Dirac operator 

DZi'j) = E [(^ - ^?)a/3 + + 7j),^ + + 4r) 6,,6^p. (79) 



^^A variant of Wilson fermions, including a "chirally twisted" mass term, was proposed in refs. |452|, 1453] 1454] . 
^•^A "chirally improved" variant of Wilson fermions was proposed in ref. [473] , 

^''Although this value has to be found as a part of the simulation, this does not create particular problems in 
practical applications. 
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The resulting two-point function is 

J-^ (27r) E^sm (p^) + (m(p)2) 

which is similar in form to the naive propagator, the difference being that now the mass depends 
on the momentum: 

m[p) = m-\ sin (81) 

(note that we have reinstated dimensionful units in this expression). When ap^ — )■ for all /x, 
m{p) ^ m, and the expected continuum result is recovered. Conversely, at the other corners of 
the Brillouin zone, m{p) diverges. Hence, in the continuum limit the unphysical doublers get an 
infinite mass, decoupling from the action. 

3.2.3. Gauge theories with fermionic matter 

We conclude this brief introduction to lattice gauge theories with the case of fermionic matter 
coupled to gauge fields. For simplicity, we consider the case in which there are Uf fermion flavors 
having the same mass m. In the presence of gauge interactions, for Wilson fermions the Dirac 
operator is given by 

D^pi'j) = E [(^ - ^^)a(s U,m+f., + {r + 7j)„^ Ulm-f.,j\ + + 4r) S,,S^p. (82) 

For the sake of clarity, wherever the context is unambiguous, from now on we suppress the Dirac 
indices and the spacetime coordinates in and in the fermion fields. In the most straightforward 
formulation, the path integral of the theory is given by 

^ = / (U^U.ij)] md^Lij)] ( n#'.O0 ) e-'-^'^-^^\ (83) 



where / is the flavor index, running from 1 to nj, and 5* is the gauge action (67). Performing 
explicitly the integration over the fermionic variables gives 

^ = / j (det Dw)"^ e-^ (84) 

in which the determinant of Dw to the power of appears. 

Among fermionic observables, we shall consider only zero-momentum correlation functions of 
isovector meson operators, which take the form 

Crr'(t) = X] (^(^'(t,f)r^'=(t,x)y (^'{0,0)T'i:\0,6)^y (85) 

where F and F' are two Euclidean Dirac matrices, (0, 0) is the conventional origin of the lattice 
and / 7^ k. Integrating the previous expression over the fermion fields yields 

Crrit) = - X] (tr (Ft(D^)-i(a;, 0)F'(D^)-i(0, x)) > , (86) 
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with = (0, 0) and x = {t, x) and tr being the trace over Dirac indices. More exphcitly, in the 
path integral formulation this expression reads 

Crr'it) = -^J (^dU.ij)^ (det Dw)"^ tr (rt(D^)-i(x, 0)r'(Z}^)-^(0, x)) e"^. (87) 

3.2.4- Quenched approximation 

Fermionic observables like Crr' require the evaluation of the fermionic determinant. This 
determinant can be expanded in fermionic loops. As an expansion parameter, one can use for 
instance the fermion mass. This gives rise to the so-called hopping parameter expansion. The 
leading order in this expansion simply consists in neglecting all fermionic loops, which means 
setting det Dw = 1- This infinite fermion mass limit defines the quenched approximation. In 
practical terms, working in the quenched approximation means neglecting the back-reaction of 
fermions on the gauge fields. 

In lattice QCD, in order to obtain an accurate numerical result it is crucial to evaluate the 
full fermionic determinant. However, in large- iV QCD, at fixed quark mass fermion loops become 
less and less relevant. In fact, the quenched large- limit coincides with the large- iV limit in 
the theory with dynamical fermions. Hence, if we are interested only in the large-A^ result, the 
calculation can be performed in the quenched theory (note that, obviously, this does not hold for 
the evaluation of finite- corrections). 

A phenomenon that has been observed in lattice QCD is the delayed onset of unquenching 
effects: the quenched calculation proves to work even in a regime in which one would expect 
significant contributions from fermion loops |116] . This feature can be seen as another indication 
that the physical strong interaction is close to its large- A^ limit. 

3.3. Monte Carlo calculations 

In a SU(A^) gauge theory, the evaluation of the path integral on a lattice of size S (in lattice 
units) involves the evaluation of d = 4(A^^ — 1)S integrals. Such a large number of integrals 
are impractical to be performed using grid methods. A stochastic approach, based on the ob- 
servation that the path integral measure is reminiscent of the Boltzmann measure in statistical 
mechanics — and, hence, at fixed bare parameters only a subset of possible values of the variables 
will dominate the integral — is preferable. Moreover, grid methods are affected by a systematic 
error that is 0(1/5"/"'), where s is a number that depends on the approximation used by the 
given grid method (for instance, s = 4 for the popular Simpson method). Hence, when d is large, 
the error becomes unavoidably of order one. Monte Carlo methods provide the wanted stochastic 
approach: in a Monte Carlo calculation, field configurations are generated according to the path 
integral measure, which means that configurations recur according to the weight of their contribu- 
tion to the path integral. Thanks to this property, the vacuum expectation value of an observable 
can be computed as a simple average over the configurations generated during the Monte Carlo 
simulation. Moreover, the error (which in this case is statistical and not systematic) can be kept 
under control, since it scales as 1/a/A'x, where A^^ is the number of generated configurations. 

The theory behind Monte Carlo calculations is based on Markovian processes. Consider a 
system that evolves through a sequence of states over discrete time. We indicate a generic state 
as Cm and the ensemble of all states (or configurations) as {Cm}- The evolution is determined 
by a probability Pnm to transition from C„ to Cm ay any given time. The dynamics is said to be 
Markovian if the configuration realized at time t only depends on the configuration realized at time 
t — 1, and not on the configurations realized at previous times. Under some technical assumptions 
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that we do not specify here, one can prove that an asymptotic probabihty distribution characterizes 
a Markovian dynamics. This is called the equilibrium distribution. Monte Carlo algorithms are 
recipes to construct Markovian dynamics that have the path integral measure as the equilibrium 
distribution. In general, several different Markovian dynamics can generate the same equilibrium 
distribution. For some Markovian dynamics that satisfies the detailed balance relation 

p p-S{C^) _ p p-S{C^) /ooN 

where S{C) is the action evaluated on the configuration C, the equilibrium distribution pm is 
given by 

g~'5'(Cm) 

in which it is easy to recognize the Boltzmann equilibrium distribution. 

The problem of generating an ensemble of configurations dominating the path integral and 
approximating the observables with controlled precision becomes then the problem of defining an 
appropriate Markovian dynamics to our theory. Once we have done this, we can start from any 
arbitrary state and let the system evolve. After discarding a sufficient number of configurations 
at the beginning of the chain, the remaining ones would be distributed with the right statistics 
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An algorithm can be characterized by the efficiency with which it explores the configuration 
space. A good measure of the efficiency is the correlation time. In general, for a given observable 
O, 

(C>(i)C(t + r)) oce-*/"o, (90) 

where the correlation time To depends on both O and the chosen Markovian dynamics. At fixed 
observable, provides a measure of the efficiency of the algorithm: the smaller To, the faster 
the configuration space is explored. Topological observables are notoriously more difficult to 
decorrelate than local observables. An algorithm that decorrelates fast topological observables 
can be considered efficient in general terms. 

For SU(2) pure gauge theory, an efficient algorithm is the one proposed by Creutz |477j and 
later refined by Kennedy and Pendleton |478] . This algorithm belongs to the wider class of the 
heath-bath algorithms, i.e. of those algorithms in which detailed balance is obtained by generating 
the new configuration according to its Boltzmann weight. To further decrease the correlation time 
and increase ergodicity, it is possible to supplement the heath-bath with an overrelaxation step, 
in which the link variables are changed in such a way that the action remains constant [479^ I48U] . 

For SU(A^), a Kennedy- Pendleton or overrelaxation update can be performed by a sequence of 
updates on different SU(2) subgroups. The idea of updating the SU(2) subgroups was proposed 
by Cabibbo and Marinari j481j . Although the Kennedy- Pendleton algorithm with a Cabibbo- 
Marinari cycle works well for SU(3), at very large N cycling over the SU(2) subgroups may 
become inefficient. To overcome this potential issue, it has been proposed to supplement the 
Cabibbo-Marinari update with an overrelaxation step over the whole SU(A^) group [4821 1483] . 

After a sequence of thermalized updates, the expectation value of an observable can be com- 
puted as the average of the values of the observable evaluated over the various configurations of 



^^A correct identification of the transient requires an analysis of the numerical stability of the observables against 
the length of this cut. 
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the Markov chain. For instance, if 0{i) is the value of the observable O at the i-th configuration 
in a Markov chain of length Nk, we define 




(91) 



(O), the vacuum expectation value of O, is found as 



{O) 




Nk- 



(92) 



Onk is a controlled estimate for (O), in the sense that 



A{0) 



) 
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(0)2 - O 
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(93) 



{O) 



K 



Therefore, in a numerical simulation, (O) can in principle be obtained with the desired accuracy, 
by tuning the length of the Markov chain. The statistical uncertainty AO can also be quantified. 
However, due to the correlation of the configurations, simple Gaussian statistics is not applicable 
in a straightforward way. Nevertheless, Gaussian statistics can be applied after the values of the 
observable have been averaged over bins of size equal to or larger than tq- Practically, this means 
that the number of independent estimates of O is not Nk-, as one would naively think, but rather 



Another issue that needs to be taken into account in order to provide a reliable estimator for 
{O) is the bias. Technically, one says that an estimate is biased if it does not agree with the 
analytical value when the latter is computable. While any bias will disappear for Nk oo, a bias 
could appear for naive estimates at finite Nk- There are standard methods to remove biases, the 
most popular ones being jack-knife and bootstrap. Although the details of those methods will not 
be discussed any further (we refer to ref. |484j for a pedagogical introduction), it is important to be 
aware that extracting a numerical value and the corresponding error for an observable measured 
in numerical simulations using Monte Carlo methods requires a careful analysis, the reason being 
that the data are correlated. This also applies to fits of Monte Carlo data [485J. Modern lattice 
calculations, including those discussed in this review, use robust procedures to estimate errors on 
observables and fitting parameters. The particular procedure used in each calculation is generally 
discussed in the corresponding original publication, to which we refer for details. 

3-4- Lattice simulations ofSl]{N) gauge theories 

Monte Carlo calculations in lattice gauge theory aim at computing numerically values of phys- 
ical observables. Typical quantities that are computed are connected correlation functions of two 
operators at zero momentum: 



where Nt is the number of sites in the Euclidean-time direction of the lattice, and the operator 
is averaged over all the other three spatial coordinates, to project onto zero-momentum states. 
We have assumed that the Euclidean-time direction is compact, and that periodic (antiperiodic) 



Nk/to- 



^(^) = F E ((^IW^2(t + r)) - {0\{t)){0,{t + r))) , 



Nt-l 



(94) 



t=o 
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boundary conditions along it are imposed for bosonic (respectively: fermionic) fields. For asymp- 
totically large Nt, 



C{t) oc e"'"", (95) 

where m is the mass of the lowest-lying state that connects Oi\0) to C2IO). If Oi are traces of closed 
loop operators carrying well-defined J^'-' quantum numbers (note that in this context J refers to 
the dihedral group, to which the group of spatial rotations is reduced by the lattice structure), 
C(r) will identify glueball states. If the Oi operator is a fermion bilinear with well-defined J^*", 
then C will be saturated by mesonic states. Finally, if the Oi are Polyakov loops wrapping around 
a compact spatial direction (with periodic boundary conditions), then the propagating states will 
be torelons, which will have a mass that is proportional to the string tension. Note that, in 
practice, extracting a signal with sufficient accuracy for a meaningful determination of the mass 
can be challenging. This problem can be overcome, by building various operators with the same 
quantum numbers that can be used in a variational approach (see e.g. ref. |486j ). The variational 
calculation also allows to extract masses of excitations in the given channels. 

Another important remark is that the quality of the numerical data does not decrease with 
N: simple large- counting arguments show that the ratio noise over signal is constant in for 
pure gluonic correlators, while improves as l/\fN for fermionic correlators. 

Thermodynamic properties can be studied in a finite temperature setup, in which the lattice 
has Nl X Nt sites {Ng is the number of sites in each spatial direction, while Nt is the number of 
sites in the temporal direction), with Ng ^ Nt, and periodic (antiperiodic) boundary conditions 
for bosonic (fermionic) fields along the compact direction, whose size is related to the temperature 
T by T = l/[aNt). The deconfinement temperature is identified by looking at the peak of the 
susceptibility of the order parameter, which for a deconfining transition in the Yang- Mills theory is 



the Polyakov loop^® The scaling of the position (in /3) of this susceptibility as A^s — ?■ 00 allows to 
identify the critical value of /3 (which is a function of Nt and of the order of the phase transition). 
The framework that enables us to perform those studies is the theory of finite size scaling. 

Numerical studies of a lattice gauge theory with Monte Carlo techniques involve the following 
steps: 

1. on a lattice of fixed size S and at a fixed value of the lattice coupling (3 (and /c, in the 
presence of dynamical fermions), evaluate numerically the vacuum expectation values of 
operators corresponding to physical observables; 

2. at the same lattice couplings, perform numerical simulations on lattices of larger volume, in 
order to extrapolate to the thermodynamic limit; 

3. repeat these two steps at various couplings in order to determine the value of the observables 
in the continuum limit. 

In the programme, there are two extrapolation processes. Both of them involve fitting the Monte 
Carlo data to analytical behaviors whose leading order in the correcting parameters is known. For 
instance, on a lattice of linear size L, in the chiral limit the mass of the pseudoscalar meson 
receives finite size corrections that, at the leading order, are proportional to e~^^^ |487j . Fitting 
the lattice data on various sizes with this Ansatz in the region in which the asymptotic behavior 



^^Strictly speaking, the Polyakov loop is a well-defined order parameter only on an infinite lattice, because there 
can not be any phase transition in a system with finitely many degrees of freedom. For this reason, one normally 
considers the modulus of the average Polyakov loop in a configuration as a (pseudo-)order parameter, and identifies 
the location of the critical point by monitoring its susceptibility, when the system parameters are varied. 
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is reached (which has to be determined as a part of the simulation) allows one to extract m^r in 
the infinite volume limit. 

Asymptotic freedom dictates that the continuum limit of the theory is reached for /5 — )■ oo. In 



fact, at the lowest order, perturbation theory^ predicts that 



1 122lLr 

a = —e'TiN^^ (96) 

(up to subleading corrections), where Ai^^t denotes the dynamically generated mass scale in the 
lattice regularization scheme. This formula implies that the lattice spacing goes (exponentially) 
to zero when /3 — ?■ oo. The variation of an observable with a is predicted by the Callan-Symanzik 
equation. The existence of a well-defined continuum limit implies that for two observables of the 
same mass dimension Oi and O2 

lini 2^ = gi, (97) 

where Oi = a~'^^Oi and di is the mass dimension of Oi. At the leading non-trivial order in a, near 
the continuum limiiP^ 

J = ^ + 0(aV), (98) 
O2 ^2 

where denotes an observable with mass dimension 2 |490[ 1491] . This expression implies that, 
asymptotically, lattice corrections are quadratic in the lattice spacing. Values of /3 for which 



observables fulfill eq. (98) are said to be in the scaling region. 

We stress again that, in order for the extrapolations described above to be meaningful, the 
system must be in the correct regime. For instance, in a small volume, deconfinement might arise 
and as a result the numerical data extracted in this phase are not simply related with their infinite 
volume limit. Similarly, the limit corresponding to strong (bare) lattice gauge coupling (/3 — )■ 
0) is heavily dominated by discretization artifacts, whose properties (like, e.g., the existence of 
confinement |2] and the finiteness of the mass gap |492j . which were proven analytically in the early 



days of lattice QCD) are not necessarily relevant for the continuum theory Nevertheless, strong 
coupling expansions are a useful theoretical tool in lattice gauge theory, and are characterized by 
a finite convergence radius |498l 1496] . For the Wilson discretization of SU(iV) Yang-Mills theories 
in four spacetime dimensions, it is known that, for N > 5, the range of (weak) couplings, which is 
analytically connected to the continuum limit, is separated from the strong-coupling regime by a 
strong, first-order bulk transition — which is signalled by a discontinuity in the expectation value 
of the plaquette, and which is not related to any symmetry breaking pattern — , while for N < 5 
a crossover connects the strong- and weak-coupling regimes |499j . In order to probe the region 
of couplings analytically connected to the continuum limit, the simulations have to be performed 



^'^For SU(iV) gauge theories, in the lattice scheme the Symanzik /3-function, which determines the variation of a 
as a function of the coupUng go, is known to three loops |488l 1489] . 

^*We are assuming that boundary conditions do not introduce corrections that are proportional to a, which is 
not always the case, but it is true in most simulations — including, in particular, those that use (anti-)periodic 
boundary conditions. 

^^Indeed, in the limit of strong bare lattice gauge coupling, even compact U(l) lattice gauge theory is confin- 
ing [493) — while it is in a Coulomb phase at weak coupling [4941 14951 HMl I497j . 
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at values of bare lattice 't Hooft coupling Aq smaller than the critical value corresponding to the 
bulk transition leading to the strong-coupling regime — which, for — )■ oo, has been numerically 
estimated to be at 0.3596(2) [500] Q 

SU(iV) Yang-Mills theory dynamically generates a scale. The whole physical spectrum can be 
expressed in terms of this scale. Hence, in order to meaningfully compare theories at different N, 
a scale needs to be fixed. This could be for instance the string tension a (i.e. the asymptotic 
slope of the confining potential at large distances) or the deconfinement temperature Tc. Once 
the choice has been made, large- iV arguments predict the scaling with of all other quantities 
relevant for the infrared dynamics. In particular, if a well-defined large-A^ theory exists, all 
spectral quantities should have a finite large- limit. From perturbation theory, we expect that 
the leading finite- corrections are of order 1/A^^ for the gauge theory, and of order 1/A^ for the 
theory with dynamical fermions. Taking further the perturbative argument, one would expect 
that pure gauge observables can be expressed in a power series in l/N"^, while in the presence 
of dynamical fermions the power series is in 1/A^. In this language, the proximity of SU(3) to 
SU(oo) means that the series converges for A^ = 3. Moreover, the large- A^ approach is useful to 
describe QCD if for a comprehensive set of observables a reasonable approximation (to the order 
of few percents) can be obtained by retaining only few leading corrections, with the quality of the 
approximation systematically improving when higher-order corrections are added. Note that, in 
the continuum limit, the quantity of reference, that one uses for comparing the results in theories 
with a different number of colors A^, does not play any role: different quantities may be affected 
by finite- A^ corrections with different coefficients, but each of them tends to a well-defined value 
in the N —> oo limit. 

We conclude this section with a brief discussion about the question, whether the continuum 
and the large- A^ limits commute. As pointed out in ref. |503j . in general the interchange of these 
two limits may be non-trivial [504[ I505j — especially if there exists a set of degrees of freedom, 
whose number does not grow with A^, but which nevertheless have a strong effect on the dynamics 
at the cut-off scale at any finite A^. However, as discussed in refs. |503[ I506j . for the gauge theories 
that we are presently interested in, one can safely assume that the continuum and the large-A^ 
limits commute. In other words, if we want to study the theory at infinite A^, we can either take 
first the continuum limit at fixed A^ and then the large- A^ limit, or take first the large- A^ limit 
at fixed cut-off a and then the continuum limit. In the latter approach, the lattice spacing is 
kept fixed across the various A^, by simulating the various theories for /3 such that the value of 
a physical quantity (e.g. the string tension or the critical temperature) has a predefined value in 
units of the lattice spacing a. While performing the large- A^ limit at fixed lattice spacing should 
be seen as an intermediate step towards getting continuum large- A^ physics, this approach can 
prove convenient in calculations that are particularly demanding from the computational point 
of view. In addition to various examples in the continuum limit, in sec. [5} we will discuss some 
results for which the large- A^ limit has been taken at a fixed lattice spacing in the scaling region. 

4. Factorization, loop equations, and large-N equivalences 

Besides the phenomenological implications and the connections with string theory discussed in 
sec. [2| large- A^ field theories and statistical models exhibit many further interesting mathematical 
properties: the expectation values of products of physical operators factorize, up to 0(1/A^) 



^°It is interesting to note that this value is close to some estimates, worked out with (truncated) analytical 
expansions based on the large- limit, which were already obtained in the early 1980 's [5011 1502j . 
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corrections (see subsection 4.1 below), which suggests an analogy with the classical limit of a 



quantum theory (subsec. H^), and indicates the suppression of fluctuations for — )• oo. This 



led to conjecture that the large- dynamics may be determined by a master field (subsec. 4.3). 
The factorization properties also imply that one can formulate a closed set of equations for the 
expectation values of gauge invariant operators, which are presented in subsec. 4.4 For the 



lattice formulation of the Yang-Mills theory, Eguchi and Kawai discovered that these equations 
reveal the volume independence of the theory in the large- iV limit, so that, in four spacetime 



dimensions, the theory can be reduced to a single-site model of only four matrices (subsec. 4.5), 
provided that center symmetry remains unbroken. Since the latter condition is not satisfied at 
weak couplings (which are relevant to take the continuum limit), different variants of the original 



model have been proposed in the literature: these include the quenched (subsec. 4.6) and the 



twisted (subsec. 4.7) versions of the Eguchi-Kawai model, its generalization with dynamical adjoint 



fermions (subsec. 4.8) or with double-trace deformations of the Yang-Mills action (subsec. 4.9). In 
parallel to these theoretical developments (and to the related numerical studies), a complementary 



approach has been pursued in a series of works (discussed in subsec. 4.10) exploiting the partial 
volume reduction of the original Eguchi-Kawai model, down to the minimal lattice volume, in 
which center symmetry is unbroken. Finally, in subsec. 4.11 we discuss how large-iV volume 
independence and certain correspondences between different large- theories can be interpreted 
in terms of orbifold equivalences. 

4.I. Factorization 

Many of the mathematical simplifications characteristic of large- field theories and statis- 
tical models are immediately made manifest by the combinatorics of large- A^ counting rules. In 
particular, in the computation of expectation values of products of (appropriate) gauge-invariant 
operators Oi in large- A^ gauge theories, the latter imply that the leading (in A^) contributions 



come from disconnected terms, because they are associated with the maximum number of color 
traces, and, hence, with the largest number of independent color indices to be summed over: 



{Oi02) = {Oi){02)+0{l/N). 



(99) 



Eq. (99) shows that the large- A^ limit bears analogies with two different limits in quantum field 



theory. On the one hand, it can be interpreted as a thermodynamic limit |118j : for a system 



characterized by a finite correlation length, eq. (99) is analogous to the cluster decomposition 



property of statistical, volume averages in a large volume V, which holds up to 0(1/1^) corrections. 
An interesting facet of this analogy is that both A^ and V are related to the number of degrees of 
freedom of the system (so that the integration measure for the generating functional, or for the 
partition function of the system, depends on A^, or on V). 

4-2. A classical-mechanics description for the large-N limit 



A different interpretation of eq. (99) is motivated by the analogy with the classical limit of a 



quantum field theory [27j. As it is well-known, a classical system arises explicitly in the h ^ 
limit of a quantum theory, when one studies the behavior of a basis of coherent states. The latter 
form an overcomplete basis, which allows one to write all operators solely in terms of their diagonal 



^^Such operators include, in particular, local gauge- invariant purely gluonic operators, fermionic bilinear opera- 
tors, and Wilson loop operators. Examples of operators which, on the contrary, do not satisfy eq. (99) have also 



been pointed out in the literature [1181 1507j . 
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elements, and have vanishing overlaps in the /i — )■ hmit. Together, these two properties lead 
to the factorization of expectation values of products of operators, and the deterministic nature 
of the limit becomes clear from the fact that, for h ^ 0, coherent states are characterized by 
simultaneously vanishing uncertainties in conjugate variables. In this limit, the classical phase 
space can be defined as the manifold of the coordinates labeling different coherent states, the 
quantum Hamiltonian (or, more precisely: the set of its diagonal matrix elements — also called its 
symbol — in the coherent state basis) can be mapped to its classical counterpart (which is minimized 
by the solutions of the classical equation of motion), and the classical Lagrange function is then 
obtained by Legendre transform. 

This construction can be repeated for the large- limit (more precisely: for the — 



limit) of a family of statistical systems or of quantum field theories, by generalizing the familiar 
Heisenberg group to an appropriate coherence group |508j . generated by suitable "coordinates" 
and "momenta". Choosing a coherence group that, for each value of N, acts irreducibly on the 
Hilbert space JifN of the corresponding theory, and a reference state |0)^, one can construct 
generahzed coherent states |a)^ by acting with elements of the coherence group on |0)^. The 
irreducibility condition allows one to express any bounded operator as a linear combination of 
elements of the coherence group |509j . For gauge theories, the coherence group is generated by 
Wilson loops, possibly decorated by a (chromo-)electric field insertion. 

Assuming that the correspondence between linear operators in and their symbols on 
coherent states is injective, one can consider those operators Ojsf whose elements in the coherent 
state basis (appropriately normalized) have a smooth large- iV limit: 

3 U.n (100) 

and introduce an equivalence relation among coherent states, \a) ~ \(3), defined by the requirement 



that, for all operators satisfying eq. (100), one has: 



lim ((a|O|a)^-(/3|O|/3)^) = 0. (101) 



If representatives of distinct equivalence classes defined by this relation have exponentially sup- 



pressed overlaps for N ^ oo, and if the 7i/N operator satisfies eq. (100), then one can prove that, 
for N — )• oo, the original theory reduces to a classical mechanics theory defined on the coadjoint 
orbit of the coherence group, and H/N tends to the corresponding classical Hamiltonian [27j . 

We conclude this subsection with an important observation: while the discussion presented 
above shows that it is possible to construct a suitable mapping of a large- quantum field theory 
to a classical-mechanics system, it is important to remark that, in general, this does not imply 
that fundamental fields of the large- quantum theory are described by the classical equations of 
motion derived from their Lagrangian density. As discussed, e.g., in ref. [22], the reason is that 
the generating functional of the system features an exponent, that receives 0{N^) contributions 
not only from the action, but also from the measure, due to the integration over 0{N^) gauge 
field degrees of freedom. For a discussion on somewhat related concepts, see also ref. |510] . 



^^Following ref. [27], we refer to a family of statistical systems or of quantum field theories, rather than to just 
a statistical system or a theory with one parameter, in order to remark that also the very structure of the theory 
(including, in particular, its number of degrees of freedom) can depend on the parameter. 
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4-3. Spacetime independence and the master field 

The construction outlined above can be carried out to find an explicit large-iV solution for 
certain classes of models (e.g., vector models |511j and single- matrix models [HI I512j ). but not 
for the most interesting case of gauge theorie^^^] — with the exception of one-plaquette lattice 
models |514[ 1515] . While, at first blush, this latter case may seem completely unphysical, it turns 
out that the factorization properties expressed by eq. ([99]) also have potential implications for the 
(lack of) spacetime dependence of the large- theory. A first, intuitive discussion of this feature 



was already expounded in ref. [18]: for Oi = O2, eq. (99) shows that quantum fluctuations are 
suppressed in the large- limit, hence it is reasonable to expect that, for N ^ 00, the path integral 
is dominated by a unique gauge configuration (or, more precisely, gauge orbit). The latter was 
interpreted as a classical field by Witten [18], and later dubbed "master field" by Coleman [T9] . 

Explicit results based on the master field approach were derived for some simple models in 
refs. |516[ 1517] , while an algebraic equation for the master field was proposed in ref. |518] , in the 
form of a quenched Langevin equation. Poincare invariance of vev's then implies that the master 
field should also be Poincare invariant — possibly up to gauge transformations, meaning that there 
exists at least a gauge in which it is Poincare invariant. While, strictly speaking, the interpretation 
of the master field as a classical, c- valued field is not correct |118j . a more rigorous treatment can be 
formulated in terms of non-commutative probability theory (see refs. |519[ 15201 152H 15221 1523[ 1524] 
and references therein), as illustrated in refs. |525[ I526[ 1527] . 

4.4- Loop equations 

As discussed above, spacetime independence of the master field is related to the large-A^ 



factorization of expectation values of products of gauge-invariant operators, expressed by eq. (99). 
Another important consequence of the same equation is that it implies that, for A^ — )■ 00, one 
can derive a closed set of Schwinger-Dyson equations for expectation values of physical operators, 
such as traces of Wilson loops along a contour C in a gauge theory [9] , 



^TrWiO = iTr 



Pexp (^ig ^ A''^{x)dx^ 



(102) 



(where V denotes path-ordering), which allow one to reformulate the theory in a gauge-invariant 
way |528] . Such geometric reformulation of gauge theories in terms of "loop calculus" [22] maps the 
familiar objects appearing in gauge theories, like (non-)Abelian phase factors, covariant derivatives 
and field strength to loop functionals, path derivatives and area derivative, respectively — see 
ref. [22] for a pedagogical introduction to the subject. Loop calculus deals with the Hilbert 
space of square-integrable functions which describe closed loops C (up to reparametrizations), 
and functionals thereof, J-'{C). For the latter, the operations of area and path derivative can 



be defined^^ by considering the variations obtained by deforming C through the addition of an 
infinitesimal loop or, respectively, a backtracking infinitesimal path at one of its points. 

Taking the area derivative (sometimes also called "keyboard derivative") of a Wilson loop at 
a point corresponds to inserting the field strength at that point [529]: 

^ ^-JtW{C) = ^Tr [F^p{x)W{C)] . (103) 



6aai3ix)N ' ' A^' 



^•^In this context, the difhculty of gauge theories is related to the representation of their fundamental fields — see, 
e.g., ref. |513] for a discussion. 

^''We restrict our attention to functionals of the Stokes' type, for which the area derivative is well-defined, i.e., 
independent of the shape of the infinitesimal loop. 
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Acting with the path derivative on a Wilson loop decorated by the insertion of a local operator O 
at a "marked" point x has the effect of replacing O with its covariant derivative (evaluated at x): 



^l^Tr[0{x)W{C)] = ^Tt[{D,0{x))W{C)] 



(104) 



(on the other hand, for functionals of the Stokes' type, the path derivative is identically vanishing 
at "regular", i.e. non-marked, points along a loop). 

The loop calculus outlined above can be used to derive the equations describing the invariance 
of the following functional integral: 



VA—TT{T''W{C))exp 



dt d-^x C 



(105) 



(in which the base point of the loop C is taken to be x, i.e. T° is inserted at the point x) under 
(functional) variation of the components of the gauge field, A''^{x): the "source" terms obtained 
when the variation is applied to T^W{C) must cancel against the "equation-of-motion" terms 
obtained applying it to the exponential. Upon contraction of the color indices according to the 
Fierz identities of the color algebra, this results in the equation [9]: 



6 



"Sao^^ix) \N 



^i:iW{C)^ = \jdyJ^''\ 



x-y) 



Trl^(a,,)W(C,,,; 



1 



N 



TiW{C) 



where C^^y and Cy^^ denote two complementary portion; 
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of the loop C. 



(106) 



Eq. (106) holds for every value of N, however it is not a dosed equation in loop space, because. 



generally, the product of traces appearing on the right-hand side is not a linear combination of 
single-trace loop operators. In the large-A^ limit, however, factorization allows one to write such 
term as the product of the expectation values of two single-trace operators (and to discard the 



last term in the first line of eq. (106) as subleading), leading to [S]: 



6 



^W(C)^ = \jdyj^'>\x - y) (^lw(a,,)^ (|^W(C,,.; 



which is a closed equation in loop space 



Eq. (107) is a geometric equatior 
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+ 0(l/iV^), 
(107) 



with deep implications]^ In particular, the construction 
of an iterative solution, starting from the expansion of the path-ordered Wilson loop W{C) in a 
series of cyclically ordered Green's functions, reproduces the set of planar diagrams which give 
the non-vanishing contributions to W{C) in the large- limit. 



^^Note that, if C is a self- inter secting loop, the J-distribution does not necessarily force Cx,y or Cy^x to be trivial. 
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Strictly speaking, eq. (107) should be formulated taking an appropriate, gauge-invariant loop renormaliza- 



tion procedure into account; this does not pose any particular problem, as it is possible to prove that the loop 
renormalization can be expressed through a purely multiplicative factor: this holds both for smooth |335j and for 

self- intersecting loops [530) . 

^^As discussed in ref. [55], the geometric aspects of eq. (107) can be exhibited by expressing it in terms of the 
Levy operator |531j . 

^^Note, however, that solving eq. (107) for all loops would not provide the complete data about the theory: 



in particular, physical quantities such as the mass spectrum and scattering amplitudes are encoded in connected 
correlation functions of gauge-invariant operators, which vanish if exact factorization holds, and hence cannot be 



captured by eq. (107) 
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At the non-perturbative level, an explicit solution of eq. (107) is not known; in ref. |532j . it was 
shown that a confining area-law Ansatz for asymptotically large loops is consistent with eq. (107), 
but the corresponding string tension could not be computed. Interestingly, the Nambu-Goto 
bosonic string action |533t 1534] (which describes confining Wilson loops in terms of fiuctuating 



surfaces, with an action proportional to their area — see the discussion in subsec. 5.1), which has 



been a candidate effective string model for the infrared QCD dynamics since the 1970's |535[ 15361 



1537] . is not consistent with eq. (107). A formal string solution to eq. (107) was discussed in 
refs. [5381 [539]: it takes the form of a fermionic string model, describing the dynamics of point-like 
Majorana spinor fields ("elves") on the surface bounded by the loop; the effective string tension 
was found to be related to the bare elf mass by a scaling law with critical exponent 12/11, but a 
complete solution of the theory was not found. 

4-5. The Eguchi-Kawai model 



One can also formulate equations, analogous to eq. (107), on the lattice |540] 15411 1542] : this 
led Eguchi and Kawai (EK) to discover the surprising property of volume reduction |543] , which, 
if the necessary conditions are satisfied, provides a concrete realization of spacetime independence 
at large iValready alluded to by the idea of a translationally invariant master field. In short, the 
statement of EK reduction is that, if: 

1. factorization of the vev's of physical operators, eq. ( [99] ), holds, and 

2. the global Zat center symmetry of the system (which, for the Euclidean time direction, has 



already been discussed in subsection 2.5 in the context of finite-temperature gauge theories) 
along each of the four Euclidean directions is not spontaneously broken, 

then the Schwinger-Dyson equations satisfied by vev's of (topologically trivial) Wilson loops in 
the large- SU(A^) Yang-Mills theory on the lattice are independent of the physical hypervolume 
of the system. As a consequence — assuming that these equations have a unique solution — also 
the physical observables of the theory are independent of the system hypervolume. Note that the 
second condition mentioned above arises from the fact that, due to gauge invariance, vev's of open 
Wilson lines of finite length are vanishing in the large- volume theory, while, in the reduced- volume 
theory, if the linear size of the system equals the line length, they vanish only if the global 
center symmetry (which tends to U(l)^ in the A^ — > oo limit) is unbroken. 

As a consequence, if the conditions above are realized, one could study the large-A^ lattice 
theory by reducing it down to a single-site lattice model, i.e. to a model of only four matrices 
(describing the parallel transporters along the four Euclidean spacetime directions — see section |3]). 
The importance of this observation is not purely academic: if the two conditions stated above 
are satisfied, then volume independence of the large- A^ theory could open up the possibility of 
studying its non-perturbative dynamics analytically, by reducing the equations of the quantum 
field theory in an infinite spacetime to the equation of the theory reduced to a single point, i.e. 
to the Schrodinger equation of ordinary quantum mechanics |544j . 

4-6. The quenched momentum prescription and the quenched Eguchi-Kawai model 

Even though a theoretically appealing idea, it was almost immediately realized that large- A^ 
volume reduction for pure Yang-Mills theory along the lines of the original EK proposal could not 



work The reason is the failure of the second necessary requirement listed above: by reducing 



^^This was already acknowledged in a note added to the original paper |543j . mentioning the result of the work 
by Bhanot, Heller and Neuberger |545] — who found evidence for the spontaneous breakdown of center symmetry 
in the regime of couplings connected to the continuum limit, i.e. in the perturbative regime — and similar findings 
by Wilson and by Peskin. 
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the lattice sizes, center symmetry gets spontaneously broken (at least in the range of couplings 
relevant for approaching the continuum limit). This can already be seen at the perturbative level: 
denoting the the eigenvalues of the [/^ matrices in the EK model as exp(26'^) (with 1 < a < N), 
the one-loop effective potential experienced by the 6*^ phases turns out to be attractive for all 
spacetime dimensions D > 2, leading to spontaneous breakdown of center symmetry |545t 1546] . 
The Monte Carlo simulations reported in ref. |547] confirmed this also non-perturbatively. 

A first, possible solution to preserve center symmetry was proposed in ref. |545] : the quenched 



EK model, ^ The quenched EK model is based on the idea of computing expectation values in 
the reduced model at "frozen" values of the ^JJ's, and then performing an average over such values 
according to a suitable probability distribution for the 6'^'s, which is chosen in such a way, that 
center symmetry is explicitly enforced. 

The way volume independence could arise in the EK model, the way it fails due to the spon- 
taneous breakdown of center symmetry, and the way this problem is avoided (at least at the 
perturbative level) in the quenched EK model can be understood via a suitable mapping of the 
degrees of freedom associated with different four-momenta in the original planar theory in a large 
volume to the matrix entries of the one-site model |553t 15541 15551 1556[ 1557] . This is based on the 
following, elementary observation: for a generic theory (either one with global or with gauge invari- 
ance under a U(A^) symmetry group) with fields in a two-index representation of the group, like, 
e.g., the adjoint representation, the fundamental lines in a generic Feynman diagram expressed 
in double-line notation are closed and non-intersecting. Hence, momentum conservation at each 
vertex of the diagram is automatically satisfied, if one associates a generic four-momentum to 



each fundamental line (with index a), and takes = ~ Pjl to be the momentum associated 



with the propagator obtained from the (oppositely oriented) lines of indices a and b. In view of 
this observation, one can prove that, by replacing each occurrence of idf^ in the original action of 
the continuum theory with the adjoint action of the matrix: 

P^ = diag(p;,pJ,...,p;r), (108) 

(or, equivalently, by replacing the finite translation operator exp(a9^) with exp(ia[P^, ■]) in the 
lattice theory), one obtains propagators which are equivalent to those of the original theory [5541 
1555] . This prescription, which goes under the name of "quenched momentum prescription" [554] , 
removes the spacetime dependence of the original matrix field from the theory, turning the action 
of the original theory into a function of a constant matrix and of the diagonal matrices P^, so 
that, for example, the vev of a generic observable O could be obtained by computing, first, its 
value at fixed momenta P^, and then integrating over the distribution of the P^'s. The theory 
obtained from the quenching momentum prescription has the same planar limit as the original 
one, at least at any finite order in perturbation theory |554] . 



Note that, by identifying each id^ with the corresponding P^ matrix defined in eq. (108), 
applying the quenched momentum prescription corresponds to representing the group of spacetime 
translations in terms of the degrees of freedom of the U(l)^ diagonal subgroup of the internal 
U(A^) symmetry group (this holds for either a global or a gauge symmetry). 

For gauge theory, restricting our attention to pure Yang- Mills theoryp] the quenched momen- 
tum prescription reads: 

A^{x) = exp{iP^x'')A^exp{-iP^x''), iD^, ^ P^^ + A^; (109) 



^"Related approaches were also discussed in refs. |548L 15491 1550] . while a Hamiltonian version of the quenched 
EK model was discussed in refs. |55H I552| . 

'^^ Generalizations of the quenched EK model to include quarks (in the Veneziano limit) were discussed in refs. |558l 
[55m [550] . 
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when applied to the lattice formulation of the theory, using = exp(mP^), this yields: 

2N 

l<a</3<4 

which (by a change of variables: Ua — Uo-D]^) is equivalent to the Eguchi-Kawai model 

2N 



(110) 



(111) 



l<o</3<4 



Note that the classical solutions minimizing S are arbitrary diagonal unitary (or special unitary, 
in the case of SU(A^) gauge group) matrices. 

The partition function and expectation values of the original lattice model are obtained by 



computing the corresponding quantities in the reduced model, described by eq. (Ill ), at fixed P^'s, 
and then averaging over the distribution of the -D^'s. In taking the latter step, it is convenient to 
express the Haar integration measure over each in terms of the eigenvalues exp(mp^), which 
yields a Vandermonde determinant A(Z}^): 
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k=l 



27r 



(112) 



(where, in the case of a special unitary gauge theory, the phases of the eigenvalues of each ma- 
trix are constrained to sum up to an integer multiple of 2tt). The dynamics of the EK model 
in the weak-coupling regime can then be investigated by studying the effect of quantum fiuctu- 
ations around a classical solution in a leading-order perturbative computation. This requires an 
appropriate gauge-fixing (e.g. to Feynman gauge) and the corresponding Fadeev-Popov determi- 
nant, so that the integration measure for the -D^'s (in D spacetime dimensions) turns out to be 
proportional to |545[ 1546] : 
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(113) 



As anticipated above, eq. (113) shows that a uniform distribution of the p^'s is a stationary point 



of the argument of the exponential, but one which corresponds to a minimum of the statistical 
weight for D > 2. As a consequence, quantum fiuctuations lead to the collapse of the p^'s, and, 
hence, to spontaneous breaking of center symmetry, which, in turn, invalidates the correspondence 
of the reduced model to the large-volume theory. 

A possible way to prevent spontaneous symmetry breaking (at least perturbatively), consists 



in modifying the integration measure for the UJs so that the dependence on the PJs cannot 



be eliminated. In ref. |554j it was proposed to do this by requiring the eigenvalues of the lattice 
covariant derivative U^D^ to be equal to those of D^, for each /x, i.e. by imposing the constraints: 



(114) 
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An alternative possibility consists in modifying the action |549) . 
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which are exphcitly gauge- invariant This leads to the integration measure: 



(115) 



which includes integration over the degrees of freedom of the unitary matrix that maps to 
U^Dfj^ by a similarity transformation. Shifting — )■ U^Djj^, the integration over the f/^'s can be 
immediately performed, leading to the reduced action: 
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J2 ReTr \v^D^V^VpDf,VlV^DiV^VpDlV^ 



l<Q</3<4 



(116) 



and, finally, one is left with an integration over the p^'s with a uniform distribution, yielding the 
quenched EK model introduced in ref. |545] . Thus, at the perturbative level, the effect of the 



gauge-invariant constraints in eq. (114) is to replace the non-uniform weight for the p'^'s of the 



original EK model, eq. (113), which was the cause of spontaneous breaking of the center symmetry. 



with a uniform distribution for the p^'s, which explicitly enforces center symmetry. 

Although early numerical tests found evidence confirming the validity of volume reduction in 
the quenched EK model also at the non-perturbative level |56H 15621 1563[ 15641 1565[ [566], these 
conclusions have been recently disproven by the high-precision study reported in refs. |567[ 1568] , 
which also discuss analytical arguments for this failure. The subtle dynamical mechanism respon- 
sible for the breakdown of volume reduction in the quenched EK model can be understood by a 



careful inspection of eq. (114), which plays a pivotal role in the model, and that, when expressed 



for the redefined f/^'s appearing in eq. (116), takes the form of a polar decomposition of the four 



matrices: 



U„, = V„D.,Vl 



(117) 



The crucial observation is that eq. (117) fixes the eigenvalues of the matrices only up to per- 
mutations |55H, 1569] . In general, configurations of the reduced model, in which the eigenvalues 
of the f/^'s differ only by permutations, yield different values of the reduced action defined in 



eq. (116). Hence, if dynamical fiuctuations of the V^'s can produce tunneling events between such 



configurations, the quenched system will choose the energetically most favored configuration(s), 
introducing non-trivial correlations between the [/^ components along different directions ("mo- 
mentum locking"), and resulting in non-uniform sampling of the gluon momenta. In fact, such 
phenomenon does occur: at fixed coupling, there exist configurations mediating eigenvalue per- 
mutations, whose energy cost is not enhanced in the N oo limit. The factorization condition, 
which is necessary for volume reduction to hold, is generally violated in the presence of correlations 
between the f/^'s in different directions. This can be seen by considering, for example, Wilson 
lines of the form: Ma/s = {1 / N)TT{UaUi3) , with a ^ P, whose value on a typical momentum-locked 
configuration is a complex number with modulus 0(1) and with a generic, momentum-dependent 

phase. Averaging over the momentum values, (M^/?) will generally vanish, whereas ^M^/jM^^ 

in which the momentum-dependent phase drops out, remains finite. 

4-7. The twisted Eguchi-Kawai model 

A different variant of the EK model was proposed by Gonzalez- Arroyo and Okawa: the twisted 
EK model |570tl571] . It is based on the observation that the perturbative behavior of the EK model 



Alternative approaches, based on non-gauge-invariant constraints, were discussed in refs. [5501 1556] . 
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can be altered, in a way that avoids the spontaneous breaking of center symmetry, by modifying 
the boundary conditions with appropriate twist factors Zap in the center of the group |572l 1573] . 

Similarly to what happens for the quenched EK model, the twisted EK model can be written 
in a volume-independent form. This is done by representing the group of spacetime translations 
of the original theory in terms of an Abelian subgroup of SU(iV). The key observation is that, 
for a theory with fields in a representation of the group which is blind to the action of center 
transformations (i.e. a representation with zero A^-ality), such as the adjoint representation, it is 
possible to represent translations in a D-dimensional spacetime in terms of D traceless SU(A^) 
matrices, satisfying the 't Hooft-Weyl algebra j572j : 

T^Ta = Zai3'roXi3 = exp{27iinap/N)TaTfi, a,/3 e {1,2,. . .D} (118) 

(where is an antisymmetric D x D matrix, whose entries are integers modulo N), i.e., com- 
muting only up to an element of the center of the group (see also, e.g., ref. [574] ) . Starting from 
the lattice formulation of SU(A^) Yang-Mills theory, the twisted EK model can be obtained by 
defining the following products of strings of matrices at each site x: 

v{x) = l[r^^, (119) 

by replacing the U^{x) link variables with: 

U^{x) ^ V{x)U^V\x) (120) 
and finally by doing a change of variables: — > t/^Fj^. This leads to the following action |570ll571j : 

^TEK = -y Wp'^<U^UpUiul) + YL.c]. (121) 

l<a</3<4 

For suitable choices of the Zap twist factors, it is possible to prove that, at the perturbative level. 



the model described by eq. (121) satisfies the same Schwinger-Dyson equations as the theory 
defined in an infinite volume in the large- iV limit. In particular, since the classical solution is 
given by [/„ = Fq, in the weak-coupling limit an open Wilson line from the origin to a generic 
point X fluctuates around V{x) (up to a Zat factor, which depends on the shape of the line). One 



can prove that any two SU(A^) matrices satisfying eq. (118) with a twist factor different from 1 are 
traceless. Using this fact, it is possible to show that a sufficient condition for the traces of open 
Wilson lines in the reduced model to vanish is that V{x) commutes with all of the F^'s. Since, 
for the theory in a lattice of finite volume, the only "open" Wilson lines whose trace can have a 
non- vanishing expectation value are those winding around the lattice (like, e.g., Polyakov loops), 
the latter condition reduces to: 

XaUap/N e Z. (122) 

In four spacetime dimensions, a simple solution for the case when is a perfect square {N = L^, 
L e No) is given by the symmetric twist: rz^^ = L for all a < (3. This also reveals that, in general, 
finite- corrections in the twisted EK model can be interpreted as finite-volume corrections, with 
the four-dimensional volume scaling like = N"^. 



''^Here, we assume D to be even. 
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If the conditions that ensure the equivalence between the twisted EK model and the theory 
defined on a lattice of finite (large) volume are satisfied, expectation values in the twisted EK 
are trivially related to those in the theory with trivial boundary conditions (i.e., no twists): for 
example, vev's of Wilson loops of area A in the (a, (3) plane get simply multiplied by a z;^^^ factor. 

As mentioned above, for the simplest symmetric twist, the classical solutions of the twisted 
EK model are of the form: 



(123) 



where the are a set of D traceless, N x N special unitary matrices satisfying eq. (118); such 
configurations correspond to the absolute minimum that the action in eq. (121) can take, and 



are hence called "twist eating" configurations |573] . As discussed in ref. |575j . in two spacetime 
dimensions the r„ matrices can be identified with the "shift' 
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and "clock'' 
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u = exp(27ri/A^), 



(125) 



matrices, which satisfy the little 't Hooft algebra: 

SnCn = expi^lni / N)C N S N ■ 



(126) 



In four dimensions, analogous matrices can be obtained from tensor products of these matrices 
(taking to be a perfect square, N = L"^, with L G No, and replacing A^ with L in the definitions 
of the shift and clock matrices)]^ 



(127) 



A systematic classification of the solutions was given by van Baal in ref. |576] (see also ref. |577] ). 
which showed the existence of non-trivial configurations which could survive in the large- A^ limit. 

The distribution of eigenvalues of the matrices in the (unphysical) strong coupling limit is 
uniform over (four copies of) the U(l) circle, hence center symmetry is not broken in this limit. 
In the opposite limit, the configuration corresponding to the classical solution is invariant under 
I^L C Zjy. As a consequence, the expectation values of the four Polyakov loops in the twisted EK 
model are vanishing both at strong coupling and in the classical limit. To understand how the 
continuum limit is approached, a leading-order weak-coupling expansion around the twist-eating 
solution was first discussed in ref. |571] : it revealed that, while the propagators coincide with 
those of the lattice theory, the vertices are generally modified by momentum-dependent phase 



^Strictly speaking, this construction yields the twist eaters in a representation in which the twist tensor is skew- 
diagonah n = iLfj2 ® l2- However, n can then be brought to the standard form by an SL(4, Z) transformation. 

48 



factors, and the overall phase of a given diagram is related to the leg ordering. In the large- 
limit, this leads to strongly oscillating factors, which suppress all non-planar diagrams |571l 1578] . 
In a continuum formulation |579j . the presence of momentum-dependent vertices, with relative 
suppression of diagrams which differ by non-cyclic permutations of the legs, can be interpreted in 
the context of field theories defined in non-commutative spaces |58Ut I581[ 1382] (see also ref. ^583j ). 
which are relevant for certain low-energy limits of M- and string theory |584[ I585t 15861 1587j . The 
twisted EK model provided a way to regularize theories defined in such spaces and to study them 
non-perturbatively |588i 1589] I590[ l59T] , ^ In fact, this formulation can also be carried out at finite 
values of N; in contrast to the quenched EK model (which only captures the dynamics associated 
with planar graphs) the twisted EK model is well-defined order by order in an expansion in powers 
of 

Non-perturbative Monte Carlo studies of the twisted EK model have been carried out since 



the 1980's; for Yang- Mills theories in four spacetime dimensions,^ early works include refs. |571l 
15971 15981 15991 16001 leOTl [6021 16031 MM- These articles found numerical evidence that the model, 
with the simple symmetric twist described above, correctly describes the physics of the large- 
N theory in a large volume, and reported results for quantities like T^/Al, the deconfinement 
temperature in units of the lattice A-parameter. Recently, however, some accurate numerical 
studies disproved these claims |605l 16061 1607] . showing that, although the center symmetry in 
the twisted EK model is preserved both in the weak- and strong-coupling limits, it does get 
spontaneously broken at intermediate couplings, and the width of the range of couplings in which 
this occurs grows with A^, making the investigation of the correct continuum physics at large A^ 
challenging — if possible at all. The spontaneous breakdown of center symmetry at intermediate 
couplings appears to be due to the system getting stuck in metastable phases, characterized by the 
fields fluctuating around center-symmetry breaking configurations (with action 0{N) above the 
twist-eating configurations, and separated from the latter by effective potential barriers of order 
0(A^^)), that become arbitrarily long-lived for A^ — )• oo. 

A possible solution to this problem was recently proposed in ref. |608] : the idea is to lift 
the action of the center-symmetry breaking configurations 0{N'^) above the twist-eating ones, 
by changing the twist to: = kL for all a < P, with k = 0{L) (in particular: k > L/4) 
and co-prime with A^. Non-perturbative investigations of this new formulation of the twisted EK 
model were initiated in ref. [608], and encouraging new results have been recently reported in 
ref. j609j . An example is shown in fig. |3| taken from ref. |609] . in which the value of the string 
tension (i.e., of the asymptotic slope of the confining potential at large distances) extracted from 
simulations in this new version of the twisted EK model at A^ = 841, and suitably extrapolated 
to the continuum, is compared with those reported in ref. |610j , from a conventional lattice study 
in a large volume: the result obtained from the twisted EK model simulation is consistent with 
the extrapolation of the latter. 

However, it is worth mentioning that the numerical study of ref. [606] . which also discussed 
the case of a twist tensor n = 0{N), concluded that the critical value of 't Hooft coupling, 
below which center symmetry is not spontaneously broken, tends to zero like 1/A^ when A^ is 
increased. Similar conclusions were obtained in another numerical study j607j . that used a twist 
tensor with /c = (L + l)/2. To completely clarify the viability of the twisted EK model with 
a non-minimal twist to investigate the behavior of large- A^ Yang-Mills theory, further numerical 
studies are needed. 



^^Alternative regularizations for field theories defined in non-commutative spaces are reviewed in refs. 159211593] . 
^^The twisted EK model has been studied numerically also in two [594. .595. ,591] and in six dimensions |596j . 



49 



CO 



< 



0.56 
0.55 - 
0.54 - 
0.53 - 
0.52 
0.51 
0.5 



SU(Ar) I H 

TEK iV=841 I — ^ 
0.515(3) +0.34(l)/iV2 



0.02 0.04 0.06 0.08 0.1 0.12 

Figure 3: A comparison of the results of simulations of the twisted EK model performed in ref. [609J using the 
twist formulation discussed in ref. |608j (green symbol), with those from lattice simulations performed in a large 
volume (in the conventional formulation of the theory, without exploiting the volume-independence properties), 
taken from ref. [610j (red symbols), and their extrapolation as a function of the number of colors (blue line). The 
quantity which is plotted is the ratio (extrapolated to the continuum limit) between the Aqcd parameter in the 
MS scheme, and the square root of the string tension, i.e., of the asymptotic force between static color sources at 
large separations. 

4-8. Volume reduction with adjoint fermions 

A different method to preserve center symmetry in the EK model, which has been proposed 
in the hterature, consists in modifying the theory by adding one or more flavors of dynamical 
massless Majorana fermions in the adjoint representation of the gauge group, and with periodic 
boundary conditions pllj,^^ When the theory is compactified on a small spatial torus of length 
R, a perturbative calculation of the effective potential Kfffor the phases of the eigenvalues exp{i9a) 
of the Polyakov line along that direction shows that the effect of fermions can (over-) compensate 
the symmetry-breaking terms coming from gluons |321l 16151 1616] : 
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where f{x) = [x mod 27r) — vr. In the naive EK model [nj = 0), the minimum of the potential 
is obtained when the term between the square brackets (which is positive semi-definite, because 



'^^This idea is related to circle compactifications and supersymmetry-preserving deformations in supersymmetric 
gauge theories [SHI [ST51[5Ti] . 
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< TT for any x) vanishes: this corresponds to the case when all eigenvalues collapse the same 
value, inducing spontaneous center-symmetry breaking. By contrast, the Uf = 1 case yields a flat 
one- loop effective potential (this case corresponds to A/" = 1 supersymmetric Yang-Mills theory), 
while for 1 < < 6 the inclusion of adjoint fermions in the theory has the effect of favoring vacua 
in which center symmetry is preserved]^ Analogous formulae hold in the case when two or more 
directions are compactified to small tori. 

While the inclusion of adjoint fermions in the model may appear artificial, it should be noted 
that, in the large- limit, the orientifold equivalence relates QCD with adjoint fermions to QCD 
with Dirac fermions in a two-index, symmetric or anti-symmetric representation |617[ 16181 \Q19\ 
1620] . Given that in real- world QCD the quark fields are in the fundamental representation of the 
SU(3) gauge group, which is equivalent (up to charge conjugation) to the two-index antisymmetric 
representation, the latter theory can be regarded as a natural generalization of QCD to the large- 
A^ limit. Incidentally, we mention that a quenched lattice study of the orientifold equivalence 
(in a large volume) was reported in ref. |621j : it was found that the masses of vector mesons 
corresponding to quarks in different representations tend to compatible values in the chiral and 
large- A^ limits, in agreement with the theoretical expectations. 

Following the proposal of ref. |611j . various works investigated the EK model with adjoint 
fermions and variants thereof, both analytically and numerically |622l 16231 16241 16251 16261 16271 
I628[ I629[ 16301 1631] . The results of the most recent non-perturbative studies via lattice simulations 
indicate that EK volume reduction with adjoint fermions works as expected, both with Uf = 1 and 
Uf = 2 Dirac flavors: as an example, ref. [626J found clear numerical evidence that center symmetry 
in the model with two flavors is preserved, for a finite (and rather large — see the discussion in 
ref. [623]) range of bare quark masses, over an interval of couplings which appears to extends all 



the way to the continuum limit, see fig. 4^ For the case with one adjoint quark flavor, a similar 
study was reported in ref. |625j . 

After establishing that volume reduction in the EK model with adjoint fermions does work, 
the natural next step in this agenda of lattice studies consists in investigating how cost-effective 
these simulations are. 

4-9. Volume reduction in Yang-Mills theories with double-trace deformations 

A related variant of the EK model was proposed in ref. |632] : it leads to a one-loop effective 
potential of the type appearing in eq. ( |128[ ) by deforming the ordinary lattice Yang-Mills action 
adding products of traces of Polyakov loopsj^ with suitably chosen (in particular: sufficiently 
large) positive coefficients, which suppress the weight of center-symmetry breaking configurations 
in the path integral. One advantage of this approach is that it allows one to reduce the extent of 
an arbitrary number of sizes of the system. For example, in order to preserve center symmetry 
when only one direction is compactified, the deformation term may be of the form: 

[N/2i 

J] J2 an|tr(L"(x))p (129) 

X n=l 

(where the summation is done over the points of a hyperplane orthogonal to the compactified 
direction), while more complicated deformations, which also include products of loops along dif- 
ferent directions, can be added, when two or more directions are compactified to small sizes. The 



■^^For n/ > 6, the theory is no longer asymptotically free. 

^"Related ideas have been considered in studies of the phase diagram of Yang-Mills theory at finite tempera- 
ture 
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Figure 4: Lattice simulations of large- iV gauge theories with Uf = 1 and Uf — 2 flavors of dynamical adjoint 
fermions in a small volume indicate that center symmetry is preserved in a wide region of quark masses, whose 
width remains finite in the continuum limit. The figure, taken from ref. |626| . is a sketch of the symmetry realizations 
in different regions of the space of the simulation parameters. In particular, the regions in which center symmetry 
is preserved in all the four directions are denoted by Zn*. The quantity on the vertical axis (b) is the inverse of the 
bare lattice 't Hooft coupling, so that (by virtue of asymptotic freedom) the continuum limit is obtained by taking 
b — )■ oo — while the region corresponding to small values of b is the strong-coupling region of the lattice theory. The 
quantity on the horizontal axis (k) is a simulation parameter related to the fermion mass: in particular, the latter is 
infinite for k — (corresponding to the pure Yang-Mills theory) , while it decreases when n is increased towards the 
dashed region, reaching zero at the critical value Kc- The figure, based on information obtained from simulations 
with Uf = 2 dynamical quarks in the adjoint representation of the gauge group in a small volume, shows that, in 
the presence of dynamical adjoint fermions, the spontaneous breakdown of center symmetry is avoided, for a rather 
wide range of quark masses, in a region that extends to the continuum limit. 



deformation terms also modify the expectation values of the observables, however it is expected 
that they only do so by 0{1/N) corrections, which become negligible in the large- limit. 

Note that the deformation terms appearing in eq. (129) are (1) non-local and (2) non-linear 
functions of the U^(x) fields on the lattice, and that (3) their number grows 0{N) (this latter 
requirement is necessary, in order to prevent partial breakdown of the 7jn symmetry down to a 
subgroup, if is a divisor of N). Even though these three features make the numerical simulation 
of the deformed model considerably more demanding than that of ordinary Yang-Mills theory, it 
should be noted that dedicated numerical algorithms for the simulation of this model are already 
available |636] . and Monte Carlo investigations are presently in progress |637j . 
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4-10. Partial volume reduction 

Finally, a different approach to EK volume reduction has been proposed by Kiskis, Narayanan 
and Neuberger in ref . |638] , and successfully carried out in a series of works, which investigated the 
large- iV limit of various observables, in two, three and four spacetime dimensions 



IMnl Em EM B, EH ElSl ElSl ElZ] (see also refs. [331 EIH] for reviews). The idea underlying 
this approach consists in reducing the lattice size down to the minimum value, for which center 
symmetry is preserved: this is possible because the critical size Ic at which the symmetry gets 
spontaneously broken (the "inverse deconfinement temperature") is a physical quantity with a 
finite, non-vanishing limit for N ^ oo. For this reason, in this approach one only achieves a 
partial volume reduction, and the number of lattice variables that have to be simulated still grows 
like (1/a)^ when the lattice spacing a is reduced towards zero to approach the continuum limit p] 
An advantage of this approach, however, is that it does not rely on any "trick" (whose validity 
may possibly fail at the non-perturbative level — perhaps through some non-trivial mechanism); 
moreover, its conceptual simplicity makes it straightforward to study numerically. 

In particular, following this approach, in ref. |638j it was shown that, in three spacetime 
dimensions, the distributions of eigenvalues for Wilson and Polyakov loops agree when extracted 
from simulations on a large volume lattice and on a partially reduced lattice. The analysis of the 
four-dimensional case was carried out in subsequent work, reported in ref. |482j , which studied the 
details of the phase diagram of the lattice theory, as a function of the lattice size in physical units: 
it was found that there exists a cascade of transitions, corresponding to the breaking of center 
symmetry in one, then two, three, and eventually in all four directions. This analysis was refined 
in ref. |640] . by addressing the loop renormalization |335] with a suitable smearing procedure, to 
get a well-defined continuum limit. In parallel with these works, the chiral symmetry realizations 
were studied in refs. |639l 1643] ; the results confirmed the vahdity of partial volume reduction at 
large N: for finite lattices of linear extent larger than 1^, chiral symmetry gets spontaneously 
broken, and a chiral condensate, independent of the lattice volume, appears. Finally, a series 
of works studied the confining potential and the string tension in the partial volume reduction 
approach |645[ 16461 1647j . In particular, this was done for the four-dimensional case in ref. |647j . 



from which the results shown in fig. 4.10 are taken: the plot shows the behavior of the confining 
potential \^ as a function of the distance r (the numerical values for both quantities are shown 
in the appropriate units of the lattice spacing a), as extracted from simulations with = 47 at 
inverse 't Hooft coupling b = 0.348, on a lattice of linear size L = 6a in all directions (which is 
small, but sufficiently large to preserve center symmetry). The potential, which was computed up 
to a distance of nine lattice spacings (i.e., one and a half the linear size of the lattice), shows the 
behavior that is expected in a large volume. 

4.11. Volume independence and large-N equivalences from orbif old projections 

Having discussed how the EK equivalence between theories defined in different spacetime vol- 
umes arises at large N (and how it can fail), we conclude this section presenting its interpretation 
in the broader context of a general class of equivalences among large- theories. Various types 
of correspondences between theories, which are different at every finite value of A^, but become 
coincident in the large- A^ limit, were already discovered during the early 1980's |649l 16501 1651] . 
but it is only during the last decade, that systematic methods to construct such equivalences (and 
to precisely identify the conditions for their validity) have been worked out. 



^^Note that, volume independence for lattices of linear sizes larger than 1^ is equivalent to temperature indepen- 
dence of physical operators in the confining phase |297| . 
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Figure 5: Computation of the confining potential V, as a function of the distance between static color sources 
r, using partial volume reduction. The figure shows the results obtained in ref. |647| for V (red circles), up to 
distances of nine lattice spacings, i.e. one and a half the linear size (L) of the lattice used in those simulations. A 
three-parameter fit (dashed red curve) to a Cornell-type potential is also displayed. 



One of the first papers pointing out tfie interpretation of EK volume reduction in large- iV lattice 
gauge theories in terms of these equivalences was ref. [503]. The discussion that we present here 
mostly follows the developments of this subject through a series of articles by Kovtun, Poppitz, 
tinsal and Yaffe [Hill EIHl ES21 EMI ESS ESni ESZ! • For historical reasons related to their 



origin in the context of string theory |658t 1659^ I66UI I661j , these equivalences are called "orbifold'' 
equivalences. 

In this context, an orbifold equivalence is defined as a mapping between a "parent" theory, and 
a "daughter" theory, which is obtained from the former by removing ("projecting out") the degrees 



of freedom which are not invariant under the action of a discrete global symmetry group, More 
precisely, in the large- iV limit these equivalences relate vacuum expectation values and correlation 
functions of a subset of gauge-invariant observables in the parent theory and in the daughter theory 
(possibly up to a rescaling of couplings and volume factors), provided that the global symmetry 
of the parent theory used to construct the orbifold projection is not spontaneously broken, and 
that possible global symmetries of the daughter theory, which interchange its gauge group factors, 



^ In string theory, equivalences of this type arise when considering the low-energy dynamics of D-branes defined 
in spacetimes which are orbifolds, i.e. manifolds admitting points whose neighborhood is locally diffeomorphic to 
a quotient of a real vector space by a finite group. 

^•^The degrees of freedom which are not invariant under the discrete symmetry group used for the projection are 
said to form the "twisted sector" of the theory. Note, however, that here the adjective "twisted" has no direct 
connection with our discussion of the twisted EK model. 
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are not spontaneously broken either. 

At the perturbative level, this equivalence is based on the observation that planar graphs of 
the parent and daughter theories coincide at all orders, and that generic correlation functions of 
gauge-invariant operators (from which one can extract physical observables, including the mass 
spectrum) obey the same set of closed equations. To extend the validity of the correspondence to 
the non-perturbative domain, the authors of ref. |652] focused on U(A^) gauge theories (possibly 
coupled to Hs species of massive scalar fields and n/ species of fermion fields, both of which 
in the adjoint representation of the gauge group) regularized on a lattice, and on the Migdal- 
Makeenko-like equations for gauge-invariant operators (which can be written as traces of closed 
Wilson loop^^ and their correlators. Consider first the pure Yang-Mills case (n^ = nf = 0): 
assuming that N is not prime (so that N = kn, with both k and n integers larger than 1) one 
can introduce an orbifold projection under a H = subgroup of the global U(A^) group of 
coordinate-independent gauge transformations, by requiring that a generic $ field variable of the 
parent theory be invariant under: 

$ ^ 7$7"\ (130) 

where: 

7 = r2 1„, = diag(l, w, w^, . . . , o;^""^), u = exp{2iTi/k). (131) 

The construction can be easily generalized to an orbifold under a subgroup (if N = k'^n), by 



introducing d matrices ja, and requiring invariance of $ under the transformations in eq. (132), 
for each 7^ separately. Matter fields can also be included in a straightforward way: for non-zero 
Hg and/or Uf, the global symmetry of the parent theory is enlarged to U(A^) (8>U(ns) ® U(nj), and 
each matter field can have a generic charge under the 7^ projection, so that the orbifolding 
condition reads: 

$ = exp {2mra{^)/k) 7a$7a (132) 

Under this projection, the gauge degrees of freedom of the parent theory are mapped to block- 
diagonal matrices, with k"^ blocks of U{n) matrices, and the gauge invariance group of the daughter 
theory is U {N)^ . The action for the daughter theory can then be written in the form of a standard 
Wilson lattice action (possibly coupled to the matter field terms), with a sum over the "sites" of 
a lattice, and with a trivial rescaling of couplings and volume factors. 

This mapping can be extended to all gauge-invariants observables of the parent theory, which 
are "neutral" under the orbifold projection, and, finally, one can derive the loop equations for the 
daughter theory. If the Z^ symmetry in the daughter theory is not spontaneously broken, these 
equations are equivalent to those in the parent theory, provided that the latter does not break 
the discrete symmetry used in the orbifold projection. For — )■ 00, these equations become a 
closed set of loop equations, and, in the strong coupling/large mass limit, they can be solved 
by an iterative algorithm, which generates the large- lattice strong-coupling expansion for the 
corresponding observables. 

Note that the construction of an orbifold duality, that we outlined above, provides an explicit 
example of how the spacetime-dependent degrees of freedom of the daughter theory can be obtained 
from a projection of internal degrees of freedom in the parent theory. 

The analysis of ref. |652j was later extended beyond the strong-coupling domain in ref. [653], 



which used the large- A^ coherent-state methods described in ref. [27] (and reviewed in subsec. 4.2) 



^^As discussed in ref. [652 , operators involving adjoint scalar or fermion fields can also be expressed through 
closed Wilson loops (without clumsy insertions of matter field decorations along loops made of only of gauge link 
variables), by generalizing the lattice to include one further dimension for each matter species. 
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to demonstrate that the symmetry conditions in the parent and daughter theories are both nec- 
essary and sufficient for (i.e., completely characterize) the validity of large- orbifold dualities. 

Using these tools, the validity of the orientifold equivalence, relating the planar limits of 
QCD with an adjoint Majorana fermion (i.e., M = 1 supersymmetric QCD) and of QCD with a 
Dirac fermion in the two-index antisymmetric representation j617[ I618[ I619[ l62Uj . was discussed 
in ref. |616] . In particular, it was shown that the equivalence between these two theories can be 
interpreted in terms of a "daughter-daughter" orbifold equivalence: both theories can be obtained 
from a common parent theory, namely S0(2A^) M = 1 supersymmetric QCD, by applying orbifold 
projections based on two different Z2 groups. Projecting out the degrees of freedom which are 
not invariant under J = 102 yields the M = \ supersymmetric U(A^) theory as a daughter, 

whereas a projection using a "graded" variant of the same operator, J(— 1)^, leads to the U(A^) 
theory coupled to a fermion in the two-index antisymmetric representation of the gauge group. 
Furthermore, it was also pointed out that the validity or invalidity of the planar orientifold equiv- 
alence (in the neutral sectors of both daughter theories), being intimately related to realization 
of the discrete symmetries used in the projections, depends crucially on the dynamics, and, as 
an example, it was shown that the orientifold planar equivalence may fail, when the theories are 
compactified on a smal^ spatial torus: this regime can be reliably studied with weak-coupling 
expansion methods, which show that charge conjugation gets spontaneously broken — and, as a 
consequence, the orientifold equivalence fails. 

As mentioned above, large- EK volume independence can also be interpreted in terms of an 
orbifold equivalence: as discussed in detail in ref. ^llj, in this case the parent theory (defined, 
say, on a D-dimensional hypertorus of linear size L) can be mapped to the daughter (defined on 
a smaller hypertorus of size L', with L an integer multiple of L', i.e. L = kL'), by projecting out 
the degrees of freedom that are not invariant under the translations by integer multiples of 
L' in each direction — or, equivalently, removing the Fourier components that are not quantized in 
units of 2tt /V . This leads to a correspondence between theories defined in two different volumes, 
and, for a parent theory on a lattice of spacing a, yields the EK model as the daughter theory, 
if L' = a. Conversely, it is also possible to map a theory defined in a smaller volume to one on 
a larger volume: this can be done by associating some of the internal degrees of freedom of the 
parent theory to the spacetime degrees of freedom of the daughter. For example, the U(A^) EK 
model with A^ = k'^n can be mapped to ordinary U(n) Yang- Mills theory on a k'^ lattice, following 
the construction that we outlined in the discussion of orbifold projections earlier in this subsection, 
which amounts to imposing a set of constraints of the form: 

Ua = exp{2ni6ai3/ k)ji3Uaj^^ , (133) 



with the 7q's defined analogously to eq. (131). The action (as well as the other neutral-sector 
observables) of the EK model is then mapped to the action of the lattice gauge theory in the 
enlarged volume (up to a trivial rescaling of the couplings), and the two theories are equivalent, as 
long as center symmetry is not spontaneously broken; furthermore, the inclusion of matter fields 
in various representations is also possible (see ref. [611] for details). 

An interesting application of combined orbifold equivalences consists in the possibility of study- 
ing the large-A^ limit of QCD with fermions in a two-index symmetric or antisymmetric repre- 
sentation in a large volume, by ffist mapping this theory to QCD with adjoint fermions via the 
orientifold equivalence (which holds in a large volume), and then studying the latter theory in a 



'^^Here, "small" means "small with respect to the inverse of the characteristic dynamical scale at which the theory 
becomes strongly interacting", i.e. small as compared to 1/Aqc£,. 
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small volume (possibly reducing it to a matrix model, or to a single-site lattice model), using the 
fact that, for adjoint QCD, the orbifold equivalence relating the theory in different volumes hold 
all the way down to infinitesimally small system sizes. 

Orbifold equivalences relating large- iV theories defined in systems of different sizes also have 
interesting physical implications for the phase diagram of these theories: in particular, volume- 
reducing or volume-expanding projections, that change the size of a Euclidean gauge system only 
along one direction, increase or reduce the temperature of the system by a factor k. If center 
symmetry is not spontaneously broken (i.e. if the system is in its confining phase), then the 
correspondence of physical observables in the neutral sectors of the parent and daughter theory 
implies that, at large A^, such observables are temperature-independent in the confining phase, 
Another potentially interesting application of orbifold equivalences for the study of the phase 
diagram of QCD has been proposed in ref. |663] and discussed in refs. |664l 16651 16661 16671 16681 
I669[ 1670] . and is based on the observation that the large- equivalence of theories based on 
orthogonal and unitary gauge groups may allow one to get information about the latter in the 
finite baryon density regime (in which lattice simulations are hindered by a severe computational 
sign problem |671t I672[ 1673] ). by performing numerical simulations of the former, for which the 
sign problem is absent. Recently, the problem of baryons in the context of orbifold equivalences 
has also been discussed in ref. |674] . 

Finally, one further important application of orbifold constructions is in the context of lattice 
supersymmetry |675[ I676[ I677i I678i I679i I680i I68H I682j — for further details, we refer the reader to 
the review |683] and to the references therein. 



5. Large- A/^ results from lattice simulations 



In this section, we review the main large- A^ results for physical observables from lattice simula- 
tions. Unless otherwise stated, these are obtained from simulations on large lattices (i.e., without 
using the properties of large- A^ volume reduction). First, in subsection 5.1, we discuss the nu- 



merical results relevant for four spacetime dimensions; then, we briefly review the results that 
have been obtained for large- A^ theories in three spacetime dimensions in subsec. 5.2 Finally, 



we conclude this section with a summary of results from studies in two spacetime dimensions in 
subsec. 



15.31 



5.1. Results in four spacetime dimensions 



As mentioned in subsection 2^, many phenomenological implications derived from the large- A^ 
counting rules are based on the assumption that QCD be a confining theory in the large- A^ limit. 
Since confinement is a non-perturbative phenomenon, the validity of this assumption should be 
assessed with non-perturbative methods. 

The first studies addressing this issue were reported in refs. 1684] . comparing the SU(2), 
SU(3) and SU(4) Yang-Mills theories. A possible way to assess whether these theories are confining 
consists in computing the two-point correlation function of zero-transverse-momentum, purely 



gluonic, spatial string operators winding around the lattice (torelons) along a periodic direction 



of size L, and studying its behavior as a function of the torelon separation r: at large r, the 



^^Related topics have been discussed in the recent paper |662j . in which the temperature independence of a 
three-dimensional Yang-Mills theory with adjoint fermions was used to map it to a two-dimensional gauge theory, 
which was then studied using non-Abelian bosonization methods. 

^^The formation of string-like objects in gauge theories is an idea that dates back to almost forty years ago [685] , 
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correlator was found to decay exponentially, oc exp[— m(L)r], where m(L) represents the mass 
of the lightest torelon. For long torelons, m{L) tends to become linear in L (indicating that the 
torelon is characterized by a non- vanishing linear energy density, or string tension, a), so the mass 



of an infinitely long torelon diverges, implying confinement,^^ The results from refs. [311 I684j 
showed that confinement persists in all the theories studied. 

Another important problem, that was investigated in the first numerical studies of large- 
gauge theories, is the relation between the scale and the gauge coupling. A key statement in 
large- A^ arguments is that the limit N = oc exists if we keep g^N fixed. Equivalently, we might 
say that, perturbatively, the dynamically generated scale Aqcd is the same for all A^, at fixed g'^N. 
To define non-perturbatively the large- A^ limit, it is useful to look at the dependence of Aqcd (or of 
a related, dimensionful, dynamically generated scale, such as the string tension a) on g'^N in the 
non-perturbative regime. By determining non-perturbatively how the lattice spacing a varies, as 
a function of the bare lattice gauge coupling, ref. |499j showed that the square of the bare lattice 
coupling scales roughly like 1/A^ {'t Hooft scaling), as expected from perturbative arguments. 

In fact, even stronger evidence for 't Hooft scaling can be obtained, by using a "mean-field 
improved" definition of the bare coupling |687[ I688j , obtained by dividing the square of the bare 
coupling by the average plaquette (normalized to 1 in the weak-coupling limit). This improved 
definition of the lattice coupling comes from the observation that, expanding the relation between 
lattice fields and continuum fields as: 

U^{x) = exp [iagoA^{x)] = 1 + iagoA^{x) - ^a^glAl^x) + ... (134) 

reveals that, in addition to the terms reproducing the continuum Yang-Mills action, the lattice 
action also includes infinitely many vertices with multiple gluons, which are pure artifacts of the 
lattice discretization. At the classical level, these terms are irrelevant in the continuum limit, 
because they are multiplied by higher powers of the lattice spacing a. At the quantum level, 
however, the tadpole diagrams arising from the contraction of fields in these terms turn out to 
be proportional to (a power of) the lattice cut-off n/a, so that, in fact, in the continuum limit 
the artifacts are only suppressed by powers of go. To alleviate this problem, one can introduce a 
tadpole- or mean-field-improved coupling, which is obtained by integrating out the high-frequency 
modes of the gauge fields, and reabsorbing the resulting constant into a redefinition of the cou- 
pling |687tl688] . In particular, a convenient definition of the rescaling factor for the lattice coupling 
/3 is given by the average plaquette so that the corresponding improved coupling can be inter- 
preted as the effective coupling experienced by the dynamical variables in a uniform background 
field. The scaling of the mean-field improved bare lattice coupling with A^ (at a fixed physical 
scale) was studied in refs. |499[ 1689] and later in more detail in ref. |61U] . from which the plot 
shown in fig. [6] is taken: the figure shows the mean-field improved 't Hooft coupling as a function 
of the lattice spacing, determined non-perturbatively in units of the inverse of the string tension 
square root: the collapse of results obtained from different SU(A^) Yang-Mills theories (for A^ in 
the range from 2 to 8) confirms that, indeed, a physically meaningful large- A^ limit is the one at 
fixed 't Hooft coupling — also at the non-perturbative level. 



^^Note that, as discussed, e.g., in ref. |686) . this criterion for estabhshing the confining behavior of a theory is 
conceptuahy well-defined, and free from the potential ambiguities between confinement and screening, that may 
hinder alternative definitions of the phenomenon. In addition, it also provides a viable operational way to study 
the phenomenon in lattice simulations. 

^^An alternative definition of the rescaling factor is based on the average value of the link variables in Landau 
gauge. However, the computation of the average plaquette is numerically much easier. 
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Figure 6: Relation between the mean- field improved 't Hooft coupling and the lattice spacing (in units of l/\/a), 
as determined non-perturbatively in ref. |610) . in Yang-Mills theories with a different number of colors, from 2 to 
8. 

While the bare (possibly improved) lattice gauge coupling can be interpreted as a physical 
coupling at distance scales of the order of the lattice spacing, it is also interesting to investigate 
how the coupling runs at lower energies — in particular, because, by definition, lattice gauge theory 
describes continuum physics for phenomena whose characteristic energy scales are well below the 
intrinsic scale set by the inverse of the lattice spacing. A possible scheme choice is provided by 
the Schrodinger functional (SF) |690[ 16911 1692] : in a nutshell, the idea is to extract the running 
coupling g{L) at a given momentum scale proportional to 1/L, by studying the effective action of 
the gauge system with fixed boundary conditions at the opposite ends of a direction of size L (and 
periodic boundary conditions along the other directions, which are taken to be of large extent). 
Denoting the field configurations on the two boundaries as Q and C/, the SF is defined as the 
probability amplitude for evolution from the state to the state Cf. In the Euclidean setting, 
this can be written as: 

Z[Ci,Cf]= f D[/exp(-5) = exp(-r) (135) 

JC,,Cf 

where the path integral is done at fixed boundary conditions, and F = F[Ci,C/] denotes the 
corresponding effective action. For and C/, one takes the lattice fields to satisfy inhomogeneous 
Dirichlet boundary conditions, depending on a dimensionless real parameter vj, so that a physical 
running coupling at the length scale L can be obtained, by comparing the derivative of the classical 
action and of the full effective action, with respect to the t] parameter, 



This allows one to trade the explicit computation of the effective action, for the computation of expectation 

59 




Figure 7: Left panel: The running coupling in the SF scheme for SU(4) Yang-Mills theory, as determined in 
ref. [693], assuming ^ ~ 420 MeV. The black and blue dashed curves are the one- and two-loop perturbative 
predictions. Right panel: Dependence of the AqcdZ-v/ct ratio (after conversion to the MS scheme) on the number 
of colors, as obtained from computations of the physical coupling in the SF scheme in SU(2) 694 , SU(3) [69^ 
and SU(4) Yang-Mills theory |693| . The plot, taken from ref. |693| . also shows the comparison with the large- iV 
extrapolation (solid straight line) obtained from the study of the mean-field improved lattice coupling in ref. 1610] ; 
the region bounded by the dash-dotted lines corresponds to 68.3% confidence level. 



In ref. |693j . this approach was used to compute the running couphng in the SF scheme in 
SU(4) Yang-Mills theory (and to discuss a large- extrapolation, by comparison with analogous 
results for the theories with two |694j and three colors [695j): the simulation results showed that 
the running coupling is in very good agreement with the two-loop perturbative /3-function down 
to momentum scales of the order of a few hundreds MeV, and that the ratio of the Aqcd scale (in 
the modified minimal subtraction scheme) over the square root of the string tension has a smooth 
dependence on with the value in SU(3) theory already very close to the extrapolated large- 

limit (see fig. [T]). 

Very recently, a similar type of comparison has been performed for SU(iV) theories with 
= 2, 3 and 4 colors and with two flavors of dynamical fermions in the two-index symmetric 
representation of the gauge group |696j . The motivation for this study comes from the sugges- 
tion |697[ I698[ I699j that these theories are among the potentially interesting candidates for walking 
technicolor models of dynamical electro-weak symmetry breaking, which recently have been the 
subject of many lattice studies (see, e.g., refs. |700l 17011 17021 1703] for reviews). For this reason, 
the authors of ref. |696] focused on the investigation of the /3-function, and on the behavior of the 
mass anomalous dimension 7^ (which relates the pseudoscalar renormalization Zpg to the length 
scale L via: Z^^iV) oc L"'*'™) as a function of the coupling, finding remarkable similarities between 
the theories with two, three and four colors, as shown in fig. [8j 

Coming back to the problem of the heavy quark-antiquark potential in large- pure Yang-Mills 
theory, the first lattice studies [311 14991 (6841 [?04J also showed that, besides the term linear in L, the 
masses of long torelons also include a correction proportional to the inverse of their length, with 
a coefficient which can be evaluated analytically, which is known as the Liischer term |705[ I706j 
(see also ref. |707] ): 

m(L) = ctL - ^ + . . . (136) 
where the ellipsis denotes subleading terms, suppressed by higher powers of 1/L. The left panel of 



values of the operators obtained by deriving the lattice action with respect to 77. 
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Figure 8: Mass anomalous dimension in SU(7V) gauge theory with N = 2, 3 and 4 colors and two flavors of 
dynamical fermions in the two-index symmetric representation, as a function of the renormalized 't Hooft coupling 
in the SF scheme, from ref. |696) . The dashed green line denotes the leading-order perturbative prediction in 
N ^ oo limit. 



fig. [9] shows a lattice calculation of torelon masses in SU(6) gauge theory. The Ansatz in eq. (136) 
describes the numerical data down to lengths of the order of 3/y/a. 

As discussed, for example, in refs. |709[ 1710] . the Liischer term can be interpreted as a Casimir 
effect (namely: as a quantum, finite-size effect), which arises from massless fluctuations of the 
flux tube along the transverse directions, and signals that the low-energy excitations of the flux 

For an overview of recent studies 
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tube can be described in terms of a vibrating, bosonic string, 
of this subject, see, e.g., the slides of the presentations at the "Conflning Flux Tubes and Strings" 
workshop held in ECT*, Trento, Italy, in 2010 [7T3] . 

While the existence of string-like excitations is a generic feature of all conflning gauge the- 
ories [m inSl inSl inZl EHI inni HM ESU, the possibility that a particularly simple effec- 



^^These massless fluctuations can be interpreted as the Nambu-Goldstone modes associated with the spontaneous 
breakdown of Lorentz-Poincare symmetry, due to the formation of the flux tube. 

^^As discussed in refs. |7111I712] . it is interesting to note that an effective string model for the gauge theory also 
suggests the existence of a deconfinement transition at finite temperature. 
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Figure 9: Left panel: Torelon mass, in lattice units, as a function of the length of the loop (in units of the string 
tension). Note that the dashed line is not a best fit to the data, b ut ha s been obtained by extracting the string 
tension at the point for which = 5 and inserting this value in eq. ( 136 ). The figure is taken from ref. Right 
panel: Ratio of the tension of the string with A^-ality k = 2 over the fundamental string tension. The Pisa group 
data are taken from ref. [689] . while Oxford data are taken from ref. |708j . The long dashed line is a continuum 
extrapolation of the Oxford data, obtained on symmetric lattices (see ref. |708j for technical details). 



tive string model for the low-energy dynamics of large- QCD could become exact is suggested 
by theoretical arguments |722j . and has been studied in a number of lattice works — see, e.g., 
refs. laH 11991 MM EOl M UM [121] and references therein. 

Later, the lattice analysis of the confining potential between static probe sources in SU(A^) 
Yang-Mills theories was extended to color sources in higher representations, and to excited states. 
In particular, a series of works |689[ l7U4t I7U81 I725[ I726[ 1727] investigated the behavior of the 



potential between static sources in representations of various A^-alities (i.e., with different trans- 
formation properties with respect to the Zat center of the gauge group). Part of the motivation 
to investigate the confining potential for sources in representations of different A^-ality stems from 
the observation that, in supersymmetric SU(A^) gauge theories, the associated string tensions ak 
(for stable strings) are related to each other by the sine-formula |728[ 1729^ 1730^ 1731] : 



(Tfc sin(A;7r/A^) 
(Ti sin(7r/A^) 



(137) 



where ai denotes the fundamental string tension. This relation can be compared with the one 
that can be obtained from Casimir scaling |732] . which, for the totally antisymmetric irreducible 
representation of A^-ality k, implies: 



CTfc _ k{N - k) 



(138) 



In the large- limit at fixed fc, both eq. ( 137 ) and eq. ( 138 ) are consistent with the expectation from 



factorization: Okjox — )■ k. However, lattice results show that at finite N the energy of A;-strings 
turns out to be lower than k times the energy of a fundamental string, and, in particular, appear to 



favor sine scaling |689ll725il726j or to lie in between the predictions of eq. (137) and eq. (138) |704l 
1708] . A summary of the situation is presented in the right panel of fig. |9j the different conclusion 
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is not due to the numerical data, which are compatible, but to the extrapolation to the continuum 
limit. To resolve the apparent discrepancy, more simulations closer to the continuum limit need 
to be performed. On the other hand, for different irreducible representations of the same A^-ality, 
showed that the lattice results are broadly consistent with Casimir scaling, for SU(3) Yang-Mills 
theory [TS^ as well as for iV > 3 [TUS] . 

Recently, a lot of work has been devoted to the study of excited string states in large- Yang- 
Mills theories (see ref. [21] and references therein for a thorough discussion of the topic): this is 
important for understanding the nature of the effective action describing the low-energy dynamics 
of confining flux tubes. A large number of lattice studies (recently reviewed by Mykkanen in 
ref. [710j ) show that the main features of the dependence of the torelon mass on its length, or, 
equivalently, of the ground-state quark- ant iquark potential on the distance between the color 
sources, are consistent with the hypothesis that the flux tube fluctuates like a bosonic string 
described by the Nambu-Goto action |533[ 1534] , which is simply proportional to the area of the 
string world-sheet Due to the simplicity of the string action, the Nambu-Goto spectrum can be 
computed exactly: for an open string of length L with fixed ends, it reads |736t 1737] : 




En{L)=aL\l + — \n-——\, n G N (139) 



(where D denotes the number of spacetime dimensions, and a is the string tension). Expanding 
the lowest-lying energy around the cjL^ — >■ oo limit yields (besides the classical, linear term aL) 
the Liischer term for the open string: 

E,{L) =aL- ^{D - 2) + 0(1-^). (140) 
Similarly, for a closed string, the Nambu-Goto spectrum reads (see, e.g., ref. |738] ): 



; 47r / D-2\ 47r2(n/-n^)2 , , 

E„„.,(L) = aL^/1 + — i^m + nr- -^j + ^J,^, (141) 

(where ni and are non-negative integers denoting the number of left- and right-moving modes, 
respectively), which, when expanded around the long-string limit in four spacetime dimensions. 



reproduces eq. (136) for the ground state. 

During the last decade, various theoretical works tried to determine the form of the subleading 
(in an expansion in terms of the inverse of the string length) terms which the effective string 
action describing a confining flux tube should include, using arguments related to open-closed 
string duality |739] and to Lorentz-Poincare symmetry and its non-linear realization |740[ 17411 
I742[ I743[ I744[ I745[ 1746^ I747j (a different approach was proposed by Polchinski and Strominger in 
ref. |748j — see refs. |749[ I75UI 1751] for a discussion of the relation between the two approaches). 
For a string of length L in four spacetime dimensions, this shows that the string energy levels 
deviate from the Nambu-Goto spectrum at order 1/L^ for an open string |742] . whereas for a 
closed string the deviations occur at order for excited states, but only at order 1/L^ for the 
ground state |742j . 



^•^It is well-known that (unless the number of spacetime dimensions is D — 26) the Nambu-Goto action is not 
a consistent string action at the quantum level, due, for example, to the Weyl anomaly, and to the existence of 
tachyonic states [7341 1735] , Nevertheless, these problems are not relevant for the case of an effective, low-energy 
model, as discussed in the present context. 
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From the numerical point of view, the spectrum of closed strings in four spacetime dimensions 
was studied in ref. |752] . using a sophisticated variational technique involving about 700 operators, 
in SU(3), SU(5) and SU(6) Yang-Mills theories. The results showed general agreement with the 
Nambu-Goto spectrum, even down to surprisingly short values of L, of the order of the inverse 
square root of the string tension — although some discrepancies were also reported. An example of 
these results is shown in fig. 10 As discussed in ref. |743] . the precision of these very challenging 
lattice computations is still insufficient to clearly distinguish deviations from the Nambu-Goto 
spectrum in the range of L values within the convergence domain of the expansion in inverse 
powers of L. 

Having ascertained that non-Abelian gauge theories are confining in the 't Hooft limit, i.e., 
that a linearly rising potential develops between a static quark and antiquark, the next problem to 
be addressed non-perturbatively, via lattice simulations, is determining the values of the hadron 
masses at large N . In the hterature, there exist computations of the masses of glueballs [311 1499i 
I708[ 1753(1754] . of mesons |755[ I756( 17571 1758] and baryons |759] — all in the quenched approximation. 

In general, the dependence of glueball masses on turns out to be very smooth: typically, for 
each physical state, the lattice results obtained in Yang-Mills theories at different values of A^ can 
be nicely fitted to a constant plus a term linear in l/N"^ (as expected theoretically); within the 
precision of the numerical results, this holds all the way down to A^ = 3 — or even N = 2, in some 
cases. Results in the continuum limit for the lowest-lying states were presented in refs. |499[ 1708] . 
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Figure 11: Left panel: Extrapolation to the continuum limit of the lowest-lying glueball spectrum in SU(4) gauge 
theory. Right panel: Extrapolation of the lowest-lying glueball spectrum for N — >■ oo, using only the leading 
finite- correction. Both figures are plots of results published in rcf. |708j . 

In particular, the calculation of ref. |7U8j leads to the following results: 

mn++ , , 2.1(1.1) 

= 3.28(8) + , (142) 



mo++. 2.7(2.0) 

= 5.93(17) , (143) 



^-4.78(14) + ^. (144) 



The quality of the data and the extrapolation to the continuum limit is shown in fig. 11 



The most recent computation of the glueball spectrum (performed at fixed lattice spacing) was 



reported in ref. |754] . and is summarized in fig. 12, showing the masses of the ground state (and 



of some excited states) in the different channels associated with the irreducible representations of 



the cubic group and distinguishing states that are "even" and "odd" under parity and/or charge 
conjugation. For comparison, the figure also shows the mass of the ground state and of the first 
radial excitation in the scalar channel, and the ground state in the tensor channel, taken from an 
earlier work |7U8] . In addition, a lattice study of the relation between mass and spin of the various 
glueballs in SU(3) and SU(8) Yang-Mills theories has also been carried out, in ref. |753] . 

Similarly, a smooth dependence on N has also been observed for the meson spectrum |755l 17561 



1758] : this is shown in fig. 13, reporting preliminary results from ref. |758j . in which the masses of 
various meson states in the chiral limit are shown — together with the pion and p decay constants 
(rescaled by \/N). Symbols of different colors refer to different values of A^, and the band denotes 
the extrapolation to the 't Hooft limit. 

These results are consistent with the studies previously presented in refs. |755[ 1756] . which 
confirmed that the values of the pion and p masses in the large- limit are close to those in the 
real world: for example, combining original results with those of |755j . ref. |756] reported: 

lim ^ = 1.79(5) (145) 



The computation is done on a hypercubic lattice, at finite lattice spacing. 
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Figure 12: The spectrum of glueballs in the large- iV Hmit of SU(A^) Yang- Mills theory, as determined in ref. [754] . 
The plot shows the results (at a fixed, finite value of the lattice spacing) for the ground state and some of the 
excitations, in various channels corresponding to different irreducible representations of the cubic group, and 
different values of the parity and charge-conjugation quantum numbers. The figure also shows a comparison with 
the results obtained in ref. |708j for the ground state and for the first excited J^^ = 0++ glueball, and for the 2++ 
ground state. 



as an estimate for the ratio in the continuum hmit. This resuh is consistent with the value 
(approximately equal to 1.75) obtained from the experimentally measured mass of the p meson, 
and from a reasonable phenomenological estimate of a. Furthermore, the dependence of the p 
mass on the pion mass appears to be approximately quantitatively consistent with the holographic 
computations discussed in ref. |158] . 

We note that, by contrast, ref. |757] reported very different results, indicating that the mass 
of the p meson in the large- limit could be about twice as large as in the real world. A possible 
explanation of the discrepancy between the results of ref. |757j and those of refs. |755t 1756] in 
terms of an unexpectedly large finite- iV correction (the simulations presented in refs. |755t I756j 
were limited to < 6, whereas those in ref. |757] were done for A^ = 17 and 19) has now been 
ruled out by the results of the most recent study |758j . which also included simulations for A^ up 
to 17, and found consistency with refs. |755l 1756] . Perhaps, a more likely interpretation of the 
difference between the results in ref. |757] and those in the other three works |755l 1756] 1758] could 
be the one already suggested in ref. [21], which noted that the computation presented in ref. |757j 
was based on the evaluation of quark propagators in momentum space, and was carried out for 
only a few small momenta, so it may have suffered from contamination from some excited state. 



^^The value of the string tension a can be estimated from the analysis of Regge trajectories in experimentally 
observed meson states, and from studies of charmonium and bottomonium spectra |7601 176T1 1762] . which suggest 
values in the range between (400 MeV)^ and (450 MeV)^— see refs. [7531 [TM] for a discussion. 
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Figure 13: Preliminary results from a quenched computation of the mesonic spectrum at large N |758j : the figure 
shows the masses (in units of the square root of the string tension) of various meson states, and the decay constants 
of the pion and the p meson (divided by ^/N), as obtained from simulations for different numbers of colors, in the 
chiral limit. The values extrapolated to — ^ oo are shown by the horizontal bands. 



Finally, a recent work |759] presented a quenched computation of the baryonic spectrum for 
large values of N, up to 7; this study was restricted to odd values of N — so that the baryons are 
fermions, like in the real world. Once again, the results confirmed the theoretical expectations for 
the dependence of physical quantities on the number of colors (including, in particular, baryon 
masses approximately linearly increasing with N), up to numerically modest relative finite- 
corrections. In addition, this study also revealed that states of larger spin are heavier, and, more 
precisely, the dependence of the mass of a baryon on its spin is consistent with a rotor spectrum. 



eq. (15), as predicted in refs. |161t I162j . This can be seen in fig. 14, showing the linear relation 
between various mass splittings for baryon states of different spin, on each other. 

The topological properties in large-iV Yang-Mills theories have been studied numerically in 
various works |499l 17651 17661 17671 1768] . and are summarized in the review In particular, 

ref. [499j considered SU(A^) Yang-Mills theory with iV = 2, 3, 4 and 5 colors, and found that, in 
agreement with the theoretical expectation |260J, the number density d of instantons in the Yang- 



Mills vacuum is exponentially suppressed when increases — see eq. (25) — and that the density 
of instantons of small size p at fixed scales compatibly with: 

d{r) oc p^'^-^ (146) 

as one can predict using perturbation theory (which should become reliable for p — )■ 0). In 
addition, ref. |499] also found that the topological susceptibility tends to a non-vanishing value 
for — 7- oo: 

^ = 0.376(20) + (147) 
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Figure 14: A recent quenched computation [759] of the baryon spectrum at large N reported evidence for a rotor- 
type spectrum in the masses of different states as a function of their spin, as predicted in refs. 161,1162 . The plot, 
taken from ref. |759j . shows differences between the masses of states of various spin, plotted against each other, in 
QCD with N = 5 (left panel) and N = 7 colors (right panel). 



Fig. 15 shows the resuhs from two subsequent, similar studies \767\ I768J . which found resuhs 
consistent with each other and with ref. |499j . In particular, ref. |767] reported simulations for 
= 3, 4 and 6 colors, yielding: 

^ = 0.0221(14) + (148) 
while the authors of ref. |768j studied the theories with = 2, 3, 4, 6 and 8 colors, obtaining: 

= 0.382(7) + —j^^ with: x /d-o.f. = 1.3, (149) 

using the real values for the topological charge computed using a twisted plaquette operator (open 



symbols in fig. 15) — and very similar results with a definition of the topological charge truncated 



to integer values (filled symbols in fig. 15). To summarize, these studies indicate that the large- 
limit of the topological susceptibility is non-vanishing, and around (170 MeV)^, rather close to its 
value in the SU(3) theory (which is about (180 MeV)^). 

The relation between the topological charge of the gauge fields (evaluated on "smoothed" gluon 
field configurations) and the chirality of the quark fields, as probed by the low-lying eigenmodes 
of the overlap Dirac operator |464] , was investigated in ref. |766j , which found that the numerical 
results for these two quantities become consistent when N becomes large, and that the low-lying 
eigenmodes of the Dirac operator tend to lose their chirality properties and to become less and 
less localized for — > 00, suggesting that instantons do not survive and are not responsible for 
chiral symmetry breaking in this limit. 



As discussed in subsection |2.4[ the topological susceptibility is closely related to the dependence 
of non-Abelian gauge theory on the 6'-term, whose lattice investigation in the large- A^ limit has 
been carried out in various works (including, e.g., refs. |767l 1769] ). and is reviewed in ref. [32] — 
see also refs. |770[ 17711 1772] and references therein for analogous studies in the SU(3) theory. It 
should be pointed out that the formulation of the theory with a topological ^-term on a Euclidean 
lattice suffers from a sign problem which prevents its direct numerical simulation, due to the 
complex nature of the action. A possible way to circumvent this problem consists in studying 
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Figure 15: Left panel: Results for the topological susceptibility (in units of the square of the string tension), from 
the simulations performed in ref. |767j . and their extrapolation to the large- limit. Right panel: The fourth 
root of the topological susceptibility (in units of the square root of the string tension) in SU(2), SU(3), SU(4), 
SU(6) and SU(8) Yang-Mills theory at zero temperature, as determined in ref. |768j using a real- (open symbols) 
or integer-valued (filled symbols) definition of the lattice topological charge operator. The corresponding quadratic 
fits, as a function of l/N'^, are also shown. 



the Taylor expansions of the physical quantities of interest around ^ = 0. In particular, following 
this strategy, the authors of ref. |769j studied the 6'-dependence of the string tension and of the 
mass gap in SU(3), SU(4) and SU(6) Yang-Mills theory, finding results suggesting that such 
dependence vanishes in the 't Hooft limit, in agreement with the theoretical expectation that the 
relevant parameter at large should be 6/N. 

In addition to these studies at zero temperature, in the literature there exist a number of 
works about large- iV lattice gauge theories at finite temperature |773t I774[ I775[ 1776^ I777[ I778[ 17791 
I780[ I78H I782[ 17831 1784[ 1785] . Typically, in these simulations the temperature (which, in Euclidean 
thermal field theory, is just the inverse of the system size along a compactified direction: T = 1/L) 
is varied by changing the value of the lattice spacing a, which, in turn, means changing the param- 
eter 13 = 2N/gl = 2N'^/Xo in front of the Wilson gauge action, i.e. the inverse of the bare gauge 



coupling. As mentioned above, in order to probe the region of couplings analytically connected 



to the continuum limit, the numerical simulations have to be run at sufficiently small values of 



^^A different approach, based on varying the temperature by changing the number of lattice sites along the 
compactified direction, is also possible |786| . however it has the practical disadvantage that it does not allow one 
to vary the temperature continuously. 
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the bare lattice 't Hooft coupling Aq. However, it is remarkable that interesting implications for 
SU(A^) thermodynamics can also be obtained from lattice strong-coupling techniques |356] . 

In four spacetime dimensions, the simulation results show that all SU(A^) Yang-Mills theories 
have a deconfinement transition at a finite critical temperature, which — when expressed in physical 
units — is in the range between 250 and 300 MeV. The transition turns out to be of second order for 
the SU(2) gauge theory [78 7[ 17881 1789j , while it is a weak first-order one in the SU(3) theory |790l 
1791] and a stronger first-order one for all > 4 |773l 17741 17751 1782] . By comparison, note that 
in QCD with Uf = 2 + 1 fiavors of dynamical quarks, for physical values of the quark masses 
the deconfinement transition is actually a crossover which takes place in a temperature range 



between 150 and 170 MeV [ 792^ 1793^ 1794^ [7^J , and appears to be compatible with experimental 
evidence |804] . Further details about the lattice investigation of the QCD equation of state can 
be found in refs. |805l 1673] . 

In large-A^ pure Yang-Mills theory, the deconfinement critical temperature Tc can be deter- 
mined unambiguously, and the most recent, highly precise results show a mild dependence of this 
quantity on the number of colors (see fig. 16, taken from ref. |785] ). which can be summarized by 
the relation: 

^= = 0.5949(17) + with: xVd.o.f. = 1.18, (150) 



obtained in ref. |785] using the results for all values of from 2 to 8. 

As mentioned above, the latent heat Lh associated with the deconfinement transition is finite 
for all A^ > 3, and scales proportionally to A^^ in the large- A^ limit |775] 1781] : in particular, 
ref. |781] reported: 

in full agreement with the result obtained in ref. |775] for the same quantity, which reads: 
0.766(40) - 0M{1.Q0)/N\ 

In addition to the finiteness of the latent heat, the first-order nature of the transition for A^ > 3 
is also related to the finiteness of the surface tension (i.e., energy per unit area) •jw associated with 
interfaces separating different center domains, which, in real-world QCD, might be of experimental 
interest |806] (although the relevance of these objects for physics in Minkowski spacetime has 
been questioned [332J). In particular, the probability for formation of a "bubble" of surface A 
separating two different phases near the critical point is proportional to exp {—'jw-A/T) — so that 
tunneling events in a finite volume are exponentially suppressed by the finiteness of 'jw Note 
that, from the point of view of numerical simulations, this makes it practically necessary to study 
the deconfinement transition at large A^ on lattices of relatively small volume, in order to have a 
sufficient number of tunnelings during the Monte Carlo history. On the other hand, however, the 
very first-order nature of the transition also makes it sufficient to consider small lattices, since it 
suppresses finite-volume effects (for a more detailed discussion of finite-volume effects in large- A^ 



^ In QCD with dynamical quarks of finite mass, there exists no bona fide order parameter for the finite- 
temperature transition, since center symmetry is exphcitly broken by the existence of quarks, and chiral symmetry 
is exphcitly broken by their finite mass. 

^^Previous studies |7961 17971 17981 1799] reported somewhat larger values (up to around 190 MeV) for the crossover 
pseudo-critical temperature, which sparked off some heated but constructive debate in the community |8001 1801] . 
However, it seems that this disagreement can be explained in terms of the uncertainties related to the extrapolation 
to physical quark masses and to the continuum limit |802] . and, more recently, the same collaboration presented 
a new study [795] . based on a highly improved lattice formulation for the quark fields |803] . in which they found 
results compatible with those obtained by the other group [7921 1793] . 
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Figure 16: Dependence of the deconfinement critical temperature (in units of the square root of the string tension) 
on the number of colors in SU(A^) Yang-Mills theory. This figure, taken from ref. |785j also shows the interpolation 
of the numerical results described by eq. (150 1. The inset shows a zoomed view of the simulation results. 



gauge theories at finite temperature, see also ref. |780] ). In tlie limit of very high temperatures, 
the surface tension separating two different domains labelled by Polyakov loops with values 1 and 
exp(27riA;/iV) in SU(iV) Yang-Mills theory can be computed perturbatively [8071 EOHl |809l ISTnl 
I81H 1812] — see also ref. |813] — with the leading-order result: 



7^^ 



4k{N - k)7r^T^ 



3V3A 



[1 + 0(A)]. 



(152) 



In particular, this shows that the largest of these interface tensions (for k = \_N/2\) is 0{N^) 
in the 't Hooft limit. As discussed in ref. |774j . for temperatures close to it is reasonable 
to expect that the tunneling from the center sector with Polyakov loop equal to 1 and the one 
with k = \_N/2\ goes through an intermediate tunneling to the center-symmetric phase, with a 
probability related to the interface tension 'j^'^ between the confining and deconfined phases, so 
that: 7;^'^ = 7|^'=/2 = 0{N^). The A^-dependence of 7^^'^ was studied in ref. [773], with the 
result: 



0.0138(3) 



0.104(3) 

iV2 ' 



with: xVd-o.f. = 2.7. 



(153) 



As discussed in ref. |810] . the surface tension 7^!^'^ associated with walls separating different 
center domains in the deconfined phase is also trivially related to the "dual" string tension a 



characterizing the area-law decay of large 't Hooft loops |248j in the deconfined phase: aT = 7^' 
The behavior of 't Hooft loops in large- A^ Yang-Mills theories at finite temperature was studied 
numerically in refs. |776] 1779] . which computed the values of a corresponding to 't Hooft loops 
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Figure 17: Left panel: The bare Polyakov loops in twelve different irreducible representations of SU(6) Yang- 
Mills theory, when their free energies are rescaled dividing them by the quadratic Casimir of the corresponding 
representation, fall on the same curve; the figure shows the results obtained from simulations on a lattice of spacing 
a = 1/(5T), plotted as a function of the Wilson gauge action parameter /3, with the corresponding values of the 
temperature, in units of Tc, indicated along the upper horizontal axis. Right panel: The renormalized Polyakov loop 
in the fundamental representation extracted from numerical simulations in the deconfined phase of SU(4) Yang- 
Mills theory (open symbols) , in comparison with the temperature dependence predicted perturbatively [8141 1815j 
at one- (solid line) or two- loop order (dashed line). Note that the Polyakov loop is normalized so that it tends to 
1 in the infinite-temperature limit, and that the leading-order perturbative correction for the renormalized loop is 
positive. Both plots are taken from ref. 



carrying k units of flux (to be denoted by a^), and found clear evidence for nearly perfect Casimir 
scaling: 

akOck{N-k). (154) 

While Casimir scaling is predicted by perturbation theory, at least at the lowest orders — see 
eq. (152) — , refs. |776[ I779j found that ratios of dual string tensions are consistent with Casimir 



scaling not only in the weakly coupling regime, but also down to surprisingly low temperatures, 
with only very small deviations close to T^. 

Interestingly, similar results have also been obtained in a recent study |784] , which investigated 
Polyakov loops in different representations and for different numbers of colors, from 2 to 6. This 
work found that the free energies associated with bare Polyakov loops in different representa- 
tions (of various A^-alities) are nearly exactly proportional to the eigenvalue of the corresponding 
quadratic Casimir, even down to temperatures close to while the renormalized Polyakov loops 
agree with perturbation theory |814[ 1815] at high temperatures, but receive large non-perturbative 
contributions at low temperatures (see fig. 17). Analogous results had been previously obtained 



for the SU(3) theory |349[ 1816] . In addition, ref. |784] also found that the value of the renormal- 
ized fundamental Polyakov loop for T — )■ (defined with a common normalization to 1 in the 
weak-coupling limit at high temperatures, and with a consistent choice of renormalization scheme 
for all groups) appears to be very close to 1/2, for all the gauge groups investigated. 

The equation of state (i.e., the temperature dependence of the free energy density) in SU(A^) 
Yang-Mills theory at large N has been independently investigated on the lattice by three different 
groups in a series of recent studies |777l I778[ 17801 1781] 1782] I783j . These works obtained consistent 
results, yielding a clear physical picture: the main finding is that, in the deconfined phase, the 
equilibrium thermodynamic quantities (pressure p, trace anomaly A, energy density e and entropy 
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Figure 18: Left panel: At temperatures close to the deconfinement transition, the quadratic dependence of the 
trace anomaly on T is exhibited by the linear dependence of the dimensionless ratio A/T'^ per gluon degree of 
freedom (with an appropriate overall normalization) on (from ref. |781j ). Right panel: The same effect — and 

the same consistency of results obtained in theories with a different number of color charges — can also be seen from 
the approximate temperature-independence of the quantity A/T^ per gluon degree of freedom (from ref. |783) ). 



density s — which are related to each other by elementary thermodynamic identities) per gluon 
degree of freedom have essentially the same dependence on T/T^ in all SU(A^ > 3) Yang-Mills 
theories. This means that, at least as far as it concerns equilibrium properties, the thermodynamic 
behavior of the "physical" theory with iV = 3 colors is the same as that of the theory in the 
large- limit. This result gives some confidence that the infinite- iV approximation underlying 
all holographic models of the QCD plasma |41Ut 1411^ 1412^ I413i 1414^ I415j should not be a major 
source of systematic uncertainty. In fact, a quantitative comparison between the lattice results 
for the large- equation of state and some holographic models was carried out in ref. |781] . 
considering, in particular, the equation of state constructed from the bottom-up model discussed 
in refs. [l32l iSl Ell Ei (see also refs. [mi [8201 [8211 [822] for related work) and the relation 
between the entropy density and the finite 't Hooft coupling for the = 4 supersymmetric 
Yang-Mills theory]^ as determined in ref. 



A further, interesting piece of information obtained in refs. |781l 1783] concerns the temperature 
dependence of the trace anomaly A. Up to temperatures of a few times (see fig. 18), this depen- 



dence seems to be quadratic — and possibly due to some contribution of non-perturbative origin,^ 
While this characteristic temperature dependence might be accommodated in phenomenological 
models ^M, [SSDl ISSD [H331 [SHU [S351 jHSni [S3Z], its physical origin at the fundamental level 
remains to be understood. At higher temperatures, the lattice results exhibit a smooth approach 
to the perturbative predictions. 

Other interesting physical quantities at finite temperature were investigated in ref. |775] : the 



^^As noted in ref. |823l . a similar type of comparison (but with lattice QCD results for = 3 and with dynamical 
quarks) was also carried out in ref. |824j . finding a similar value for the 't Hooft coupling at which the M — ^ 
theory best "mimics" QCDjA ~ 5.5. 

™As discussed in subsec. 2.5 the existence of non-perturbative effects in deconfined gauge theories is not surpris- 
ing. Such effects can be investigated numerically, via lattice simulations of dimensionally reduced effective theories: 
this was done, e.g., in refs. |8261 1827] for the SU(3) theory, and in ref. [828] for the generalization to the large-A^ 
limit. 
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Debye mass rriD, which characterizes the gluon screening in the deconfined phase, and which, at 
the leading order in a weak-couphng expansion, can be expressed as |321] : 



mD = ^^T, (155) 
and the string tension associated to large spatial Wilson loops, cr^, which, due to the strongly- 



coupled nature of the long-wavelength modes of the plasma (see the discussion in subsec. 2.5), 
is non-vanishing at any temperature. In particular, in ref. |775] the Debye mass was extracted 
from the connected correlator of Polyakov loops; the results showed the independence of mD 
from (except close to Tc), and also its (approximate) proportionality to the temperature in 
a temperature range 1.5Tc ^ T < 2.5Tc. The linear relation between mo and T is compatible 
with the perturbative prediction (assuming that the coupling is only mildly varying in the narrow 
temperature range considered), although there is no reason to believe that, at those temperatures, 
the coupling should be sufficiently small to validate a leading-order perturbative result. As for 
the spatial string tension, the simulations performed in ref. |775j indicated that, when is large, 
as is approximately constant for < T < whereas it is discontinuous at T = Tc (with 
lim^^^- (Ts{T) < lim^^^+ (Ts{T), by about 15%), and increasing with T in the deconfined phase. 

The temperature dependence of the topological susceptibility Xtopoi at large A^ was studied 
in refs. |768l 1838] . It was found that Xtopoi is nearly temperature-independent (and finite in the 
't Hooft limit) for < T < but is strongly suppressed (and, likely, vanishing for A^ — oo) in 
the deconfined phase. 

Conversely, the authors of the recent work |839] investigated the dependence of the critical 
temperature Tc on 6 in the SU(3), by monitoring the Polyakov loop susceptibility in simulations at 
imaginary values of 6 (to avoid the sign problem), and then performing an analytical continuation 
to real 6. They showed that Tc is decreasing when 6 varies from zero to a finite, real value, and, 
following ref. |283] , they obtained: 

at large A^, having combined their results with large-A^ values for the topological susceptibility, 
for the latent heat and for the critical temperature from refs. [32l I768[ 1775] . 

5.2. Results in three spacetime dimensions 

Laige-N gauge theories in three spacetime dimensions (3D) have been investigated in sev- 
eral lattice studies, with many findings qualitatively similar to those obtained in four spacetime 
dimensions (4D). 

SU(A^) Yang-Mills theories in 3D are characterized by a dimensionful gauge coupling: g'^ has 
the dimensions of an energy, so that a natural dimensionless parameter for organizing pertur- 
bative expansions for physical processes at the momentum scale k is the g^/k ratio (see, e.g., 
refs. |840l 1326] and references therein). They are asymptotically free at high energy, and super- 
renormalizable; in addition, one can also prove that the Coulomb potential in 3D is characterized 
by a logarithmic dependence on the distance — a behavior which, in particular, would already be 
sufficient to imply color confinement. However, perturbation theory is not reliable in the limit of 
large distances, in which these theories become strongly coupled. The computational tool to study 
them in the non-perturbative regime is, thus, lattice simulations — in which the (3 parameter ap- 
pearing as a coefficient in front of the Wilson gauge action is defined as (3 = 2N/ {ag^) = 2N'^/ (aAo), 
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where Aq denotes the bare (lattice) 't Hooft couphng. Like in 4D, SU(A^) lattice gauge theories in 
3D are characterized by a strong coupling regime at small /3, and by a weak-coupling regime at 
large /3; in 3D, the two are separated by a crossover, which seems to turn into a third-order phase 
transition for N ^ oo |841j . Recently, a similar study has also been performed for the phase 
structure of 3D S0(2A^) Yang-Mills theories 



The numerical study of large-A^ Yang-Mills theories in 3D was initiated during the 1990's, 
and the main results are summarized in refs. |842l 18431 1844] : in particular, at zero (or low) 
temperature, all 3D SU(A^ > 2) gauge theories are confining, with a non-perturbatively generated, 
asymptotically linear quark-antiquark potential V{r) at large distance r. In particular, the string 
tension a (i.e., the asymptotic slope of the potential at large r) is related to the bare 't Hooft 
couphng Ao by |844] : 

niQ7^™ 0-1200(29) 
— = 0.19755(34) — . (157) 

More recently, the numerical determination of the string tension in 3D Yang-Mills theories was 
pushed to very high precision in ref . [845J : the goal was to test the predictions formulated in a series 
of works by Karabali, Kim and Nair |84^ 1847] I848[ 1849] . in which the vacuum wavefunction and the 
string tension of 3D Yang-Mills theories were computed analytically in a Hamiltonian approach. 
The results of previous lattice studies |843l 1844] had shown some discrepancies with respect to 
these predictions, but could not give a conclusive answer about their validity for N ^ oo. On the 
contrary, the authors of ref. |845j presented an accurate extrapolation to the large- limit, which 
revealed a discrepancy of six standard deviations with the predictions of ref. |848] . The interesting 
analytical program discussed in refs. |846[ 18471 18481 1849] is presently ongoing, and refinements of 
the string tension predictions have been recently proposed |850] ; however, it is fair to say that its 
theoretical foundations are still under debate |851] . 

Inspection of the numerical results for V{r) at intermediate distances reveals that the potential 
also includes a 1/r contribution, with a coefficient that is compatible with the Liischer term 
prediction: this means that, similarly to what happens in 4D, also in 3D the low-energy dynamics 
of confining fiux tubes can be described in terms of (massless) fluctuations of a bosonic string |710l 
I7T51 IM31 18521 18531 [85il 18551 18561 1857] . Again, this is not specific of the large-iV limit only, 
but rather appears as a generic feature of all confining models in 3D: the Liischer term (or the 
broadening of confining fiux tubes with their length |858l 1859] I86U] 1861] 1862] 1863] 1864] , which is a 
related imphcation of the same effective picture) is also observed in high-precision studies of SU(2) 
Yang-Mills theory [H65l MEl EHTl EMI EHHl ElQl EST], of 3D lattice models with local invariance 
under a discrete gauge group [8571 EnZl EUl MM MM MM MM MM MM MM MM EHQl EHI] and 
even in random percolation models (with an appropriate, topological definition of Wilson loop 
operators) 



In fact, these studies are now reaching a level of numerical precision, at which it becomes 
possible to investigate subleading (in 1/r) corrections to the confining potential, and possibly 
compare them with analytical computations that predict at which order the correct effective string 



model should deviate from the Nambu-Goto string |742] 1743] 1744] I745"l 1747] . As an example, fig. 19 



taken from ref. [856], shows the torelon spectrum obtained from simulations in SU(6) Yang-Mills 
theory, in comparison with the closed-string energy levels predicted by the Nambu-Goto model, see 
eq. ( 141 ). As shown by this figure, in which both the torelon length and the energies are measured 
in appropriate units of the relevant energy scale (set by the square root of the string tension), a 
somewhat surprising feature is that the agreement with the Nambu-Goto model appears to hold 
down to very short distances, at which a priori there is no reason to expect that the approximation 
of the confining fiux tube as a uni-dimensional string should still be valid. In fact, the deviations 
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Figure 19: Ground state and excitation spectrum of a closed confining flux tube, winding around a spatial size of 
the lattice (torelon) in SU(6) Yang-Mills theory in 3D, from ref. |856| . The symbols show the lattice results for the 
energies of different states, as a function of the torelon length, while the lines denote the corresponding predictions 



from the Nambu-Goto string model, i.e., eq. (141) for 13 = 3. In particular, the solid lines refer to the predictions 



using the continuum definition of the momentum, while the dashed lines are obtained from the lattice momentum 
(in a free theory) — which differs from the former at scales of the order of the inverse lattice spacing. 



from the curves predicted by the Nambu-Goto model are mostly accounted for, by replacing the 
continuum definition of the momentum with its lattice counterpart (the two become significantly 
different only at scales comparable with the inverse lattice spacing). 

The spectrum of 3D Yang-Mills theories at large (including several excitation states) was 
studied via lattice simulations in various works |753[ 18421 18431 18441 1884] (and compared with a 
variant of the Isgur-Paton model |885j in ref. |886j ): the physical states are hadrons (glueballs), 
which can be classified according to the irreducible representations of the SO (2) group in the 
continuum (or to the corresponding subduced representations restricted to the square group, if 
the theory is regularized on a cubic Euclidean lattice), as well as a "mirror" parit}''^^ and charge 
conjugation. The simulations show that these theories have a finite mass gap (which remains 
non- vanishing in the large- iV limit), and that ratios of different masses are almost independent of 
— up to small, 0(1/A^^) corrections. 

A lattice study of correlation functions of gauge-dependent quantities (gluon and ghost prop- 
agators) in the Landau gauge in 3D Yang-Mills theories was carried out in ref. |887j . This work 

^^Note that, for any even number d of spatial dimensions, the "usual" parity transformation, corresponding to 
inversion of all spatial coordinates, is nothing but a rotation, because — 1 G SO((i) for d even. 
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compared the results obtained in SU(A^) theories with a number of colors ranging from 2 to 6 
and in G2 Yang-Mills theory, finding essentially no dependence on the rank of the gauge group 



72 



The properties of 3D SU(A^) Yang- Mills theories at finite temperature are qualitatively similar 
to those found in 4D: there exists a deconfinement transition at a critical temperature T^, at which 
color-singlet hadronic states give way to a plasma of deconfined particles, and the value of (in 
units of the square root of the string tension) remains finite in the 't Hooft limit 



Tc , , 0.880(43) , , 

= 0.9026(23) + (158) 



The transition is of second order for "small" values of A^, i.e. for SU(2) and SU(3) gauge theory; 
for SU(4), identifying the order of the transition has been a challenging problem |893j . but the 
most recent and accurate studies indicate that it is a weakly first-order one |894j . like for > 
5 |892[ 1895] — in agreement with the intuitive picture of a more and more abrupt transition, due 
to an increasing imbalance between the free energy in the confined phase (which is 0(1)) and the 
same quantity in the deconfined phase, in which it scales like 0{N'^). 



The large- iV lattice computations performed in refs. |896l 1897] have shown that in the confining 
phase the equation of state can be adequately modelled in terms of a gas of almost non-interacting, 
massive glueballs, as obtained by summing the known glueball states |843i I844j — except very close 
to Tc, in which heavier states (whose density can be modelled in terms of a simple bosonic string 
model, following an approach similar to the one discussed in ref. |886j ) appear to become relevant. 
Except in a narrow region of temperatures close to Tc, the equation of state is approximately 
independent of A^. On the contrary, in the deconfined phase the pressure (like the other equilibrium 
thermodynamics properties) is nearly perfectly proportional to the number of gluon degrees of 
freedom (A^^ — 1), although the approach to the asymptotic value given by the Stefan-Boltzmann 
limit: 

(where ({3) denotes Apery's constant, approximately equal to 1.20205690316 . . . ) is rather slow: 
at temperatures around 7 T^, the pressure is still about 15% smaller than the limiting value in 



eq. (159). This indicates that, at those temperatures, the deconfined plasma is still far from a 



gas of non-interacting gluons (for the SU(3) theory, analogous results have also been obtained 



in ref. |898j ). Fig. 20, taken from ref. [897], shows a plot of the trace of the energy- momentum 
tensor A in units of (which equals the derivative of p/T^ with respect to InT) per gluon d.o.f., 
as a function of T/Tc. the perfect agreement between the simulation results for different gauge 
groups shows that, in the deconfined phase, the equilibrium thermodynamic quantities are exactly 
proportional to A^^ in the large- A^ limit. An interesting feature emerging from these results is that, 
like in 4D, also in 3D A is proportional to in this range of temperatures. Finally, the solid 
yellow curve in the figure shows the prediction from a generalization of the holographic model 
proposed in refs. |817l 1818] for 3D Yang-Mills theory. 

5.3. Results in two spacetime dimensions 

QCD toy models in two spacetime dimensions (2D) are interesting theoretical laboratories, 
in which, by virtue of a limited number of physical degrees of freedom, various problems can 



Analogous results have also been obtained from the comparison of the SU(2) and SU(3) theories in four 
spacetime dimensions [5551 [HMl 155111 IMT] . 



77 



Trace of the energy-momentum tensor 



0.18 



0.16 
0.14 
0.12 - 

0.08 
0.06 
0.04 
0.02 




I I I 1 1 I 1 1 I I I I I I I I 1 1 I I I I I I I I 1 1 I 1 1 I I I I I I I I 1 1 I 1 1 I I I I I 1 1 I 1 1 I I I I I I I I 1 1 I 1 1 I I I I I 1 1 I 1 1 I I 



I I I I I 1 1 I 1 1 I 



I I I I I I I I 1 1 I 1 1 I I I I I 1 1 I 



SU(2) 
. SU(3) 

• SU(4) 

• SU(5) 
. SU(6) 

— holographic model, for a = 3 / 2 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



0.5 1 1.5 



2 2.5 3 3.5 4 4.5 
TIT 



5.5 



6.5 7 7.5 8 



Figure 20: Temperature dependence of the trace of the energy-momentum tensor in units of per gluon from 
simulations of the deconfined phase of 3D SU(A^) Yang-Mills theories, from ref. |897j . The plot shows the simulation 
results for = 2 (brown symbols), 3 (black), 4 (green), 5 (blue) and 6 (magenta), and their comparison with the 
expectation (yellow curve) obtained from a generalization of the holographic model proposed in refs. [8171 ISlBj to 
3D Yang-Mills theory. 



be studied analytically (or semi-analytically) in the large- iV limit. As we discussed in sec. |2| 
classical examples include the determination of the meson spectrum in the 't Hooft model p], 
the discovery of a phase transition separating the weak- and strong-coupling phases of the lattice 
regularization of the theory with the Wilson gauge action |56l and the analysis of the 



spectral density of Wilson loops in the continuum [59j (which continues to attract attention to 

this dav pomroTOTOMimiTiMroigTU] ). 

Recent examples of analytical results in 2D models of QCD at large include, e.g., those 
obtained in the study of form factors in the principal chiral model [91H I912[ 1913] — a theory which 
has also been studied numerically |914j (see also ref. j915] ). 

Finally, another problem that has been recently investigated in the context of 2D large- iV QCD, 
is the one of baryonic matter at finite density |916l I917[ I918j . In particular, in ref. |917j it was 
found that, at finite baryon density, the ground state of the system breaks translation invariance, 
with the formation of a chiral crystal; related topics have been discussed in refs. |919l 1920] . 



^•^Note, however, that this phase transition is not physical, and its nature and very existence depend on the 
lattice discretization that is used |899|. 19001 1901] — see also ref. [902] for a discussion. 
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6. Conclusions 

The idea put forward by 't Hooft in liis seminal article [Ij proved to be among the most 
fruitful ones of the last few decades in theoretical elementary particle physics, and paved the 
way to a number of developments in different directions. The implications of the large- limit 
for non-Abelian gauge theories have been studied in fields ranging from string theory to QCD 
phenomenology, and in many cases have even generated whole new research areas. Along with this 
rich variety of developments, however, also came an increasing specialization and fragmentation of 
the different subfields — a trend presently common to all scientific disciplines (or, more generally, 
in nearly all branches of human knowledge). In our view, this excessive specialization (together 
with a tendency towards increasing precarity in the scientific job market) may have harmful 
consequences for the long-term progress of science, as it naturally limits the possibility of fruitful 
exchanges of ideas between scholars working in different, albeit closely related, areas. 

This state of affairs gave us our first motivation to write the present review, whose primary 
purpose consists in offering a broad, unitary, panoramic view of large- A^ gauge theories. While 
describing the general aspects of this topic in section |2| we tried to emphasize the main ideas and 
the fundamental physical concepts, leaving a technical account of the details out of our discussion. 
We tried to be as pedagogical as possible, precisely in order to facilitate the comprehension of 
the various problems, studied by different research communities, as clearly as possible. In doing 
so, on several occasions we were, unfortunately, bound to be incomplete, or not as accurate as 
we wished. There are many topics, which we could only touch upon, or mention very briefly: 
in apologizing for our omissions, we hope that, in this review, the interested readers can find at 
least a reference to an article or to a review presenting a more detailed discussion of some specific 
subject — but, at the same time, we also encourage them to read the sections that may look not 
directly relevant for the topic they are interested in. 

Our second motivation for writing this manuscript was, admittedly, more directly related to 
our own research work, and consisted in presenting the current status of lattice studies of large- A^ 
gauge theories. These studies (both the analytical and the numerical ones) are now reaching a 
level of maturity, and can give precise answers to many theoretical questions, which are relevant 
in a large variety of contexts. We hope that this manuscript serve as a vehicle, to inform the large 
community of researchers interested in topics related to 't Hooft 's limit about the recent progress 
in the lattice field. We found it important to devote an entire section (sec. [s]) to a pedagogical 
introduction of the lattice formulation of field theory, in order to facilitate the interpretation of the 
results reviewed in sec. [5] also for non-experts. In doing so, we hope that we succeeded in clarifying 
the main virtues and strengths of lattice computations, and in clearing out some misconceptions 
or prejudices which are, unfortunately, not so uncommon among non-practitioners. 

The overview of lattice results presented in sec. [5] shows the convincing numerical results that 
have been obtained for large-A^ gauge theories during the last fifteen years, and highlights the 
progress in areas which, until recently, appeared to be out of the reach of lattice calculations. 
Highly excited states of confining flux tubes, hadron spectroscopy, thermal observables (to name 
but a few) are all quantities for which lattice computations can now provide very precise and 
accurate results — while the first few dynamical simulations at A^ > 3 are beginning to appear. 
On the theoretical side, the last decade has also witnessed very significant developments on the 
topics related to the time-honored idea of large- A^ volume independence, for which, in particular, 
there has been a lot of progress in understanding the transition from the strong- to weak-coupling 
regimes. At the same time, mathematical tools inspired by orbifold equivalences in string theory 
led to deep conceptual advances in our understanding of large- A^ volume independence and large- A^ 
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equivalences between theories with different field content. As we discussed in sec.|4| the techniques 
that have been developed in the process are also finding practical applications for some of the most 
challenging problems to be studied on the lattice — supersymmetry is one, might QCD at finite 
density be the next? 

To summarize the state-of-the-art in this field, we can say that, while many problems have 
already been successfully studied, many others are currently still open, and being actively inves- 
tigated. In view of the advances in lattice QCD computations in the last few years, we foresee 
further, dramatic progress for large- studies on the lattice in the near future. For these reasons, 
we would like to conclude this review encouraging especially the young, ambitious and talented 
researchers to join this field. 
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